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XXX. BAND DRITTES HEFT 1961 


Der dynamische Kardanfehler beim Vermessungskompab 


Von M. Klumpp 


1. Einleitung. Im Anschlu8 an verschiedene Untersuchungen von Magnus soll in der folgenden 
Arbeit die Auswirkung von Unwuchten in der Kreiselaufhangung beim Vermessungskreiselkompa8 
gezeigt werden. 

Magnus’? hat u.a. Auswanderungserscheinungen an schwingenden, kardanisch gelagerten, 
kraftefreien Kreiseln nachgewiesen, die durch die Deviationsmomente der inneren Kardanringe 
hervorgerufen werden. Hier soll nun ein von K. H. Stier? beschriebener Vermessungskompa8B 

untersucht werden, an dessen Anzeigegenauigkeit groBe Anforderungen gestellt werden. Ziel der 
Berechnungen ist, den mittleren azimutalen Anzeigefehler zu bestimmen und anzugeben, wie dieser 
von der Massenverteilung der Kreiselaufhangung abhangt. 

Bei dem VermessungskompaB besteht das nordweisende System aus einer schwimmenden Kugel, 
in der ein Rotor mit horizontaler, nordweisender Drehimpulsachse lauft. Der Schwerpunkt des 
Gesamtsystems Kugel samt Rotor liegt um die Metazenterhéhe h unter dem Auftriebsmittelpunkt. 
Da dieser Kompa8s nur ortsfest (d. h. erdfest) verwendet wird, geniigt ein Kreisel zur fehlerfreien 
Nordanzeige. 


2. Voraussetzungen der Rechnung, Bezeichnungen. In Abb. | sind Rotor und Kugel schematisch 
dargestéellt. Folgende Koordinatensysteme werden eingefiihrt: 

1. Ein erdfestes (x, y, z)-System 2, dessen Ursprung in den Schwerpunkt des Rotors gelegt 

wird, der mit dem Auftriebsmittelpunkt tibereinstimmen soll. In der Ruhelage (« = 6 = y = 0) 


Abb. 1. Schema des Vermessungskompasses, 


5 sail die y-Achse mit der nordweisenden Drehimpulsachse des Rotors zusammenfallen und die 
- lotrechte z-Achse durch den Schwerpunkt des Gesamtsystems gehen. 
*s 1K, Magnus, Advances in Aeronautical Sciences, Vol. I, S. 507, London 1959. 


2 K. Magnus, Z. angew. Math. Mech. 35 (1955) S. 23. 
~ Ku ene area d. Markscheidewesen 60 (1953) Heft 1. 
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2. Ein kugelfestes (xx, ¥x, %x)-System Xx, das fiir 1 =f =y =0 mit 2 tibereinstimmt. 
In Xx hat der Schwerpunkt des Gesamtsystems also die Koordinaten (0, 0, —h). 

3. Ein rotorfestes (xp, Vp, 2g)-System Xp, das fiir yp = 0 mit dix ubereinstimmt. 

Dabei sind «, B und y die Winkel fiir die Drehungen der Kugel um die Achsen xx, yx und 2x, 
wihrend sich der Rotor um die yg-Achse relativ zu 2x, umy dreht. In Abb. | sind die entsprechen- 
den Drehvektoren eingetragen. 

Die Winkel .«, B, y legen die Lage der Kugel im erdfesten System nicht eindeutig fest. Da aber 
die Winkel « und f schon aus konstruktiven Griinden sehr klein sind, werden einige Vereinfachungen 
méglich, die die Rechnung wesentlich erleichtern. Vor allem wird die azimutale F ehlanzeige als 
Mittelwert (y) der Ost-West-Schwingungen der Drehimpulsachse berechnet ohne Beriicksichtigung 
von Kippung (8) und Elevation («). Das stimmt mit der Ablesung iiberein, bei der die Rotor- 
achse yx und nicht deren Projektion in die Horizontalebene beobachtet wird. 

Weiterhin wird vorausgesetzt, daB simtliche Bewegungen reibungsfrei erfolgen. 


3. Ermittlung der auBeren Momente. Als auBere Momente wirken auf die Kreiselkugel das 
Schweremoment und das Fiihrungsmoment infolge der Erddrehung. 

Das Schweremoment entsteht durch die Schwerpunktverschiebung um die Metazenterhohe h 
und ist nach GréBe und Richtung abhangig von der Lage der Kugel (oder der Lage von 2x). In 
jedem Fall ist die Wirkungslinie des riickdrehenden Schweremomentes die Schnittgerade von 
(x, y)-Ebene und (xx, ¥x)-Ebene, da der Momentenvektor senkrecht auf der (z, zz)-Ebene stehen 
mu. 

Fiir kleine Winkel « und 6 — wenn also die trigonometrischen Funktionen bez. 4 und f lineari- 
siert werden kénnen — ergeben sich dann fiir das Schweremoment folgende Komponenten: 


um die x,-Achse My, —Qa, 
um die yx-Achse Mz~ —Q6, 
mit Q = Gh, wobei G das Gewicht von Rotor samt Kugel und h die Metazenterhéhe bedeuten. 
Ausgehend von denselben Voraussetzungen la8t sich das Fiihrungsmoment infolge der Erd- 
drehung, das auf den Kreisel wirkt, einfach ermitteln. Seine Komponenten sind bekannt!. Sie 
betragen 
in Richtung —a: K, = D,w* 4 cos cosy, 
in Richtung —y: K, = D,w* cosgsiny , 
in Richtung + 4: K; = D,w* sing, 
wobei D, = Bp w den Betrag des Drehimpulses des Kreisels mit der Drehmasse Bp um die yx- 
Achse und der Winkelgeschwindigkeit w bedeuten, wihrend w* die Winkelgeschwindigkeit der 
Erde und 9 die geographische Breite des MeBortes sind. 
Werden die trigonometrischen Funktionen noch beziiglich des Winkels y linearisiert (auch die 
Amplitude der Azimutschwingung wird klein sein) so hat man einfach 


K, = D,w* «cosy = Ko, 
K, = D,w* ycosp=Ky, } mit K = D,w* cos. 


4, Aufstellen der Bewegungsgleichungen. Die Drehmassen werden mit A, B und C, die Kipp- 
massen mit D, EF und F' bezeichnet. Dann haben die Trigheitstensoren folgende Gestalt: 


fiir die Kugel fiir den Rotor 
A —F —E |Ap 0 0 | 
—F B —D|=60,, 0 Bp 0 Pa 
—E —D C 0 0 AR 


Wegen der auSerordentlich sorgfaltigen Gestaltung und Auswuchtung der Rotoren fiir Kreisel- 
kompasse soll hier mit einem symmetrischen Rotor und verschwindenden Kippmassen gerechnet 
werden. AuSerdem ist die Hauptmasse des Rotors elastisch mit der Welle verbunden, wodurch 
sie sich bei iiberkritischem Lauf so einstellt, da8 die Drehung um eine Hauptachse erfolgt. 

Der Drehvektor u der Kugel hat voraussetzungsgemaB die Komponenten 


u = (a; 8,7), (1) 
1 Vgl. z. B. R. Grammel, Der Kreisel, 2. Aufl. Bd. 2, S. 94, Berlin-Gottingen- Heidelberg 1950. 
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der des Rotors ergibt sich einfach zu 


Un=Treut+y (2) 
wobei Tax die Matrix der Koordinatentransformation von 2p nach 2’, ist. In Komponenten 
geschrieben hat man 


U.~R =ACosp—ysiny, 
gescapa gt 2’) 
U,R =ASin yp + cosy. 

Die Drehbeschleunigungen sind dann 


Wea, UB i ey, (3) 
Upp =X CoS p—apsiny—ysiny—ypeosy, | 
typ B+, © 


tp = sin yp +48 cos p + ycosy—yypsiny. 
Die Bewegungsgleichungen werden iiber die Komponentendarstellung der Eulergleichung auf- 
gestellt, die mit dem Drehimpulsvektor ), dem Drehvektor u und dem Momentenvektor 9 lautet 
d'D 


wobei - die zeitliche Ableitung im kérperfesten System bedeutet. In Komponenten geschrieben 


erhalt man 
fur den Rotor : 
ARU,p — (Br — Ap) yn aR = Myr: 
Briyr = Myr ’ (5) 
Ap l,p — (Ap — Ba) Ugr UR MR; 
fiir die Kugel 
A&é—F§—E}—(B—C)B) +7 (F4 + D7) —B (Ed + DB) =M,, 
—F&é+BR—Dy—(C—A)ay +4 (Ea + DB)—jy (FB + Ey) =M,, (6) 
—E&—Dp+Cy—(4—B)ap +2 (FB + £9) —4 (Fa + Dy) =M,. 

Neben den berechneten duReren Momenten sind in diese Gleichungen die inneren Momente 
zwischen Rotor und Kugel einzusetzen. Ist Jt’ das Moment, das der Rotor auf die Kugel ausiibt, 
gemessen in x, so gelten folgende Beziehungen: 

Mr=TrxM und Nt=—M' +M,, 
wo Nt, die iuBeren Momente sind, die an der Kugel angreifend gedacht werden. Im einzelnen ergibt 
sich : 
M,.x = M, cos y — M, siny, 
Myr = M! 3 
M,p = M,siny + M, cosy, 
M, =— M;— K,. 


(7) 


Fiir die weitere Rechnung ist es sinnvoll, das Antriebsmoment M,, des Rotors gleich Null zu 
setzen, da es wahrend der kurzen Zeit einer Periode der hier zu betrachtenden schnellen Schwin- 
gungen ohne Bedeutung ist?. Mit (2’) und (5) gilt somit 


B+ p=0 
oder : 
6 + yp =konst. = w (9) 
und mit (8) M, melon: (10) 


1 Siehe z. B. FuBnote 2 von Seite 153. 
i 
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AuBerdem lassen sich durch Kombination der Gleichungen (7) und (8) zwei Beziehungen finden: 


M,=— up sin y— M,pc0sy—Q4+ K,—K,, (11) 


x 


M, = M, sin y — M,z 008 y— K,. (12) 
Durch Einsetzen von (1) bis (9) erhalt man aus (10) bis (12) die Bewegungsgleichungen 
BjB— Fa—Dj +06 =fay+EP+ Fb —Dip—E®@, (13) 
ai —FB—Ep—by+(Q+K)4—K, =—dp y+ Eap—Fip+De—Dy, (14) 
c}—DB—E& +b +Ky =ch$ + Fo?— Ff? + Day—EBy, (15) 
mit den Abkiirzungen 
a= A + Apz $ 
b = Bao, 
c= C+ Ap, 


d=C—B-+A,g, 
e = Ar + A a B 5 
f=Cc—A. 
Aus Gleichung (14) kann sofort die ,,stationdre Elevation berechnet werden. Sie betragt 
_ _Ks 
Sete e 
Wenn weiterhin unter « die Schwingungsausschlage um die Nullage a) verstanden werden, kann K, 
in (14) gestrichen werden. 


a 


5. Lésung der Bewegungsgleichungen. Die linearen Glieder der Differentialgleichungen (13) 
bis (15) stehen auf ihren linken Seiten. Bei einer Linearisierung beziiglich der Variablen o, 8 und y 
sind also die rechten Seiten gleich Null zu setzen. Die Lésung des verbliebenen linearen Teils kann 
als erste Naherung aufgefaBt werden. Diese erste Naherung ist von der Form 


v1 = Yasin (nt), 
O, = 4,sin (nt -+ é), (16) 
B, = 8, sin (nt + @)- 


Das Ziel der Untersuchungen ist, durch Bildung eines Mittelwertes fiir die rechten Seiten der 
Differentialgleichungen (14) und (15) tiber eine Periode der Naherungslésung (16) den EinfluB der 
nichtlinearen Glieder abzuschatzen um so die azimutale Fehlanzeige und die ,,dynamische Ele- 
vation®* zu ermitteln. 

Die Frequenz n in (16) kann aus der charakteristischen Gleichung des linearisierten Systems (13), 
(14), (15) bestimmt werden, die die Gestalt hat 


18 [Bac— F?c—2 DE F— BE*—aD®] + 
A! [b? B + Q (ac — E®) + (Q + K) (c B— D®) + K (a B— F2)] + 
PQ + K}(KB+ Qc) +Q(h +aK)]}+Q0(Q+K)K=0. (17) 


Diese Gleichung dritten Grades in /? 1aBt sich allgemein nicht explizit lésen, auch dann nicht, wenn 
durch verniinftige GréBenbetrachtungen einige Vereinfachungen vorgenommen werden. Natiirlich 
kann (17) numerisch mit beliebiger Genauigkeit ausgewertet werden, doch lassen sich auf anderem 
Wege itibersichtlichere Naherungslésungen finden. Die charakteristische Gleichung (17) wird lésbar, 
wenn der Kugel entweder ein Freihcitsgrad genommen wird, also z. B. B =~ =f =0 gesetzt 
wird, oder wenn die 4uBeren Momente verschwinden. Dabei geniigt es schon, entweder 0.=% 
oder K = 0 anzunehmen. 


4 


Die geringste Anderung diirfte eintreten, wenn der Kugel der Freiheitsgrad um die yx-Achse 


genommen wird (6 = i= Se ies = 0), da Winkelanderungen & und 3 y keine Kreiselmomente um die 
¥x-Achse hervorrufen. Unter dieser Voraussetzung erhalt man als Frequenzen aus der charak- 
teristischen Gleichung 


ny = |/P+eK+QtK)c+ Vie + aK +(Q+K)cP—4@+ K) (ac— EF 
2 (ac — E?) : (18) 
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Da gegeniiber b? alle anderen Glieder klein sind und sicher E? < ac ist, ergibt sich die Naherungs- 


frequenz der schnellen Nutationsschwingung, die fur das hier interessierende Ergebnis vor allem 
wichtig ist, zu 


b 
ears (19) 


Das ist die bekannte Naherungsformel fiir die Nutation eines schnellen Kreisels. 

Werden die auBeren Momente gleich Null gesetzt (Q = K = 0), so laBt sich die schnelle Nuta- 
tionsfrequenz durch verniinftige Vernachlassigungen ebenfalls auf (19) zuriickfiihren. Durch dieses 
Ergebnis wird plausibel gemacht, da8 die aus (17) folgende Frequenz nicht viel von der Naherungs- 
lésung (19) abweichen wird. 

Wird unter der vereinfachenden Voraussetzung (Q = K = 0) der EinfluB® der nichtlinearen 
rechten Seiten von (14) und (15) auf die Lésung dieser Gleichungen abgeschatzt, so ergeben sich 
die schon von Magnus berechneten mittleren Auswanderungsgeschwindigkeiten der Kreiselachse. 
Erst die 4uBeren Momente begrenzen die Auswanderungen und geben Gleichgewichtslagen, aus 
deren Verschiebung die gesuchten Fehlanzeigen ermittelt werden. 

Zur Bestimmung von Phasenwinkel und Amplitudenverhiltnis fiir die Lésung der linearen Teile 
von (13) bis (15) wird der Lésungsansatz (16) besser etwas anders geschrieben: 


V1 = 7a 8in (nt), ' 
a, = Msin (nt) + Neos (nt), (20) 
By, = Ssin (nt) + T cos (nt). 


Dabei wird der Zusammenhang mit (16) durch folgende Beziehungen hergestellt: 
y= tee. S=ty, /N2+ M@=a,, VT? + S?= By. 


Geht man mit (20) in die linken Seiten von (13) bis (15) ein, so erhalt man Gleichungen, die die 
gesuchten GréBen zu berechnen gestatten. Man findet 


oe (21) 
"(ae—o +?) 
one . bn (22) 
Er? + =; [an*—(Q + K)] 
FDnr*® \2 
Set) E+ gn—o) (23) 
— Die aree?, 
Ve (e+ +"(gR 94) 
iene + [emt 42 [a n?—(Q +x)]| 
Ba = fl 21 (Bn? ») ae (24) 
v4, [Frnt al n?—Q | 


Diese Ausdriicke miissen ebenso wie (17) der Zahlenrechnung zugrundegelegt werden. 
Die gesuchten Fehlanzeigen ergeben sich als Einfliisse der nichtlinearen Glieder der rechten 
Seiten in den Differentialgleichungen (14) und (15): 


R,=—dpy +Eap—Fay+Dp24+D/, 
R, =edB+Day—Epy+ Fe—Fp. 
Mit der Naherungslisung (16) werden nun die iiber eine Periode genommenen Mittelwerte gebildet 
20 
R, = - Ry (4%, Bry V1) 4(" t) , 
6 


Ry = aa | Rr (arPrn)d (m4), 
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womit sich dann die Mittellagen der Schwingungen naherungsweise bestimmen lassen: 
es 1 = 5 tetas 
a= ge y= RR,- 


Man berechnet die folgenden Ausdricke: 


= YE Ba 4B a 2 iat F%A Dao, 25 
oa —d“cosy+E D (cos € cos g + sin é€ sin ¢) oe cos € + vy (25) 
2 
Sela eA (cos ¢ cos p + sine sin 9) + D?4 cose ppp e cosg + F— FL (26) 
A A A 


2K OF, 
Hier sicht man am besten die Bedeutung der schnellen Nutationsschwingungen fiir die Fehlanzeige, 
da sowohl & wie y den Faktor n? haben. Diese Abhangigkeit wird auch nicht aufgehoben durch die 


Gleichungen (21) bis (24), die ja ebenfalls von n abhangig sind. 


6. Auswertung der Ergebnisse. Eine kurze numerische Auswertung soll nun Auskunft tiber die 
praktische Bedeutung der Anzeigefehler geben und helfen, mittels nachgepriifter Vernachlassi- 
gungen einfache Endergebnisse zu formulieren. Ihr werden die folgenden Werte zugrundegelegt, 
die entweder vom Institut fiir Markscheidewesen in Clausthal zur Verfiigung gestellt wurden, oder 
aber auf Schatzungen beruhen, die nach Méglichkeit abgesichert sind: 


Ap = 34kg cm’, A = 240 kg cm?, 
Br = 43 kg cm’, B = 280 kg cm?, 
a~= 2030-s+ , C = 300 kg cm?, 
G = 8,3410? kg cms, O29 10-8. 
he = 7575 cn = 


Die Kippmassen werden vorerst willkirlich mit D = EK = F = 10 kg cm* angenommen, was 
etwa 3% von Centspricht. Diese Werte gelten fiir die Bestimmung der Nutationsfrequenz und fiir 
weitere Abschatzungen. Es wird sich zeigen, daB die geforderte Anzeigegenauigkeit des Kompasses 
eine Beschrankung der Kippmasse F verlangt. 

Die schnelle Nutationsfrequenz ergibt sich aus (17) zu n = 288s-1, wahrend die Naherungs- 
formel (19) n ~ 289s! ergibt. Gleichung (19) ist also genau genug, was auch aus (17) direkt 
erkannt werden kann. 

Werden nun die Ergebnisse (21) bis (24) der exakten Liésung der linearisierten Differential- 
gleichungen (13) bis (15) weiter untersucht, so ergeben sich folgende Abschatzungen, die eine fiir 
praktisch vorkommende Falle ausreichende Genauigkeit haben diirften: 


ee ee (21’) 
tep=—E ye, (22') 
eae ee “ 
oa a) 


Diese Gleichungen entsprechen genau denjenigen, die bei einer Vereinfachung durch Vernach- 
lassigung der duBeren Momente gefunden werden; die Werte (21’) und (23’) stimmen auBerdem 
mit den entsprechenden Ergebnissen fiir 8 = 6 = # = 0 iiberein, so da die friiher durchgefiihrten 
Naherungsbetrachtungen nachtraglich gerechtfertigt erscheinen. 

Die Anzeigefehler (25) und (26) miissen noch genauer untersucht werden, da sie das eigentliche 
Ergebnis dieser Arbeit darstellen. Die mittlere »,dynamische Elevation‘ & wird so klein, da ihr 
keine weitere Beachtung geschenkt zu werden braucht. Das verhaltnismaBig groBe Riickstell- 


moment (Q + K) ist die einleuchtende Ursache hierfiir. Fiir die mittlere azimutale Fehlanzeige 
findet man 


= Wad, 1 ly/a n? 02, B? x2 , 
ee re [PE(2— 5/2) + Fle We Sey ere 26) 
pw aE 2 
Se Rega * (27) 
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Von den Kippmassen hat nur F einen wesentlichen Einflu8 auf die GréBe von y. Die Forderung 
nach einer Anzeigegenauigkeit von 1 Bogenminute fiihrt nach (27) zu der Forderung, den Betrag 
von F auf 0,55 kg cm? oder 0,2°% von B zu begrenzen, wobei die Amplitude der Nutationsschwin- 
gung mit y, = 2- 10-4, also erheblich weniger als 1 Bogenminute, sowie y = 50° angenommen 
wurde. 

Die Formel (27) fiir die azimutale Fehlanzeige 1a8t sich leichter itibersehen, wenn anstelle der 
Abkiirzungen wieder die urspriinglichen Werte eingesetzt werden. Es folgt dann 


eo BroF 2 
v= oF (A + Ap)? cos p A> G2) 


Hieraus ergibt sich: Wird der Drehimpuls Bz w des Kreisels als konstant angenommen, so hangt 
die Fehlanzeige y auBer von der Kippmasse F der Kugel noch wesentlich von der Amplitude der 
Nutationsschwingung y,, von der Drehmasse (A + Ag) des Gesamtsystems um die xx-Achse und 
von der geographischen Breite y ab. Fiir y + 90° wachst y iiber alle Grenzen, da das Riickstell- 
moment zu Null wird; jedoch ist der Kreiselkompa8 an den Polen sowieso unwirksam. 

Wie an dem berechneten Beispiel deutlich wird, muB bei hochgenauen Kreiselkompassen auch 
die Kreiselaufhangung sorgfaltig ausgewuchtet werden. Da auch bei anderen Kreiselgeraten die 
Genauigkeitsanforderungen laufend verscharft werden, wird man auch bei ihnen der Verteilung 
aller bewegten Massen gréBere Aufmerksamkeit zuwenden miissen als bisher. 


(Eingegangen am 2. Mai 1960.) 
Anschrift des Verfassers: Dipl.-Ing. M. Klumpp, Hamburg-Rahlstedt, am Ohlendorffturm 45b. 
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Analoge Systeme von Schalengleichungen 
Von W. Ginther 


1. Einleitung. Seit langem ist es bekannt, da® man die Feldgleichungen eines elastischen Pro- 
blems auf zweierlei Art formulieren kann. Das erste und iibliche Verfahren besteht darin, tiber das 
Stoffgesetz alle KraftgréBen durch die Verschiebungen auszudriicken; die Gleichgewichtsbedin- 
gungen sind dann die Differentialgleichungen des Verschiebungsfeldes. Dual dazu kann man aber 
auch das Gleichgewicht der KraftgréBen automatisch herstellen, indem man Spannungsfunktionen 
einfiihrt. Diese werden dann iiber das Stoffgesetz so bestimmt, daB ein kompatibler Verformungs- 
zustand entsteht. Welche der beiden Méglichkeiten man wahlt, hangt von der Art des gestellten 
Problems, insbesondere von seinen Randbedingungen ab. H. Schaefer! hat gezeigt, daB sich die 
beiden Gleichungsgruppen fiir das System Platte/Scheibe (ebener Flachentrager) im wesentlichen 
nur durch das Vorzeichen der Querkontraktionsziffer unterscheiden, wenn man die drei Verschie- 
bungskomponenten zu den drei Spannungsfunktionen dieses Systems in Beziehung setzt. Daf diese 
Analogie besteht, ist natiirlich kein Zufall, sondern beruht letzten Endes auf der Existenz einer 
Bilinearform, in der die KraftgréBen und die virtuellen Deformationen zu einem Skalar zusammen- 
gefaBt sind, namlich der virtuellen Arbeit der inneren Krafte. Wenn man daher solchen Analogien 
nachgehen will, mu8 das Prinzip der virtuellen Verriickungen im Mittelpunkt der Betrachtungen 
stehen; im Falle kleiner Verschiebungen und Deformationen kann an seine Stelle bekanntlich auch 
das Prinzip der virtuellen Krafte treten. Von diesen beiden Prinzipien aus kommt man ja auch 
leicht zu den dualen Variationsproblemen der Elastomechanik. — Ziel der vorliegenden Arbeit ist 
es, fiir die allgemeine Biegetheorie der schwach verformten elastischen Schale systematisch eine 
Theorie aufzubauen, in der die genannten Analogien klar hervortreten. Daf bei der Schale solche 
Analogien bestehen, ist vor allem von russischen Forschern bemerkt worden; Ausfiihrungen dariiber 
findet man z. B. bei W. S. Wlassow? und V. V. Novozhilov?, wo weitere russische Literatur angegeben 
ist; kiirzlich hat ferner H. Schaefer + die Gleichungen des geraden Kreiszylinders untersucht und 
die vollstandige Analogic beider Gleichungsgruppen festgestellt. Versucht man nun, die Analogien 
fiir eine beliebige Schale nicht nur hinzuschreiben — was nach einigen Vernachlassigungen mit den 
Mitteln der Tensorrechnung nicht allzu schwierig ist —, sondern systematisch herzuleiten, so zeigt 
es sich, daB die Statik und Kinematik der Schale als die einer sogenannten Cosseratschen Flache 
anzusehen ist, einer Flache also, die aus lauter orientierten Elementen besteht. Mit Hilfe dieses 
Schalenmodells lassen sich, wenn man das Prinzip der virtuellen Verriickungen benutzt, simtliche 
kinematischen und statischen Gleichungen der Schale so einfach und in ihren Analogien so pragnant 
formulieren, da man sie sich ohne Schwierigkeit jederzeit ins Gedachtnis zuriickrufen kann. Durch 
geschickte Wahl der Feldgr6é8en kann man dann auch noch das Elastizitatsgesetz in die Analogie 
einbeziehen und erhilt so schlieSlich die gesuchten beiden Darstellungen des Schalenproblems. 


2. Flachentheoretische Hilfsmittel. Wir beginnen mit der Zusammenstellung einiger Begriffe und 
Formeln der Flachentheorie, wobei wir uns in der Darstellung an A. E. Green u. W. Zerna® an- 
schlieBen: 

Seien g® und q@) Gaufsche Parameter auf einer Flache, und sei, von einem festen Punkte 0 
ausgehend, 


r = r(q%, q?)) (2,1) 


der Ortsvektor zu einem nicht singularen Punkte P dieser Flache hin, so kénnen wir in P ein lokales 
Dreibein definieren durch die ,,MaBvektoren‘ 
or ° 

a; = Ow => qi (i i i 2) (2,2) 
ae ore are Pome. Wiss. Ges. 8 (1956) S. 142. 
. 5. Wlassow, Allgemeine Schalentheorie und ihre Anwendung in der Technik, Berlin 1958. 
V. V. Novozhilov, The Theory of Thin Shells, Groningen 1959. AEE e Sas © 
H. Schaefer, Ing.-Arch. 29 (1960) S. 125. 
2 A. E, Green u: W. Zerna, Theoretical Elasticity, Oxford 1954. 
* Wahrend der Drucklegung dieses Aufsatzes erhielt der Verfasser noch Kenntnis von folgender Arbeit, 


die sich vor allem auf ? stiitzt und die Analogien in der technischen Biegetheorie dii 5 
pee OP ee ae g chnischen Biegetheorie diinner Schalen behandelt: 
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und den Normaleneinheitsvektor 


2 (2,3) 
der Flache (Abb. 1). Mit den Ma8vektoren a; tangieren auch die ,,kovarianten Permutations- 
vektoren“‘ 

e,; = EX 4G; (i = 1, 2) (2,4) 
und die ,,Kriimmungsvektoren“ 

b; =—0,E (a= 1; 2) (2,5) 
die Flache. Geht man nun zu neuen Parametern q() und q@) iiber, so erhalt man die zugehérigen 
Vektoren a;, €; und b; durch eine ,,kovariante Transformation‘ 


= aq* a aq* ae éq* 
On Ge oa Ci a ogi? Di Caan 


(2,6) 


(uber griechische Indizes wird von 1 bis 2 summiert!). Vektoren, die sich in ihrem Index so 
transformieren, nennen wir ,,kovariant indiziert’*. Durch ,,kovariante Projektion“, d.h. durch 
skalares Multiplizieren mit den Ma8vektoren, erhalt man 


1. die kovariante Darstellung des ,,MaBtensors‘‘(,,erster 
Haupttensor der Flache“‘): 


a), = 4,°0A = 4); (2,7) 
mit der Determinante 
@ = yy Agy — (449), (2,8) 


2. die kovariante Darstellung des ,,Permutationsten- 
sors‘*‘ der Flache: 


Cee C50 Cp 55 (2,9) 


lei} = Va fiir l =e U t) (2,10) Abb. 1. Ma®8- und Permutationsvektoren . 
und 


3. die kovariante Darstellung des ,,zweiten Haupttensors der Flache‘: 


b;, = 60 a, = — 0,Foa, = + Eo0;0r = by; (2,11) 
mit der Determinante 
b = by; bay — (by 9)” « (2,12) 
Der ,,dritte Haupttensor der Flache“ ist durch 
¢;, = bob; = Gj; (2,13) 


- definiert. 7 
Haufig ist es zweckmaBig, die Symmetrieverhdltnisse eines Tensors zweiter Stufe besonders 
hervorzuheben: Sei u;, ein solcher Tensor, dann kennzeichnen wir durch 


im if 
uy = = (i + uy ;) (2,14) 


seinen symmetrischen, durch ; 
ita oF (ui, — 443) (2,15) 


Sis 


seinen antimetrischen (,,schiefsymmetrischen“) Teil. Ihre Summe ergibt wieder den urspriinglichen 
Tensor. Danach ist 
A : x ‘es A 
GG, ey =A, yp = by GU = (2,16) 
In bekannter Weise erzeugen wir nun ferner die ,,kontravarianten MaBvektoren“ a‘ (i = 1, 2) als 
Lésungen der Gleichungen 


: , P 0 fiir 1-7, | 
emda t=, fide tala ee 


ee> 
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; = : , ee 
bei einer Parametertransformation werden sie durch die ,,kontravariante Transformation 


=; go 24 2.18 
a’=—a aq (2, ) 


in die neuen kontravarianten MeBvektoren a iiberfiihrt. Dasselbe Transformationsgesetz gilt auch 
fiir die ,,kontravarianten Permutationsvektoren“ 
= EB <a: (2,19) 


Kontravariante Projektion ergibt: 
1. die kontravariante Darstellung des MaBtensors der Flache: 


at = aoa Sa, (2,20) 
Det (ai!) =—, (2,21) 


2. die kontravariante Darstellung des Permutationstensors der Flache: 
ei! — eiog! —e'!, (2,22) 
le] = = fiir 1 -bi, (2,23) 
3. die ,,gemischte Darstellung“* des 2. Haupttensors der Flache: 
bi = b,c'a = ab, (2,24) 


aus der sich die ,,mittlere Kriimmung“ 


= be (2,25) 
und die ,,GauBsche Krimmung“‘ 
K = Det (bi) =~ (2,26) 


der Flache aufbauen. 

Die MaBzahlen dieses Tensors bj haben stets die Dimension einer reziproken Lange; fiir die 
»»Kriimmungslinien“ der Flache, also fiir diejenigen Parameterkurven, welche iiberall die Haupt- 
richtungen des Tensors bj tangieren (und daher ein Orthogonalnetz bilden), wird, mit R, und R, 
als Hauptkriimmungsradien der Flache, bei geeigneter Orientierung des Koordinatensystems 


1 1 
paren s 1 ase hee (2,27) 
dabei ist die Gaufsche Kriimmung als nicht negativ angenommen worden. 
Die Anwendung elementarer Vektoroperationen fiihrt noch zu folgenden, gelegentlich niitzlichen 
Formeln 


e;X E=a;, ey Cap REA (2,28) - 
é;;, ef = 6 Or— oF ol, 6°04 (2,29) 
Kéi—2Hbi+cd=0. (2,30) 


Sei nun \/,(. . .) das Symbol der kovarianten Differentiation beziiglich der Flachenparameter g@) 
und qg®). Es gilt dann 


V,;H =0E =— b; = — be a,, 
(2,31) 
ViG=b,E, Vie; = G X b;. 
Hieraus folgt 
Vi ty = 0, Vi on = 0 (2,32) 
(Mafi- und Permutationstensor sind also kovariant konstant) und ferner 
eXB Wi bg = 0 —>e%h Wes ba, = 0 3 (2,33) 


die Gleichungen von Codazzi fiir den 2. Haupttensor der Flache, die erfillt sein miissen, wenn ein 
eindeutiges Feld von Normalenvektoren der Flache existieren soll. 


hieee> ak 
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ry 


Auf der Fliche sei jetzt ein raumliches Vektorfeld » — v(q”, g) gegeben mit den Zerlegungen 


v=1,0 +vE=v%a,+vE, 


; (2,34) 
%,=00G;, v=veoa, v=voE; 


ee sind also, im Sinne der Tensorrechnung, der Vektor v; baw. v' und der Skalar v zusammengefaBt. 
nlich fiir die kovariant indizierten raumlichen Vektorfelder VY = ¥(q™, gq): 
Y= 4,0 +0 E=v%a,+ EF, Yn, =VjoAn, vy" =vVea™, te oR<. (2,35) 
Offenbar lassen sich auch mehrfach indizierte Vektorfelder leicht herstellen, etwa durch 
C15 = EX Om, = Ohi (2,36) 


fiir solche Felder gelten dann entsprechende Zerlegungen. — Das kovariante Differenzieren von 
Vektorfeldern geht nach (2.31) so vor sich: 


Viv (= dv) = (Vita — b:q v) a® + (V;v + bf v,) E 


= Oe ob; v) Ay + (V: vb v*) E; (2,37) 
Vit = (Virta — Oi 4) a + (Vit + 5} 0,4) EB 
= (Vim — bF u) ag + (Vi, + big 0%) E, (2,38) 


mit leichter Verallgemeinerung auf mehrfach indizierte Vektorfelder. 


3. Kinematik der Schalenmittelflache. Wir denken uns (Abb. 2) die Schale aufgelést in eine — 
zunachst endliche — Anzahl starrer Blécke, deren Mittelflachen sich im unverformten Zustand der 
Schale zur unverformten Schalenmittelflache zusammenschlieBen; sie seien dreh- und verschiebbar 
federnd gegeneinander so abgestiitzt, da in der Ausgangslage alle Federn entspannt sind. Durch 
Grenziibergang erhalten wir als kinema- 
_ tisches Modell der Schale eine Flache, die 
geometrisch mit der Schalenmittelflache 
zusammenfallt, aber in jedem ihrer Punkte 
durch die raumliche Lage des zugehérigen 
Blockes orientiert ist. Jeder ,,Punkt‘’ die- 
ser Flache kann sich mit 6 Freiheitsgraden 
bewegen und ist gegen seine Nachbar- 
»punkte“ federnd so abgestiitzt, wie es das 
(spater zu erérternde) Elastizitatsgesetz der 
Schale verlangt. Statt von einem Kontinu- 
um ,,orientierter Punkte‘‘ kénnen wir natiirlich auch von einem Kontinuum orientierter Koordi- 
natensysteme sprechen, deren Anfangspunkte auf der Schalenmittelflache liegen. Kurz: es handelt 
‘sich bei unserem Modell um eine ,,Cosseratsche Flache‘‘, nach E. und F’. Cosserat', die in einer Mono- 
graphie die Kinematik und Statik solcher, aus orientierten Punkten bestehenden Kontinuen unter- 
sucht haben (in moderner Schreibweise behandelte der Verfasser® die Grundgedanken der Cosserat- 
schen Konzeption). Dieses Schalenmodell — im iibrigen die sinngemaBe Erweiterung des Modells der 
technischen Biegelehre des Balkens — ist nun alles andere als neu: bereits bei E. und F’. Cosserat 
wird die Schale als Beispiel herangezogen, eine eigentliche Schalentheorie jedoch nicht angestrebt; 
die heute nur noch schwer lesbare Darstellung ist nicht kovariant. K. Heun? benutzt, ohne sachlich 
iiber das Cosseratsche Buch hinauszugehen, die symbolische Vektorrechnung, wodurch die Darstellung 
iibersichtlicher wird. J. L. Ericksen und C. Truesdell4 schlieBlich entwickeln, mit anderer Zielsetzung 
als wir, eine Theorie endlich verformter Schalen auf der Basis der Cosseratschen Vorstellungen. — 
Wir werden sehen, da uns gerade dieses Modell zu einer auSferordentlich symmetrischen und daher 
. tibersichtlichen Schalentheorie fihren wird. 


Abb, 2. Cosseratsches Schalenmodell. 


1 E. et F. Cosserat, Théorie des corps déformables, Paris 1909. 
2 W. Giinther, Abh. Braunschw. Wiss. Ges. 10 (1958) S. 195. i ; bpd 
3 K, Heun, Ansatze und allgemeine Methoden der Systemmechanik, Enz. d. Math. Wiss. IV, 2, 11, Leipzig 


1914. 
4 J. L. Ericksen & C. Truesdell, Arch. Rat. Mech. and Analysis 1,4 (1958) S. 295. 
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Die Lagednderungen eines ,,orientierten Punktes‘ der Schalenmittelflache sind: eine Drehung, 
die so klein sei, da® man sie durch einen Drehvektor w beschreiben darf, und eine Verschiebung w; 
es sei |v| klein im Verhialtnis zur Schalendicke t. Die Funktionen 

We w(q, q?)) = q”* Ay + @ E., o= v(q®, q?)) —= 0, a” +vE (3,1) 
sollen iiber der Mittelflache stetig und geniigend oft stetig differenzierbar sein. Unter ihnen be- 
finden sich natiirlich auch die starren Bewegungen 


SEAwW, 7=% +0 Xx (rT) (3,2) 


eines endlichen Schalenteils, mit konstanten Vektoren w, = (10), Up = U(1). Aus den Lageande- 
rungen (3.1) leiten wir die Deformationen als relative Drehungs- und Verschiebungsanderungen 
folgendermaBen her: Die absoluten Anderungen sind 


(dwd)absonut = Ag9 dg* = Vy dg*, —— (dv)absotur = Ag0 dg* = \7gv dq”; (3,3) 
ferner ist 
(ded)absotur = (dW), — (AW) absotus = (A) retain + 2 X ar. (3,4) 
Setzen wir die relativen Anderungen als lineare Funktionen der Fortschreitungsdifferentiale dq’ an: 
(ded) retariv = Xa Iq" > (A) retativ = Ea AQ” 5 (3,5) 
und beachten, da nach (1.2) auf der Schalenmittelflache 
dr = a, dq” (3,6) 
ist, so erhalten wir 
YX = VM; é-—Vwt+a;,xo. (3,7) 


Diese Deformationen sind offensichtlich kovariant indizierte Vektorfelder. Wir zerlegen sie: 
Ni = MEG +E, €, = €,0° +6, F (3,8) 
mit 
M=HX0U, YW=XeH, ee =E0q, e = E 0k. (3,9) 


Aus den 6 Lageandeiungen (3.1) entspringen also 12 DeformationsgréBen, die sich nach (2.37) aus 
den Mafizahlen der Lagedinderungen wie folgt ableiten: 


4/ = V,0'— bio, Xi = Viw + 5:4 0%; (3,10) 
&y = Viy— by vu —e,0, & = Viv + OF vg + e400". (3,11) 


Die zweite der Gleichungen (3.7) baw. die Gleichungen (2.11) kénnen nach der Drehung w bzw. ihren 
MaBzahlen aufgelést werden: 


W = (Vx 0 — Eq) o Y* (3,12) 
mit den dyadisch dargestellten Matrizen 
i due Does 
2 =—EHe The By (3,13) 
bzw. 
ee eat 1 x 
cot = e' (7yv + bf vg — &q) » eae: eB (\74¥3 — Exp) 3 (3,14) 


in (3.12) kommen also nur die GréBen ¢; = €; 0 H und der antimetrische Teil des Deformations- 
tensors wirklich vor. 

Die geometrische Bedeutung der Deformationsmafzahlen (3.10) und (3.11) liegt nach ihrer 
Entstehung auf der Hand: die y;! sind die ,,Verkriimmungen“ der Schale, die x; die Normalkom- 
ponenten der relativen Drehung. Der symmetrische Teil @,, des Tensors €;, beschreibt die in der 
Tangentialebene stattfindenden Deformationen der als Punktkontinuum aufgefaBten Mittelflache; 
seien namlich : 


a; — Or + v) = a; + Vv (3,15) 
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die MaSivektoren der verformten Mittelflache, so wird, bei Vernachlassigung von Produkten aus 
DeformationsgréBen, der neue MaBtensor: 


ai; = a;00; =a, +2 Ein (3,16) 
A 1 
Ey = > (Viti + Vit) — by 0 (3,17) 
1/2 (a;; — a@,)) ist aber gerade die genannte Deformation. Wir merken noch an, daB wegen 
Pees : Oo fir =e: 
oe eee (3,18) 
Polar lt 
fiir die Differenz der kontravarianten MaBvektoren gilt: 
a — gt _- 2 8! (3,19) 
Der schiefsymmetrische Teil é,, ist einem skalaren Winkel wy aquivalent: 
1 hep. shes. x 
paz er’ ean, i lee (3,20) 
und es wird nach (3.11) 
yp = ex Wir Ug — W ° (3,21) 


Dies zeigt, daB eri dem Unterschied zwischen der aus der Mittelflachendeformation stammenden 
»mittleren Drehung“‘ und der Cosseratschen Drehung w um die Schalennormale entspricht. 

Zur Deutung der GréBen ¢,;: Der Normaleneinheitsvektor EH’ der verformten Mittelflache geht 
durch eine kleine Drehung 9 aus E hervor: 


BE’ =E+9xE, (3,22) 


fallt jedoch i. a. nicht in die Richtung der vom Srtlichen Block mitgenommenen Normalen der 
unverformten Mittelflache 


E” =E+oxE. (3,23) 
Aus 
E’.a; =0 (3,24) 
und (3.15) folgt dann leicht 
E’ = E—(EoV/,v) a’, (3,25) 
und weiter mit (3.7) 
EK’ = E” —«,a’. (3,26) 
Daraus bestimmt sich 9 zu 
0 6" e,: (3,27) 


Die Verformungen ¢; sind also ein MaB dafiir, wie weit sich eine bei der Verformung mit- 
genommene Normale der unverformten Mittelflache aus der Normalen der verformten Mittelflache 
herausgedreht hat — es handelt sich somit um die Querschubverformungen der Schale. Sind nun 
die Deformationen x; und €; als Funktionen der Flachenparameter gegeben, so wird es im allge- 

_ meinen nicht méglich sein, aus ihnen eindeutige Lageinderungen w und v durch Integration der 
- Gleichungen (3.7) zu errechnen. Damit dies méglich wird, miissen vielmehr die Deformationen 
gewissen ,, Vertraglichkeitsbedingungen*‘ — natiirlich im Cosseratschen Sinne — geniigen, die sich 

~ von (3.7) her nach (2.31) leicht hinschreiben lassen: 


cP TaXp= 9, e*P (Yq €p + da X Xp) = 0- (3,28) 


Sind sie erfiillt, so liefert die Integration von (3.7), erstreckt iiber eine von P, (Ortsvektor 1) nach P 
(Ortsvektor r) fiihrende Flachenkurve C (laufender Ortsvektor s) die von der speziellen Wahl der 
Verbindungskurve unabhangigen Lagednderungen 


w(t) = elt) + [Xe de 
ns : (3,29) 
U(r) = U(r) + (7) X (¥ — Yo) mee [eg — (” — S) X Xa] dq”; 


mye, We 
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die ausintegrierten Teile sind die starren Bewegungen (3,2). Sehen wir von diesen ab, und beachten 
die Definition (3.5) der Deformationen, so wird 
P 


or) = { (deo) v(r) = J [(dv)retativ — (7 — $) X (dO2)retatio] 5 (3,30) 


Py 


7 


bis auf eine starre Bewegung sind die Lagednderungen in P durch die fiir diesen Aufpunkt berech- 
nete Kinemate der relativen Lagedinderungen zwischen P, und P gegeben. 

Wir ordnen nun jedem Punkte der Flachenkurve C (die insbesondere die Randkurve der Schalen- 
mittelflache sein kann) ein Dreibein von Einheitsvektoren zu: n sei ihr, in der Tangentialebene 


‘ é ds. : : ares ; 
liegender Normaleneinheitsvektor, = qs wot Tangenteneinheitsvektor, und EH wie bisher der Kin- 
heitsvektor der Flachennormalen; n,t und E sollen in dieser Reihenfolge ein Rechtssystem bilden, 
so daB n der Vektor der auBeren Normalen ist, wenn ein Flachenelement im positiven Sinne um- 


fahren wird. Wir zerlegen 
(dod) retativ aa xX ds ) (dv) relativ = eds (3,31) 

nach diesen Richtungen: 

ea is 

(dv) relativ = (y ds) n + (e ds) E + (6 ds) E; 
darin ist 
yn die Anderung der Normalkriimmung, 7p die Anderung der Torsion, 7¢ die Anderung der geo- 
datischen Kriimmung der Integrationskurve infolge der Flachendeformation, y ihre Richtungs- 
anderung in der Tangentialebene, ¢ ihre Dehnung und f = ¢, t” ihre Richtungsanderung in der 
t — E-Ebene. 


Wir setzen (3,32) in (3,30) ein, nehmen im zweiten Integral eine partielle Integration vor und 
lassen die ausintegrierten Terme fort, die wieder starre Bewegungen darstellen. 


Es bleibt schlieBlich 


iz 
(0) = [Gin +p +168) ds. 
3 (3,33) 
a 
o(r) = [ [et—@—s) (xv 2 + Xn t + Xe B) +5, (@n—y B)|| ds. 
Py 


Kinematisch bedeutet diese Umformung, daB bei der Berechnung der Verschiebungen aus den 
Deformationen die Schubverformungen # und y nicht selbstandig auftreten, sondern in ,,kine- 
matisch aquivalente Zusatzverkriimmungen“ umgewandelt werden. Wir kénnen aber nicht nur v, 
sondern auch — was fiir das Randwertproblem wichtig ist — die Anderung 0v/dn senkrecht zur 
Flachenkurve C berechnen: es ist 


ov : . - 
at Va® mit n'=noa'; (3,34) 
dies gibt, in (3,7) eingesetzt, 
ava av 
ey B ov 
Non = Eqn n > tor =eqgn*t? + wok, Bor = é,n* —wot (3,35) 


oder, wenn wir in dyadischer Darstellung die Matrizen 


CE =c,,a° a’ +e,0° E (3,36) 
(die Deformationsmatrix) und 
’=fLt— tL 
benutzen, 
ow 
ma MoE + wok; (3,37) 


® ist natiirlich aus (3,33) zu entnehmen. 
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Wir kehren nun zu den Vertraglichkeitsbedingungen (3,28) zuriick und geben noch ihre Tensor- 
darstellung an: 


e*F (7a %p1— b) v_) = 0, OPW a Xa + ba, 4h) = 0; (3,38) 
e*? (4&1 — bai gp — Car Xp) = 9; e%B (Va &g + b4 en, + €g, ¥p4) = 0. (3,39) 


An Stelle der Verkriimmungen X; wollen wir nun — im Hinblick auf spatere Entwickelungen — 
andere VerkriimmungsgréBen 


x = '*@, + E (3,40) 
bzw. 
Pi = 0, a + 0; E (3,41) 
einfiihren durch 
R= el ¥,, AX; = ey %” (3,42) 
bzw. 
P=XXE+(XicH)E, y,=EX e;+(Eoe) E; (3,43) 
in den MaBzahlen: 
Soran ee Ny ene! ee ay Oy, 0° (3,44) 
bezw. 
Cit = 1a Kis» ei, =e! on, OC, %> Xi = Oi (3,45) 


die Verkriimmungen x‘ und pe; hangen dann in ihren Mafzahlen durch 


wit —= — eld @lB Cap Oi = ea: @, xP ? wi =e 0q3 0; = Cg x” (3,46) 

zusammen (unseren x! entsprechen die (—x’') bei Schaefer!). Einen schematischen Einblick 
in die gegenseitige Zuordnung der GréBen y;, x; und pe; gibt die Abb. 3; man sieht, daB o;, der in 
der Schalentheorie gewéhnlich verwendete Verkriimmungstensor ist (unser Q;, korrespondiert 
annahernd mit q;, bei Green und Zerna?, und mit x;, bei W. Fliigge*). Wir stellen jetzt die- 
_ jenigen kinematischen Gleichungen fiir die GréBen €; und x; bzw. e; zusammen, die fiir uns von 


g 


[¢7 -7,,-@,] 


ne), ete 
72752 

| i £-a4] 
aa Tia gl? 


[2 Ae" -Or,| 
[x “ee Gey] 


Abb. 3. Schema der Verkriimmungen. 


weiterem Interesse sind. Bei Benutzung der x; ist dies das gesamte System der kinematischen 
Gleichungen; dies ist bereits ein Teil unseres endgiiltigen Systems analoger Schalengleichungen 
(wir kennzeichnen sie durch einen Stern): 


H,=e% Vg, G=Viv+a,X w (3,47)* 

bzw. 
nit = (TF, a'— bh), x = 6 (Van + by w*) 5 (3,48) * 
6p =Viy— by y—e@, & = Viv + bi vy + eq 0% (3,49)* 


1 Siehe FuBnote 4 von Seite a 
2 Siehe FuBnote 5 von Seite 160. 
3 YW. F litgge, Statik und Dynamik der Schalen, 2. Aufl., Berlin-Gottingen-Heidelberg, 1957. 


168 W. Giinther: Analoge Systeme von Schalengleichungen Ingenieur-Archiv 


und die Vertraglichkeitsbedingungen 
Vee =O e*P 74 €, +d, X x” = 0 (3,50) * 
bzw. 


NV i be 0 Var + bap xtP= 0 5 (3,51)* 


a 
eB. Eg; — bar Es) — eg X= 9, P(g &g + 0) faa) + eng neh = 0. (3,52) * 


Die erste der Gleichungen (3,52) kann nach den x! aufgelést werden: 


x! = e!4 eh (Va Epa by. Es) 4 (3,53)* 


diese GréGBen sind daher bei einem vertraglichen Deformationszustand durch die Deformationen 
€; festgelegt und konnen, wenn es gewiinscht wird, ganz aus den Vertraglichkeitsbedingungen 


eliminiert werden. 
Fiir den konventionellen Verkriimmungstensor 0,;, geben wir seine Darstellung durch die Lage- 


anderungen an: 


0:1 = 41a (Vi Oo — BF @) (3,54) 
und rechnen ferner die zweite der Vertraglichkeitsbedingungen (2,39) auf ihn um: 
e*8 (94g — bd &5, — Va &p) = 0- (3,55) 9] 
Dies zeigt, daB bei einem vertraglichen Deformationszustand der Tensor 
Ou == Op — OF 8g — Vike (3,50) 


symmetrisch ist. Der Grund dafiir ist leicht einzusehen: Driicken wir namlich die rechte Seite von 
(3,56) nach (3,49) und (3,14) durch die Lageinderungen aus, so finden wir 


On =— Vi Viv — bf V, ea bf Vi % — Vi br: +e, 8, (3,57) 
und das ist, wie eine kleine Rechnung zeigt, identisch mit 
07 =— ViVi DOE (3,58) 


mit offensichtlicher Symmetrie in den beiden Indizes. Die Tatsache, daB (3,58) in den Tensor der 
Plattenkriimmungen tibergeht, wenn die Schale zur ebenen Platte entartet, reizt zu einer anschau- 
lichen Deutung auch fiir die Schale. Hierzu greifen wir auf die Darstellungen (3,15) und (3,25) 
fiir die MaBvektoren und den Normaleneinheitsvektor der verformten Schale zuriick und berechnen 
uns den zweiten Haupttensor der verformten Mittelflache 


bis = — 0,E' 0 a; = — V7; [E — (Eo Va v) a*] 0 [a, + V, 2] 5 (3,59) © 
vernachlassigen wir wieder die Produkte von Verschiebungen, so ergibt eine kurze Rechnung 
bi = by + Vi Vi Vo E = by — Gi). (3,60) 


Der symmetrische Verkriimmungstensor @,, ist also gleich der negativen Anderung des 2. Haupt- | 
tensors der Mittelflache infolge ihrer Deformation. 


4. Statik der Beanspruchungsgréfen. Die auBeren Schalenbelastungen (Volumkrafte, Krafte an — 
den Schalenlaibungen) reduzieren wir fiir jedes Schalenelement in seine Mittelflache df hinein — 
zu Kraften pdf und Momenten qdf. Langs des Schalenrandes (wir setzen eine einfach berandete _ 
Schale voraus) mégen ferner an den Randelementen ds der Mittelfliche die Krafte dK —K ds und 
die Mamente dM = Mds angreifen. Der Rand seiin dem durch (2,3) gegebenen Sinne orientiert. 
Nach dem Prinzip der virtuellen Verriickungen ist im Gleichgewichtsfalle die negative virtuelle 
Arbeit der inneren Krafte und Momente, also hier der Beanspruchungsgréfen, bei der Defor- — 
mation gleich der virtuellen Arbeit der iuBeren Krafte und Momente. Als virtuelle Lageainderun- 


gen sind hier natiirlich die unseres kinematischen Schalenmodells anzusehen. Wir definieren nun __ || 


»virtuelle Deformationen‘* durch 


6x; = Vi (dea), de, = V7, (8v) + a, x dw, 4) 
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setzen also differenzierbare virtuelle Lageanderungsfelder voraus, und setzen fiir das Prinzip der 
virtuellen Verriickungen an 


— 0A) = — Jf o*? (Ky 0 b€s + My 0 dy) df = ff (p odv + q 0d) df 
+ dK odv + dM od, (4,2) 


womit die zu unserer Kinematik gehorenden »»Kraftbeanspruchungen“ N, und ,,Momentenbean- 


spruchungen“‘ M; eingefiihrt werden. Setzen wir (4,1) in (4,2) ein und formen nach dem Stokesschen 
Satz um, so entsteht: 


— SJ {(e%? Va Kp + p) 0 dv + [e*? (Va Mg + ay X Kz) + q] 06} df 
= f {(dK — K, dq*) od6v + (dM— M, dq*) o dw}. (4,3) 
Daraus folgen, weil im Cosseratschen Kontinuum die virtuellen Lageanderungen dv und dw 


in allen ihren MaBzahlen unabhingig voneinander gewahlt werden kénnen, die Gleichgewichts- 
bedingungen, giiltig im Bereichsinneren, 


e*P \7, Kz = — p, e*8 (\7, Mp + a,x Kz) =—q, (4,4) 
und die Beziehungen 
dK = K,, dq* , dM = M,, dq’ , (4,5) 
gultig auf dem Bereichsrand. Nimmt man einen von Parameterkurven begrenzten Bereich, ‘so 
liest man aus (4,5) die statische Bedeutung der von uns gewahlten Beanspruchungsgréfen 
K, = K;* a* + K,E (4,6) 
und 
, M;= M;,,a*+ M,E (4,7) 
ab (vgl. Abb. 4 und 5); im wesentlichen sind: die K;' Membran-Schubbeanspruchungen, die K;! 


fiir 1 =-i Membran-Langsbeanspruchungen, und die K; Querkraftbeanspruchungen; die M,;; 
Biegemomente, die M;, (J ==1) Drillmomente, die M; Momente um die Schalennormale. Die 


Ma ie mi\ 


[Mp, mr 
SESE 


_ 4 > 
ook 1] 28 
Ww sae [My 
Efe th 2 
\w ) K,3| (Mg. | 
Abb. 4, Schema der Schalenkrifte. Abb. 5. Schema der Schalenmomente. 


6-fach funktionale statische Unbestimmtheit des Gleichgewichtsproblems der Schale auBert sich 
einmal darin, daB statisch fiir die 12 Beanspruchungen (4,6) und (4,7) nur die 6 Gleichgewichts- 
bedingungen (4,4) verfiigbar sind, zum anderen auch darin, daB die Zerlegungen (4,5) der 
Randkrifte und -momente nicht eindeutig bestimmt sind, denn die homogenen Gleichungen 


K,dqg@=0, M,dq*=0 (4,8) 


haben je 3 willkiirliche skalare Lésungen. Wir beschranken uns jetzt auf den Fall, dafs die Schale 

von Belastungen durch Flachenkriafte und -momente frei, also nur am Rande durch ein Gleich- 

gewichtssystem von Kraften und Momenten belastet ist; darauf ]aBt sich, durch Abspalten einer 

Partikularlésung, der allgemeine Fall stets zuriickfiihren. Dann haben die Gleichgewichtsbedin- 

ungen 

a eX? Ky =0,  o*? [Va My + a X Ky] = 0 (4,9) 
12 
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die Form der Vertraglichkeitsbedingungen (3,28), wobei den relativen Drehungsanderungen Xi 
die Kraftbeanspruchungen K;, den relativen Verschiebungsanderungen €; die Momentenbean- 
spruchungen M;, entsprechen. Es folgt, daf ein den Gleichgewichtsbedingungen (4,9) geniigendes 
System von Beanspruchungen analog zu (3,7) durch vektorielle Spannungsfunktionen $2 und & 
dargestellt werden kann, welche den Lageanderungen und » entsprechen: 


K,=Vi2, M;=ViP +4, xQ. (4,10) 

(In nicht kovarianter Form wurden die Gleichungen (4,9) und (4,10) auf anderem Wege vom Ver- 

fasser! hergeleitet.) Wiederum kann die zweite dieser Gleichungen nach Q bzw. seinen Maf- 
zahlen aufgelést werden (vgl. (3,12) bis (3,14)): 

Q =(V,2 — M,) o U* (4,11) 


bzw. 


Qi = e* [VP + b2 O, — Mg], 


l 2 (4,12) 
Q = 2 Ores D, = Mgsl 3 
in (4,11) kommt also der symmetrische Teil des Momententensors nicht wirklich vor. 
Ist nun ein Gleichgewichtssystem von Beauspruchungen der randbelasteten Schale gegeben, 
so kénnen, entsprechend (3,29), die zugehérigen vektoriellen Spannungsfunktionen durch Inte- 
gration von (4,10) nach Vorgabe ihrer Anfangswerte eindeutig berechnet werden: 


P 
Q(r) = 219) ea K, dq , 


(4,13) 
P 
D(r) = P(r,) + Q(r,) X (r — 79) ae [M,,—(r— s) X K,] dq’. 
Die ausintegrierten Terme ° 
Q=Ar), $= Hr) + Q(m) x (r—%)- (4,14) 


welche den starren Bewegungen (3,2) entsprechen, sind die ,, Nullspannungsfunktionen“ der Schalen- | 


statik, denn sie erzeugen, in (4,9) eingesetzt, offenbar einen beanspruchungsfreien Zustand der 
Schale. Mit Hilfe von (4,10) und (4,5) kénnen wir nunmehr auch die Randwerte der Spannungs- 
funktionen aus den gegebenen Randlasten ermitteln: wir nehmen, ausgehend von einem ihrer 
Punkte, die Randkurve als Integrationsweg und erhalten unter Fortlassung von Nullspannungs- 
funktionen 


Q (7 Rand) = [aK ®(7 Rana) a4 [dM — (r — s) X dK]; (4,15) 


diese Randwerte sind also durch die fiir den Aufpunkt P berechnete Dyname der zwischen P, 
und P angreifenden Randkrafte und Randmomente gegeben — man vergleiche hiermit (3,30). 
Wir zerlegen nun wieder nach den orthogonalen Einheitsvektoren n, € und E: 


dK = K ds = (Lds)n + (T ds) t + (Qds) E, 
dM = M ds = (Mp ds) n + (Mz ds) t + (My ds) E ; 

darin sind, auf die Kinheit der Randkurvenlange bezogen: L die Membranlangskraft, T die Mem- 

branschubkraft, Q die Querkraft; Mp ist das Drillmoment, Mz das Biegemoment und My das Mo- 


ment um die Schalennormale. Ebenso wie in (3,33) ergibt sich dann, bis auf belanglose Nullspan- 
nungsfunktionen 


QUrRana) = f (Ln + Tt + OH) ds, 
P, 

| , : (4,17) 

P(rRona) = J M;,t —(r—s) X (Ln + TL+QE) +3, (Mym—MpE)|| ds. 


Statisch bedeutet diese Umformung, da8 bei der Berechnung der Randwerte der Spannungs- 
funktionen aus einem gegebenen Gleichgewichtssystem von Randlasten die Momentenanteile 
My und Mp nicht selbstandig auftreten, sondern in statisch aquivalente ,,Zusatzkrafte’* umge- 
wandelt werden; in der Plattentheorie (My =0, OE/0s = 0) ist dies die »»Kirchhoffsche Zusatz- 


(4,16) 


1 W. Giinther, Abh. Braunschw. Wiss. Ges. 8 (1956) S. 111. 
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kraft“ Q* = —0QM)/ds. Ferner kénnen wir wieder, etwa fiir di i 
a ee er, etwa fur die Randkurve, die Normalableitung 
o@ rs 
UR ae My, n* nP , 
oD 
toe = Magn* 4 Q oH, (4,18) 


oD 
Bo7-=M,n*—Q ot, 


bzw. in dyadischer Matrizenschreibweise mit 


WM = My,a* a + M,a*E: (4,19) 
o® 
Be OM Rok. (4,20) 


wobei Q aus (4,17) zu entnehmen ist. Hierzu vergleiche man die Formeln (3,35) bis (3,37). 

Unser nachster Schritt besteht darin, die Kraftheanspruchungen K; (4,6) durch 

Ni = e'*K, , K; = €;.N’, (4,21) 
in den MaBzahlen: 
Nae K!, Kit seq NY, 

‘ ; 4,22 

Nese eo? Ki K, = e,; N° oa) 
entsprechend (3,43) zu transformieren. Damit haben wir die in der Schalentheorie gebrauchlichen 
Schalenkriifte N* (Membranbeanspruchungen) und Ni (Querkrifte) eingefiihrt. Wir erhalten 
dann die folgenden statischen Gleichungen, die den kinematischen Beziehungen (3,47) bis (3,52) 
entsprechen und mit ihnen zusammen zu unserem System analoger Schalengleichungen gehéren: 


Ni =e Wg Ses M;=V;®8 +4; x Q (4,23) * 
bzw. 
Nil = gio OSS ; 
ee eae (4,24) 
Ne == 6° (\/ 5 Qe by g Q?) 2 
M;, —. V ai OD, — by, D—e;, Q ? (4,25) * 
M; = V/,® + b? ®, + ¢, 2" ; 
VeN=0,)_ 288 7, M, +a, X N°; (4,26) * 
in den MaBzahlen: 
Va NN —bLN® = 0,0 "V7, N* +8,, N° = 0; us 
e%F [74 Mgi — bx1 Mg] — ear N* =0, 
x (4,28) * 
eXB [Va Mg + bx Maz] + x p Nu p =0. 


Die erste der Gleichungen (4,28) kann nach den Querkraften N' aufgelést werden: 
N! = el4 e8 [V, M?* — by, Mg). (4,29) * 


diese GréBen sind daher im Gleichgewichtsfalle bereits durch die Momente M; festgelegt und 
kénnen, wenn es gewiinscht wird, ganz aus den Gleichgewichtshedingungen eliminiert werden. 
Fiir den endgiiltigen Anschlu8 an die konventionelle Darstellung der Schalentheorie (in der 
Analogien nur noch schwach erkennbar sind) rechnen wir auch noch die Schalenmomente M; auf 
neue Momente mi‘ um [vgl. (3,44)]: 
mi = e* [M, x E+ (M, cE) E], (4,30) 
M;, = ¢,i [E < m* + (Eo m*) E] 
in den MaBzahlen 
mi! = mi oa! = e'* eb My, 
M; = M; ° a, = Cui BI m*B 9 
M,; = M, cE =e; M"*. 
12* 
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Die Gleichgewichtsbedingungen (4,26) verwandeln sich dann in 
Vip = LN 07, Vg Nt bg NO 0; (4,32) 


Vz m*!— N! + e%'by,m? =0, 


(4,33) 
Va" + Gap [N+ bf mis] = 0, 
wie sie, ohne die Momente mi’, z. B. bei Green und Zerna angegeben sind. Die Abb. 4 und 5 zeigen 
schematisch die gegenseitige Zuordnung der in diesem Abschnitt definierten Beanspruchungs- 
groBen. 

SchlieBlich geben wir noch zwei leicht zu verifizierende Formen der negativen virtuellen Arbeit 
der BeanspruchungsgréBen an: 


—§A® = ff (N* ode, — M, 06x’) df, 


(4,34) 
— $A = ff (N% ode, +m odp,) df. 

Auch iiber diese Ausdriicke hatten wir natiirlich zu den entsprechenden Gleichgewichtsbedingungen 

gelangen konnen. 


5. Einbettungstheorie: Geometrie, Statik, Kinematik. Die bisher gefundenen Gleichungen er- 
fassen noch nicht das ganze Schalenproblem; sie miissen durch ein Stoffgesetz erganzt werden. 
Ein solches Stoffgesetz kennen wir aber nur fiir das rdumliche Punktkontinuum, némlich — im 
inear-elastischen Fall — das Hookesche Gesetz. Wir miissen daher die rdumlichen Deformationen 
durch die der Mittelflache ausdriicken, und die raumlichen Spannungen zu den Beanspruchungs- 
groBen zusammenfassen, was nur fiir diinne Schalen sinnvoll ist. Dazu fiihren wir eine Koordinate 
q°) = z senkrecht zur Schalenmittelflache ein und ordnen jedem Punkte des Schalenraumes den 
Ortsvektor 


RQ”, @)5 2) = rq, @) +2 EQ, @), —zSeS4+y (5,1) 


zu (wir ubernehmen wieder die von Green und Zerna verwendeten Bezeichnungen). Die MaBvektoren 
dieses raumlichen Koordinatensystems (die Indizes i, I, . . .;, 8, .. . nehmen wie bisher die Werte 1 
und 2 an) sind 


9: = OR =O" + 20H = a,—2b;, = 


93 =0R=E, 
woraus sich fiir den kovarianten raumlichen MaBtensor ergibt 
Sit = Gi Gr = % — 225, + Poy = Bi» | 
ts Gok (5,3) 
&3 =~EoHk=1. 


Die kontravarianten raumlichen MaBvektoren g' werden, als die Lésungen der Gleichungen 


g og, =O, 5,4 
Potenzreihen in z: (5,4) 
g =a'+2b'+2b,b%+...; (5,5) 

ferner ist 
g° = E 9 (5,6) 


so daB der kontravariante raumliche MaBtensor die MaBzahlen 


gs" 2.ghog —a 4-225)! 32! 
gs? =g ok =0, (5,7) 
g?= FoR =1 
erhalt. 
Ist 


& = 811822 — (812)® (5,8) 


ee ee ny rE ee Oa 
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die Determinante des riumlichen MaBtensors, so gilt 


hoy 1 2H KS... ; (5,9) 
h ist das Verhaltnis des Inhaltes 
Af (q™, g); 2) = (gy dq X gz dq?) oE (5,10) 


eines zur Mittelflache parallelen Flachenelementes zum Inhalt 
df(q™, q®); 0) = (a, dq! X ay dq?) o E (5,11) 


des entsprechenden Elementes der Mittelflache. 
Wir beginnen nun mit der angekiindigten Reduktion und fassen zunachst die rdumlichen 
Spannungen 7'!, ‘3 (wir nehmen an, daB 1°? iiberall verschwindet) zu den BeanspruchungsgréBen 


I — 
der Schale zusammen: Die an einem Element df =/g;; dqi dz (i +1) angreifende Schnittkraft ist 
l l 
dK = (1%8 gg + 1% E) ny df ; (5,12) 

darin ist 

l I 

Pet rt on hea 
die raumliche (!) Zerlegung des in den Schnitt hineinweisenden Normaleneinheitsvektors n der 


Flache df: 


1S i I Nee ere a (5,13) 
gu V820 \'s Vu 
was sich zu 
b= ef; Ve g* g? 
na, ——— e,,h—— (i +)) (5,14) 
z Va Veii 4 Veii 
zusammenfassen laBt.. Daher wird 
ee ae ee ee (5,15) 
Vai 
und 
1 , 
dK = e,; (hr? Ip + hr? E) dq dz. (5,16) 
Wir definieren mit Green und Zerna durch 
Fada In = oP Qs » oth = 7a gn ° a = 7%*8 —z be thas (5,17) 


. ° t *. (74 ie . 
einen (i. a. nicht symmetrischen) ,,reduzierten Spannungstensor“ o*'; dann wird 


dK = 0,;(ho%? ay + hr E) dg dz. (5,18) 
Nach (4,5) ist am Schnittelement ; rz 

ds = ya;; dq (5,19) 

der Mittelflache 
K,dg = f[dkK (i+]), (5,20) 

() 
Iso 
: Hi = tui( J bo°? dey + eau ( J 2°? de) B. (5,21) 
(2) (z 


Durch Vergleich mit (4,6), (4,21) kommt 
AE Ore LAI ae cee (5,22) 
() (2) 
in den MaBzahlen: 


Nil = fhoildz, Ni= [hrdz. (5,23) 
@ eas oe”) 
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Ahnlich findet man 


dM =z EX dK =e,¢,2ho atdgid: (i +1); (5,24) 
mit 
M, dq = f dM (i =-1) (5,25) 
(2) 
folgt dann 
M; = é:g eae ( fzeho*? a) a’, (5,26) 
(2) 
in den Mafzahlen: 
M | = e754 e1B fzho? dz, Mo 0. (5,27) 
(2) 


Vergleicht man dies mit (4,31), so erhalt man fiir den konventionellen, nach den Spannungen 
indizierten Momententensor 

m= [zh ods. (5,28) 

@) 

Die Momente M, um die Schalennormale verschwinden in dieser, von den (nicht-singularen!) 
Spannungen des raumlichen Punktkontinuums hergeleiteten Beanspruchungstheorie, ebenso 
natiirlich die ihnen dquivalenten Momente m'. Dann entartet die zweite der Gleichgewichts- 
bedingungen (4,28) in eine algebraische Beziehung: 


exp NX? = e%? bf Ma, (5,29)* 
die zweite der Gleichgewichtsbedingungen (4,33) in die bekannte Gleichung 
2 Ni = NY NY — bi met — bE mt’. (5,30) 


Ferner verlieren jetzt die Flichenkomponenten 2‘ der Spannungsfunktion 92 ihre Selbstandigkeit; 
denn nach (4,12) werden sie von der Spannungsfunktion ® abhangig: 


M=e*(V,B+ be P,| x .3))5 
und daraus, fiir eine Flachenkurve mit dem zugehérigen Dreibein n, f, E: 
of =——— sk, (5,32) 


was auch aus (4,18) entnommen werden kann. 


SchlieBlich geht (4,29) mit (4,31) in die Gleichung 
NE =o ot! N/M = Vint (5,33)* 


ber, die formal mit der entsprechenden Plattengleichung iibereinstimmt. 


Wir wenden uns nun den raumlichen Deformationen zu, die wir unter der Voraussetzung auf. 


die Deformationen der Schalenmittelflache zuriickfiihren wollen, daB deren Querschubanteile & 

verschwinden (diese Annahme ist nur erforderlich, um die Analogie zur Gleichung (5,27) herzu- 

stellen), womit die Querkrafte N‘ zu Reaktionskraften werden. Aus ¢; = 0 folgt zundchst, wie 

ohne weitere Rechnung von den entsprechenden statischen Gleichungen her iibernehmen 
énnen, 


exp xr? = erF bbe, ; (5,34) * 


diese algebraische Vertraglichkeitsbedingung entspricht der algebraischen Gleichgewichtshedingung 
(5,29); ferner wird, analog zu (5,31), 


cot = e'* [Va v + bf og] ; (5,35)* 


die Drehungen um die Tangenten an die Koordinatenlinien der Mittelflache werden also vom Ver- 


schiebungsvektor v abhiangig. Formal mit (5,32) tibereinstimmend wird fiir eine Flachenkurve 


Q 
w of =— So E. (5,36)* 


Der statischen Beziehung (5,33) entspricht auf der kinematischen Seite 
nb = ol4 8 725, (5,37) *. 
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und der in (3,56) eingefiihrte symmetrische Verkriimmungstensor 6;, reduziert sich schlieBlich auf 


Oi = On — OF eq; 
daraus folgt areas - (5,38) 


2 0:1 = On — O14 = OF B15, — bf &, (5,39) 
[was man mit (5,30) vergleichen mége], und 
E x 1 ‘ 
Gin = Qi — > (OF Eq + DF e;q) - (5,40) 
Die réumlichen Deformationen y,, sind nun durch 
ree eis) 
a eet 2 ears (8:1 — i) (5,41) 


definiert, wobei gj; der raumliche MaBtensor nach der Deformation ist. Zusammen mit (5,3) 
liefert das 


1 c tA 1 / 
ite ey (4:1 — 4) — 2 (bi; — by) + 2 - x (ei — ¢;;) . (5,42) 
Greifen wir auf (3,16) und (3,60) zuriick, so entsteht zundchst 
“~ ~ 1 / 
Yaa = En 1 2 On + 2+ > (Cir — cy) - (5,43) 
Nun ist 
cin = BE big = a” Bi, big = (aX? — 2 848) (b;4 — Gin) (big — Gia) » (5,44) 
wobei (3,19) verwendet wurde; naherungsweise wird daher, auf Grund von (5,38) 
Ci — Cy = — (BF Qa + BF ia) » (5,45) 
also 
n~ ~ 1 
Vin = Fi + 21 — 2 + > (OF Gia + OF Oia) - (5,46) 


Die gewogenen Mittelwerte von y,;, werden 


2 1] mn ? 1 
Saree dz = E— 79° > (BF Ora + bi Oia) » 


12 | aie (647) 
Ou = a | 2%ud = Cu— 3 (OF Gia + BF ea) - 
2) 
Dieses Ergebnis legt es, im Hinblick auf (5,39) und (5,40), nahe, fiir den unsymmetrischen Teil ey 
des Tensors ¢;, anzusetzen 
“ ei ee 5 (0 01a — bi ein) (5,48) 


bzw. in x’ geschrieben 
x 12 
oF bap = ipap 6.49)" 


vorausgesetzt, daB wir uns mit dieser Annahme nicht in einen Widerspruch zum Elastizitatsgesetz 
begeben, das wir erst noch aufstellen wollen. Das ist, wie sich zeigen wird, nicht der Fall, Mit 
der Annahme (5,49) wird jetzt auch noch die Drehung w um die Schalennormale abhingig vom 
Verschiebungsvektor v; eine langere Rechnung zeigt, daB 


1 1 te 
w =—y— oe" Va le — 5h (Vi + bf “| (5,50) 
1—7,k ‘ 


wird. Fiir diinne Schalén ist es stets zulassig, GréBen von der Form t- b} als so klein gegen 1 an- 
zusehen, da® man sie héchstens in der ersten Potenz mitzunehmen braucht. Dann wird 


1 N 2 
wr se? Vy (e — 54 Vs ) (5,51) 
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Einfacher und fiir praktische Rechnungen auch wohl in der Regel erforderlich ist es anzunehmen, 
daB der Tensor ¢;, symmetrisch ist: . 
JX x 
Ei] => El bzw. Ej] = 0 . 


(5,52) 


Dies fiihrt sofort auf 


wo =e? Yu Xp 5 (5,53) 


w wird zur ,,mittleren Drehung“, die in den Membrandeformationen steckt. Mit e; = 0 und (5,52) 
entartet dann unser Cosseratsches Kontinuum endgiiltig zum gewéhnlichen Punktkontinuum, in 
dem nur noch dem Verschiebungsvektor v selbstandige Bedeutung zukommt. 

Fiir die raumlichen Schubdeformationen y;, finden wir gemaB (3,26) leicht 


Vis = & @* og; - 
Sie verschwinden also gleichzeitig mit den ¢;. SchlieBSlich folgt mit (3,26), daB die Deformation 
Y33 = 93°93;—1 (5,95) 


in unserem Modell ebenfalls verschwindet. 


(5,54) 


6. Die vollstandigen Schalengleichungen. Das Elastizitatsgesetz des isotropen Werkstoffes im 
linear-elastischen Bereich lautet 


WH 7 (Sete ee ene (6.1) 


1— 7 


wenn wir, wie tblich, einen ebenen Spannungszustand voraussetzen, also dariiber hinwegsehen, 
daB dazu die kinematische Gleichung yz, = 0 im Widerspruch steht. Aus (6,1) gewinnen wir 
mit (5,17) leicht fiir den reduzierten Spannungstensor o*! die Formel 

ho! = 7s [(L—») 8'* (hg? oa!) + 7 g%? (hg' oa')| yas, (6,2) 


1—y 


in der wir nun rechts alles nach Potenzen von z entwickeln und wieder Glieder von der Form zb! | 


nur bis zur ersten Potenz mitnehmen. Man bemerkt dann, da die Unsymmetrie (5,48) des Tensors 
€;, der Membrandeformationen von dieser Naherung nicht mehr erfaBt wird, so daB dieser als 
symmetrisch angenommen werden kann. Nach Integration gemaB (5,23) und (5,28) iiber die 
Schalendicke ¢ stellen die folgenden Gleichungen das Stoffgesetz der Schale dar: 


eee 1—r oe x 12 3 1 
eu (a —v) ai al6 lay, + ate b2 ose + 5 Of ex —2 H 04,)| 
; Pre ; ; 
+ylaile+ © [ail b? on + bi 9p —2 Hai al (6,3) 
il J EF ix lp As ar A 
are reat v) aia LOag + (bx €,5 + bf €,, —2H Eqs) | 
SGA eS a ae Beet | (6,4) 
mit 
é=a*? eye, 9 =a? og. (6,5) 


Weil nun nach (5,30) der antimetrische Teil des Langskrafttensors bereits durch die Schalenmomente 
gegeben sind, bilden wir von (6,3) den symmetrischen Teil 


me Oe 
Exe + 7g (be Crp + 53 Qin —2H xa) 


Ajil 


Et : 
— ae (a —y) ai* al b 


i ee: ; F : 
+vlail e+ 1 (al! b? 0,5 + b# 9 —2 Hail al (6,6) 


Fiihren wir daher die geometrischen Tensoren 


ane irs i f - i 
Oe a reer) Sec raeM | els mee ai) gis 4 (6 — Ha!) a? = bale + hae (6,7) 
und ‘ 


an 


tLrs  jilrs eel 5 aA rs . i 
Qu = Qe'' = (6 — Hall) at + (6 — Har) ai! BH! gre 4. Br git (6,8) 
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in, so wird 
Sai Nt : ee Po. 
N =e {a v) ja al B Exe + 12 Poe Ons] aly 
2 ieee Ee SP {(1 ) ia GIB pib%bs i] il,xaB a 
12(1 —»2) v) [a'*a Cap + Exel tv [ait o+Q Exp] } : (6,10) 


il P Ail, ap 
qué + 752 Ca B 


, (6,9) 


5 < aay s 
Man beachte, daB der soeben eingefiihrte Tensor b" der Deviator des zweiten Haupttensors ist, 
in den Koordinaten der Hauptkriimmungslinien: 


* lel 1 “ i 
in at ole) ak, %): moi —0. 6,11) 


Von (6,9), (6,10) her 1aBt sich das Problem der Umkehrung des Stoffgesetzes leicht lésen. Wie 


man durch Einsetzen bestatigen kann, wird namlich, im Rahmen unserer Vernachlassigungen 


2 x 
a;,m — Oi a B inal (6,12) 


| ) 
1 Zs 
= (E+) [aie a9 NOY — Pir. ag m*?] —v [ay N— Qui, 09 m*°]} , (6,13) 


As a pete 
ta “1B 12 


ou = get +9) 


Pit, ap Art| —y 


ei] 


mit 
m=d,gm*P, N —=ay, NB, (6,14) 


Hier mége eine Bemerkung iiber die Unsymmetrie des Momententensors mi‘! eingeschaltet 
werden: Nach (6,10) ist 


ea i i Ee i 1 
aes 1) — = jaa) (tL —») 0"? (Ong — Opa) (6,15) 


wegen der algebraischen Vertraglichkeitsbedingung (5,39) wird daraus nach (6,13), unter Beachtung 
von (5,30) 
2 ; oh 
mi! — mi =H (bi Ne = bh Ne). (6,16) 
womit die Analogie sowohl zur algebraischen Gleichgewichtsbedingung (5,30) als auch zu den 
kinematischen Gleichungen (5,39) und (5,48) hergestellt ist. Durch (6,16) wird dann auch der 
Skalarteil der Spannungsfunktion Q festgelegt: in Analogie zu (5,51) wird 


1 e 
Q=Lety, (2, — GUY: 0) (6,17) 


zusammen mit (5,31) ist nunmehr &2 durch ® véllig bestimmt. Man erkennt: Um die Symmetrie- 
verhaltnisse des Momententensors beurteilen zu kénnen, miissen das Elastizitatsgesetz und eine 
Vertraglichkeitsbedingung herangezogen werden. Bei der ebenen Platte verwischen sich diese 
Zusammenhange: die Symmetrie des Momententensors kann dort allein aus der Symmetrie der 
raumlichen Schubspannungen erschlossen werden, die allerdings ihrerseits nur aus der speziellen 
Kinematik des raumlichen Punktkontinuums heraus begriindet werden kann. (In einem etwas 
anderen Zusammenhang hat D. Riidiger! die Unsymmetrie des Momententensors crértert.) Die 
tbliche Vernachlassigung 

mi! = 0 (6,18) 


entspricht der Annahme (4,52). Sie hat fiir den Skalarteil von 2 die Gleichung 
1 
QO = eB 7, Dz (6,19) 
zur Folge, in Analogie zu (5,53). 
Wir fahren in der Behandlung der Elastizitatsgleichungen fort und rechnen sie auf die Mo- 


mente M,, und die Verkriimmungen x‘! um, wobei wir uns die Einzelheiten der Rechnung ersparen 


“1 2D. Riidiger, Ing.-Arch. 28 (1959) S. 281. ~ 
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wollen. Mit den geometrischen Tensoren 
# 7 *. 
UE = Orb} 6! b; 2 bi a,,, 


il ihe Fil 
Ves =a bo oe a,5 3 (6,20) 
~ * * . _* 
Tit, bt + bai—2ath,,, 
eas a bi Le a‘! bee = ee Vai ’ 
so daB also, wenn wir symbolisch die Matrizen 
* * 
Tere 
el ae fae ps us (6,21) 
Oi Die eee 
bi 58 


der gemischten Formen des ersten Haupttensors und des Deviators des zweiten Haupttensors 
einfiihren, die Dualitaten 


Ui! (D, B) = U'!,, (B,D), Vil (D, 8) = Vil (8, D) (6,22) 


bestehen, lautet jetzt das Elastizitatsgesetz: 


a: Et me ees . ue RM 
Ni = =< la —y?) [ais al xp +12 Uap) tr(a ois V.'gprt8) , (6,23)* 
E 3 , — nw na 
M;, = Dd ea [(L—v) (ain arg %°8 + Uy,7? By 9) +9 (ayy 2 + Vi"? a 5)13 
12 fo Rae ? yi ye 
a '=— ge [a +r) (0 ae Mag tT iB Ui ap N j= (a1 M TB Fe ape ) : (6 24) * 
1 kA 
: ei pest a [1 +) (i. a1 nv’ ar Dee M,») —v (a, N+ Vine Mx~)| 
mit 4 = Gap ob, M = a*? My,. (6,25) 
Dazu kommen die Gleichgewichtsbedingungen 
Va N* — e%4 bu bi 7, M,, =0, 
e%4 ofr Va Ve Mya =F bap N= =0 ° (6,26)* 
Cap Ne + XB 4 Mg, = 0 Py 
deren Lésungen aus einer Spannungsfunktion ® abgeleitet werden kénnen: 
i F 4 8 
N" = ei e!8 [74 Vp ® + Va(b} B,)] —> ce BL 7, De, (6.20 


My =s(Vi% + Vi ®) —b, B+ HS (FV — OV) Va. 


Das unterstrichene Glied fallt weg, wenn der Momententensor als symmetrisch angenommen wird. 
Ebenso leiten sich die Deformationen, welche den Vertraglichkeitsbedingungen 
Vg %*! — 4 of BI 7gt,, =0, 
eho Tq Vp laa + bap e'? =0, (6,28)* 
Cap 5B + X86 b4 €, = 0 
gentigen, aus einem Verschiebungsvektor v ab: 
: é tea 
UL att ptB a aed tang 4 
: CA Ware al “ol =e ot DLV, 09, | (Game 
3 t 
Eil SAV str V 1%) — 5,0 ago OF V 5 OF Vi) Var- 
Das unterstrichene Glied fallt weg, wenn der Tensor ¢;, als symmetrisch angenommen wird. 


Mit seinem symmetrischen Teil und dem Momententensor ist auch der ganze Langskrafttensor 
bestimmt. Er berechnet sich zu 
ip oseek? ; x Peat gree nes. 470 
N 12 la ed) (ais alB Exp 72 Wr ap xs) +9 (ai é+ 72 Ve ap xs) 


Wi = Ul +5 (tid— 06). | (6,31) 


(6,30) 


mit 
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Eine analoge Formel kénnte auch fiir die Verkriimmungen x! aufgestellt werden; dies ist aber nicht 
erforderlich, weil das Elastizitatsgesetz hier unmittelbar den vollen Tensor liefert. 
Sehr leicht la8t sich jetzt auch die elastische Energie angeben; als Deformationsenergie 


M(é1, x *] = Qo | | os [(l—») a’ abu a3 ya%b at] bap ean df 


1 ; Ee PR ee 
rrr I] T2198) UA) ag ay HY Gap py] 29° 204 


+ 2[ —») UM, eV lag Mt df,  (6,32)* 


und als Spannungsenergie 


MMi, N"]= : | | gt tv) a%4 abu —y a8 git’) Map M,, of 


1 1 sic 
“ [J 7 UCL + Y) Gay gy —? ap a, ,] NOP 


+2[(1 +>) 0%, —v V2", Mag N’"} df. (6,33) * 
AbschlieBend wollen wir die Frage der Randbedingungen erértern. Die virtuelle Anderung der 


Formanderungsenergie ist gleich der negativen virtuellen Arbeit der BeanspruchungsgréBen, so 
da wir aus (4,2) erhalten 


Oll[v, wo] = dT [e;; (v, w), x"*(w)] = — 0AM 
= Jf e*? {Kg 0 [Va(5v) + ay X dw] + Mg o Vq(de)} df’. (6,34) 


Sind die Gleichgewichtsbedingungen fiir die Beanspruchungen erfiillt, so bleibt nur ein Rand- 
ausdruck iibrig 


dIT[v, o] = SIT [v, ]rana = f e*? (Kp 0 60 + My 06) ng ds ; (6,35) 
dlT[v, w] = p (Ky t* odv + My t* 0 dw) ds. (6,36) 
Nun ist 
cnet 
t ats (6,37) 


und wegen 
M;,oE=0, €,oE = 0 
1aBt sich (6,36) mit Hilfe von (3,12) und (4,16) in 
élT[v] = f (K o dv + M odw(v)) ds = f [K o dv — (e* o M) (V,(6v) o E)] ds (6,38) 
umformen. Ferner ist > 
eX =—n*t— tn; (6,39) 
damit ergibt sich 


‘ 0) 
dII[v] = g K + (Mp E) 0) (My B) 0 | ds (6,40) 
und schlieBlich, nach einer partiellen Integration, 
6 
dII[v] = g {|x — 2 (My B)| 2 8v—(Mp EB) ct ds (6,41)* 


mit der bereits von (4,17) her bekannten Zusatzkraft im Integranden. Nachtraglich fihren wir 
hier noch die BeanspruchungsgréBen wieder ein; es ist 


My = M on= Mygt*n’, Mz = Mot = My,t*t, (6,42) 
und ferner nach (4,5), (4,16) und (4,21) 
K = K, t® = N%n, =(N*? a; + N*E) ng. (6,43) 


Ganz entsprechend kénnen wir aus (4,34) die Randformel fiir die virtuelle Anderung der Span- 
nungsenergie (6,33) entwickeln. Das Ergebnis ist, wenn die Beanspruchungen ein Gleichgewichts- 
system bilden, 


dII[M; (®), N’* (&)] = — g {« ene Gy B) 06@ —(e E) | ds (6,44)* 


% = XK" ny = (x%? ap + 2% E) ny, (6,45) 
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und es ist nach (3,31), 3,32), (3,5) 
y =€on = Eggi? n?, e=eot=e,,0t. (6,46) 
7. Die Schalengleichungen in den Koordinaten der Kriimmungslinien; Zylinder- und Kugelschale; 
Membranspannungszustand. Wenn es irgend méglich ist, wird man als Parameterkurven auf der 
Mittelflache die Kriimmungslinien wahlen — insbesondere also dann, wenn der Schalenrand eine 
Kriimmungslinie ist oder stiickweise aus solchen besteht. Da es sich um ein Orthogonalnetz han- 
delt, kénnen wir unsere Schreibweise etwas vereinfachen; wir setzen 


ay, = %> Ang = Xe » 
und es ist 
Ries oY pees a= 4,6 (7,1) 
Bog? mer = X1 Xe ) 


Die fiir das kovariante Differentiieren benétigten Dreiindizessymbole sind 


1 Well 1 raed 1 i a 
Pus zo %1%» Di2= F gq, %» 32 = — gage 
2 bed jee | i Dae f ee 
2 2 
DiS == 2 lat» 12= 9 %, A102» Dr2= 3 |, 92% 


Im Folgenden sind weiter die geometrischen Tensoren (6,20) und (6,31) zusammengestellt, wobei 
die Vorzeichen nach (2,27) gewahlt wurden: 


wy = Uae = USA | 
11 & {1 1 
U_ 22 =a (ee). 
12 1 
U1 = Us = 0, 


aE ia 1 Mie fet: 1 
yi? =tle—-z) id eres Lees 7 
Puce 4 2 R, R, 9 U5 2 2 R, Re 9 


o \R,  R, 

12 
V Se 0: (7,3) 
22 Deda 
Vine 

a \R, Re 
722 
Fo ir, 3 = 0 
wil i Ti i 
Wee Seer? Vi = Becys 


rs Xp pid 
Oi. =i. One 
Oa, 


(entsprechend fiir die iibrigen GréBen); 


1 1 
ue pain Ww? a Ef 2 1 
aS saat 21h, Ry? 
Ni es Paes) 1 
Ww”? Sa | Vy..=— 
me ae Giae 7h OR 


sonst ist 
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Damit folgt 


we nee Ne 1 ‘ Oy e 1 1 
SS +v—e = 22 
1—r oa, a, | 21 Oh, 22 19°22 Ri n 


nS ; = (1 — v) <a ese + oy : : = (12 — 2 RT 2S (1 — v) a €a» p (U4)* 
y= is we faa a £11 501% lz— ae ; | 

M,, =— pao 0 vod — al Be | 

ae me oS 35 (L — v) oy &% x}, M3, = 12 a5" V) Oy &» x72, (7,5)* 
M,, = — vate OX |x22 4 ve mit 4 na (RE — z;) éx| . 


mit folgender Umkehrung: 


2st a te 1 1 , 
il —— eee: 1 a 22 
Et® 0,0, nae Ade Mas 12.1% (R; a) | 4 ] 
12 12 1 
1 ee a = 
ae rs Ew (1 + v) aX, & My. 9 oat *. Et (1 + ee a M2, ? (7,6) * 
12 1 e 1 1 
Dio Sa ty ee fe 11 
ze Et? a,a | Bare OR Mas + 31% R, nN |; 


1 11 Oe p22 {eesti 1 
E41 = Fam | rae +e — x) Maa 
é 1 2 


A 1 a Late (al 1 1 A 
Ee = 7; (L +4) Oy O Musee — £)(Ms—Ma)] ¥ 7,0 + 9) 24 0% N*, (7,7)* 


CPC ORY ae) 
eee 22 By nil j= 1 ) 
€9g = Fr, 2 2 N ag. N ae R, My, Zi 
Ferner ist noch 
Et Ue 7 Teel: 2 l 
Pee ears a aes \ eee eas 3 12 21 
Ton, Seger ea aera Rete ec 
(7,8) 


Et Telarc 13) Lome? 
Ng ree (iy) a pe gy gg, sm |e? 12 
(asp oat aang 8 1s le” Reo Bute 


. 


. . . . . . x. . oe ° oe 
Da wir in diesem Zusammenhang den antimetrischen Teil x‘! von x'! vernachlassigen diirfen, 
kénnen wir statt (7,8) etwas tibersichtlicher schreiben 


R Lee 1\ 39 
ER PED AG Reed Phen 


gel Bt LS A ? Teal el! AS 64 
N = E40 ele t am Flee : 


Et 
1— vr? 


ie = 


(1 —») 


Oty My 


(7,9) 


Die Gleichgewichtsbedingungen (6,26) spaltet man am besten wieder nach (4,27), (4,28) auf, 
da die Formeln sonst zu unibersichtlich werden: 


ee Ne =O 1 N, ETL AN 2N } 
pee ed oot Se TS 5Ns +5 N'=0 : 
I. 
(7,10) 
a,N? a a,N7? eS pe N22 rN 42 n"?) 


i 
+I N40 NY +N? +e N= 05 
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N' = : (0,Mo» — 0,My. + T's. My, + D2 My,—It2 M,,—Ti2 Mp») ’ 
ie (7,11) 
N* = 1 (.M,; — 2, Mg, + T'2 Myp + Dt My —T't2 Myp—T'i2 Mus) 
1% 
1 1 
OjNI a ONO (ara et alee) Was =a =0. (7,12) 


Und schlieBlich ist nach (6,29) 
1 Br fat 1 
7B me ee pt Pilg, ie a.) — 1}, 8 — dal ) 


2 


9 
a 


ee 


' i. ai 
2.00 a Its fa — dv) Ths Feely eo (0,02 — O0,) — 2 (Fy 


i 
Oy Xp 


Lael iE 
el 7 2.20 + I?, (F; V2 — ax) —I,dv— R23 (0.v; — 0,2) — 0, (F; 2] : 


1 i! 1 
1d — ave ey gels (F; 1 — dye) — Thy tae — 24 (F; *) 
o 
a1 = 4 — Th Tie + RY 
A iN : 
£49 (2, 2+ 0:%) —LT.% —L7 2%» (7,14) 


a 
= 1 ] 2 2 


Setzt man nun die durch die Verschiebungen ausgedriickten Deformationen (7,13), (7,14) in 
die Gleichungen (7,4) baw. (7,9) und (7,5) des Stoffgesetzes ein, so sind die Beanspruchungen N'! 
und M;; und nach (7,11) auch die Querkrafte als Funktionen des Verschiebungsfeldes der Mittel- 
fliche bekannt. Die Gleichgewichtsbedingungen (7,10) und (7,12) liefern dann 3 Differential- 
gleichungen fiir die Komponenten v,, v, und v des Verschiebungsvektors v. In analoger Weise ver- 
fahrt man, um die Differentialgleichungen des dualen Problems fiir die Komponenten ®,, ®, und ® 
des Spannungsfunktionenvektors ® abzuleiten; das soll hier nicht weiter ausgefiihrt werden, weil 
sich dabei weder besondere Schwierigkeiten noch neue Einsichten ergeben. 

Besonders einfach werden unsere Gleichungen fiir den geraden Kreiszylinder (Radius R), auf 
dem wir die Koordinaten q() = s, in Richtung der Erzeugenden, q®) = s,(= Rq) in Umfangs- 
richtung und qg®°) = z in Richtung der auBeren Normalen einfiihren, so daB 

= R, > oo, R,= Rh (7,15) 


wird. Samtliche Dreiindizessymbole verschwinden, und es ist nicht mehr erforderlich, zwischen 
kovarianten und kontravarianten Indizes zu unterscheiden; der Einfachheit halber schreiben wir 
sie jetzt alle als untere Indizes. Wir bekommen mit c? = ¢?/12 R? 


Et 
M1 a ee (E11 + ¥ 92 — 0? R x9) » 


Et A 1 
Nig a (be) (é + zo Rn), 


(7,16) 
Et x il 
Ney =qaal—(&e— ge Roa); 
Et 
Noo = Top (fae +¥ 8, + Rxy,)5 
Et 1 
Ma =~ ia ( 19 Maa + Re); 
Et 
M,, = — Rasa”) %129 
(7,17) 
E# 
Maa = — pq (1 — 9) ma» 


“im Ee 1 
M2. eae. (1 —??) ‘ +P yy eed 
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MN, = 0,Myo ee 02M,» ’ N, = 02M, , in 0,M,, 3 (7,18) 
1 
441 => 020.0 ae OnVq 9 ] 
ert 
%12 = — 010.0 — Roo (0,% — 0904) » 
1 
%o1 =S 0,040 ———s R 040, 9 
Yoo = 010103 (7,19) 
&1 = 01%, 


1 
Si x (01% + 090) » 
1 
€59 => O05 4 Rie . 
OM ap OgNoi =0, 
1 
0,Ni2 + 02oo “- Ri Ny =0, (7,20) 
1 
0,N, + 025 — # Nez == 0 e 


Die Differentialgleichungen fiir v,, v, und v lauten 


(22 + Pe : : a | Ro]e= 0. 
il 
ae 2}e + (a0, + ec: 


[-# Ra,(a,— ~F tata) + ale, + ( 


1 
ae A ee eee 
mit symmetrischer Matrix der Differentialoperatoren. 
Diese Gleichungen finden sich auch bei Wlassow!. Da dort mit anderen Deformationen und 
Beanspruchungsgréf8en gerechnet wird, mag diese Ubereinstimmung im Resultat als willkommene 
Kontrolle angesehen werden. Rechnet man nun das duale Problem fiir die Spannungsfunktionen 


@,, ®, und @ durch, so kommt man auf folgende SchluBgleichungen (es sind dies die in den Span- 
nungsfunktionen geschriebenen Vertraglichkeitsbedingungen) : 


(212s +3" ale) ®, + (“F"Aa)oa + [oP RA, (B10, + FAA) — eA] =O. 


(FA) Ps + (Ae? +A 77 00,) % aha me Para +22 tx a,|  =0, 


7,22 
fe Ra (a2 ese ax) = a ©, + = RG, ( +19, +2 x) > 4 ®, ae 
ee Jer(A AO. zm) +7 ®=0, 
mit ebenfalls symmetrischer Matrix der yee ialoperatoved Die Ansatze 
ee Peseta 
V, = Ay ‘R cap ape ee ie) 
Reet wr gages 


1 Siehe FuBnote 2 von Seite 160. ; 
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zur Lisung der Differentialgleichungen (7,21) bzw. (7,22) fiihren in beiden Fallen zu der »,charak- 
teristischen Gleichung*‘ é 
1—? 


A8 — 2 (2 m? — y) A8 4 + 6 m? (m? — 1) + 1/4 


C2 


— 2 m? [2 m4 + (4 — 1) m2 + 2] 2? + m4 (m? — 1)? = 0. (7,24) 


Physikalisch bedeutet dies: Verschiebungen und Spannungsfunktionen zeigen das gleiche Abkling- 
verhalten vom Rande weg, was zwar plausibel, aber den Differentialgleichungen (7,21), (7,22) 
nicht ohne weiteres anzusehen ist. Nur in unwesentlichen Termen weicht (7,24) von den charak- 
teristischen Gleichungen ab, die W. Fliigge! und C. B. Biezeno und R.Grammel* bzw. K. Girk- 
mann? gefunden haben. (Vereinfacht man das ganze Problem dadurch, daB man zum Elastizitats- 
gesetz des ebenen Flachentragers iibergeht, so erhalt man, wie H. Schaefer* gezeigt hat, Diffe- 
rentialgleichungen, die sich voneinander nur noch durch das Vorzeichen von » unterscheiden, so 
da die Analogie starker hervortritt; in den charakteristischen Gleichungen selbst, die wieder tiber- 
einstimmen, kommt y nur in der Kombination (1 — vy?) vor. 

Starke Vereinfachungen ergeben sich auch bei der Kugelschale. Ist ihr Halbmesser R, so wird 
der zweite Haupttensor 


bi) = — ‘4 ils (7,25) 


also ebenfalls kovariant konstant. Infolgedessen wird nach (6,27) der Momententensor symmetrisch 
und daher auch der Verkriimmungstensor. Aus der algebraischen Gleichgewichtsbedingung folgt 
dann ferner die Symmetrie des Langskrafttensors—kurz: alle FeldgréBen der Schale werden durch 
symmetrische Tensoren beschrieben, also 


ul Gil Et ; ee BT 
N' =WN Zee [(1 —v) at* al8 + yp at! a%?] Exp > | 
(7,26) 
M;, = M,, = ete [1 —v) ain aig + 7 aj) ay ,] 22° 5 
U il 12 (1 — v2) ta “lp il a ? 
he eH le 12 p43 ia gl it 081 
we x — Fa l + y) ai* al? —y at! @? | Men» | 
; (7,27) 
Ei = Ei = x LA + ¥) Gin 1g —Y ait Oy g] NP ; | 


ie ae “a “za 
Va( fi ‘+04 ole} ,,) = 0, e%4 eH V7, \7 5M, — x tag Ne = 0; (7,28) 


(Gleichgewicht) 
“~~ il “a gia “”~ 
Va (A! + oro B2,,] == 0), e*4 eu]. Ve ban — tap 2? = 05 (7,29) 
(Vertraglichkeit) 


Net eta ole Ve Vio a’! axe = (ai a'? + a! 7) Va PD, ? | 
. | (7,30) 


K ] 
Mu =3(Vi® + V,®) +49; 


agi eX el 8 We Ve v —— 


2 1 : . 
ai hast — —(ai® alP + gla «| Wig OE 
(7,31) 


rR 1 1 
cL (Vint Vie)+ Rar. 


1 Siehe FuBnote 3 von Seite 167. 

2 C. B. Biezeno u. R.Grammel, Technische Dynamik, Bd. 1, 2. Aufl Berli Otti i 

3 K. Girkmann, Flachentragwerke, 5. Aufl. Wien 1959. ‘ ep ee 
4 Siehe FuBnote 4 von Seite 160. 
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Artet die Kugel mit R —> co zum ebenen Flachentrager aus, so reduzieren sich die Gleichungen 
(7,26) (7,31) auf das bekannte Analogiesystem Platte/Scheibe; D(q, g®)) wird die Airysche Span- 
nungsfunktion (H. Schaefer). 


Ks ist selbstverstandlich, daB in unseren Ansatzen auch die Gleichungen des Membranspannungs- 
zustandes der Schale enthalten sein miissen. M. Lagally? hat der Frage des Membranspannungs- 
zustandes eine interessante Abhandlung gewidmet; eine nicht kovariante Herleitung aus den allge- 
meinen Schalengleichungen gab der Verfasser®. Hier soll eine kovariante Darstellung gegeben 
werden: Im momentenfreien Zustand der Schale verschwinden nach (4,29) auch die Querkrafte N’, 


so da nach (4,24) 
bag 28 = —V, 2 (7,32) 


wird. Wir nehmen jetzt an, daB die Gaufsche Kriimmung K der Schale nicht verschwindet. Dann 
kénnen wir den symmetrischen Tensor 


Piss aie & es (7,33) 
bilden, dessen Matrix die Reziproke zur Matrix des Tensors b;, ist: 
Bix b,, = ef4 ere. ai Big = (ai* at — ait a?) = by, Bia 
=¢(2Hbj— cl) = Kj = 5) (7,34) 
nach (2,30). Wir iiberschieben (7,32) mit B®! und beachten (5,34); dann erhalten wir 
Ol wh ue bie Va 2 (7,35) 


Nun ist aber, der algebraischen Gleichgewichtsbedingung (4,38) zufolge, im momentenfreien Zustand 
der Langskrafttensor symmetrisch: 


exp NXP = 0, (7,36) 


- und hieraus folgt nach (4,24) 
Vp 2° —2HQ=0. (7,37) 


So ergibt sich, wenn man noch die Codazaischen Gleichungen (2,33) benutzt, fiir Q(q(), q@)) die 
Differentialgleichung zweiter Ordnung 


o>4 oft B, Valk Vp) +2HO=0. (7,38) 


Sie wurde von Lagally als die ,,charakteristische Differentialgleichung‘‘ des Problems bezeichnet. 
Wir zeigen noch, daB die Gleichgewichtsbedingung 


Ba,iN' == 0 (7,39) 
eine Identitat wird; aus ihr folgt namlich der Reihe nach 
bg e%4 (VW, O° — bf Q) =0, 
bap eX4 Va OF —et4G,, Q = by, e*4 Va Qe == (} 9 


1 \ 
ba 5 Sa co ve a (7,40) 


1 
ory, Fac Tg Seer: a| = 0) 
e%4 Vale K tap ve jae Vor 8 
und die letzte Gleichung ist die gesuchte Identitat. 


1 Siehe FuBnote 1 von Seite 160. 
2 M. Lagally, Z. angew. Math. Mech. 4 (1924) S. 377. 
3 Siehe FuBnote 1 von Seite 170. 
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Alles, was bisher iiber den Membranspannungszustand gesagt wurde, kann natiirlich auf Grund 
unserer Analogie wortlich auf den Fall der ,,reinen Biegung“ iibertragen werden, in dem die Mittel- 
flache von Verzerrungen frei ist — ein Fall, der besonders in den Anfangen der Schalentheorie gern 
betrachtet wurde, als die differentialgeometrischen Darstellungsméglichkeiten noch nicht aus- 
reichten, um den allgemeinen Fall in einigermaBen iibersichtlicher Form behandeln zu kénnen. 

Es bleibt noch die Frage zu erértern, ob der Membranspannungszustand, der ja innerlich 
statisch bestimmt ist, bei vorgegebenem Stoffgesetz auch kinematisch méglich ist. Es ist bekannt, 
daB dies zwar allgemein nicht zutrifft, jedenfalls aber bei der Kugelschale konstanter Wandstarke, 
deren Deformationen nach (4,27) mit den Beanspruchungen verkniipft sind. Da die Momente 
M;, verschwinden, miissen auch die x‘! und nach (3,51) auch die x' Null sein, was nach (3,48) 


wo=0, w=0 (7,41) 
zur Folge hat. Wegen ¢; = 0 gibt dann die zweite der Gleichungen (3,49) 
v; — R Was v 5 (7,42) 
so daB der Tensor der Membrandeformationen 
. 1 
ey(— RVs Viv + atu?) (7,43) 


wird, also symmetrisch ist, wie es sein muB. Der Deformationszustand ist also durch die Normal- 
verschiebung v allein bestimmt. Driickt man jetzt gema® (7,26) den Langskrafttensor N‘' durch é,, 
und damit durch v aus, so erhalt man aus den allein noch tibrigbleibenden Gleichgewichtsbedin- 
gungen 


Wi Kea (7,44) 


zwei Differentialgleichungen fiir v. Wir unterdriicken die umstandliche Rechnung und schreiben 
gleich das Ergebnis hin: 


Vi(A + 3°) =9, (7,45) 
woraus 


2 
A v + xv = konst. (7,46) 


folgt, und dies ist die gesuchte Differentialgleichung, waihrend die charakteristische Differential- 
gleichung der Kugel 


AQ+22=0 (7,47) 


lautet. 


(Eingegangen am 4, Mai 1960.) 


Anschrift des Verfassers: Prof. Dr. W. Giinther, Karlsruhe, Technische Hochschule, 
Lehrstuhl fiir Theoretische Mechanik. 
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Ausbildung eines Wirbelpaares an den Kanten einer Platte 


Von E. Wedemeyer 


1. Einleitung. Bei stationarer Potentialstrémung um eine Kante treten an ihr unendlich groBe 
Geschwindigkeiten und damit negativ unendlich grofe Driicke auf, die bei einer realen Fliissigkeit 
unmoéglich sind. Wie die Erfahrung zeigt, lésen sich stattdessen von der Kante Wirbel ab, deren 
Einflu8 auf die Strémung in jedem Augenblick derart ist, daB die Geschwindigkeit an der Kante 
endlich bleibt. Bei einer Behandlung des Problems fiir ideale Fliissigkeiten sind die Wirbel auf 
eine Flache konzentriert, die von der Kante auslauft und sich spiralig aufrollt. Die Tangential- 
komponente der Geschwindigkeit (v,) macht einen Sprung in der Wirbelflache (auch Unstetigkeits- 
flache) vom Betrag y = dI’/ds = Belegungsdichte der Zirkulation. 

Die vorliegende Rechnung behandelt die Ausbildung eines symmetrischen Paares von Unstetig- 
keitsflachen an den Kanten einer senkrecht angestrémten Platte. Fiir die Unstetigkeitsflache 
an der Kante einer einseitig unendlich breiten Platte liegt bereits eine von Anton! gefundene 
Lésung vor. Im ersten Teil der vorliegenden Arbeit wird die von Anton! behandelte Aufgabe noch 
einmal mit einem etwas geanderten Verfahren behandelt, das gré®ere Genauigkeit erméglicht und 
als Ausgang fiir die Wirbelentwicklung an der endlich breiten Platte dienen konnte. 

Als Nullpunkt der Strémungsebene (z-Ebene) wird die eine Plattenkante gewahlt und die Platte 
in die negative y-Achse verlegt. Bei einer Breite der Platte 2 H ist also ihr Bild gegeben durch das 
Stiick —2 H <y < 0. (Abb. 10) Ist die senkrecht zur Platte gerichteteAnstrémgeschwindigkeit imUn- 
endlichen v,9, so ist im ersten Augenblick, solange sich noch kein Wirbel von merklicher GréBe aus- 
gebildet hat, das komplexe Potential der Strémung F =v, /2iHz-+ 2 und die konjugierte 

komplexe Geschwindigkeit v = v, —iv, = dF/dz. 
; La®t man nun H ->oo und v,,—0 gehen darart, daB vy /H = konst. = , /H, ist, so laBt 
sich das Potential fiir die unendlich breite Platte auch schreiben F = v, y2 tH, z. Dabei ist v, der 


Betrag der Geschwindigkeit im Punkt z, = —i H,/2 und H, eine willkiirlich festgelegte Lange. Die 
unendlich breite Platte wird jetzt durch das Stiick — oo < y < 0 der y-Achse 
dargestellt ; die Plattenkante, von der die spiralférmige Unstetigkeitsflache aus- ly 
geht, liegt im Punkt z = 0 (Abb. 1). 
Das Problem fiir die unendlich breite Platte ist insofern sehr viel einfacher, v 


als die zeitlicheAnderung von Form und Belegungsdichte der Unstetigkeitsflache 

in einer 4hnlichen VergréBerung erfolgt. Bei der Platte endlicher Breite stimmt 

der Vorgang am Anfang, solange der sich bildende Wirbel noch klein gegentiber 

der Plattenbreite ist, mit dem an der unendlich breiten Platte iberein. Daher 

~ kann man bei der endlich breiten Platte von der Lésung fiir die unendlich breite Prk Ge Pearle 

Platte ausgehen und die Entwicklung von Form und Belegung der Unstetigkeits- Unstetigheitefliche an 

flache fiir eine Reihe spaterer Zeitpunkte verfolgen. eh Sree peta 
Ist P ein Punkt auf der Unstetigkeitsfliche, so miissen in P dauernd fol- 


gende Bedingungen erfiillt sein: 

1. Die Normalkomponente der Strémungsgeschwindigkeit v, ist zu beiden Seiten der Flache 
gleich und gleich der Geschwindigkeit V,,, mit der sich die Unstetigkeitsflache bei P senkrecht zu 
sich selbst verschiebt: 

0, Ve (1) 


2. Der Druck ist zu beiden Seiten der Unstetigkeitsflache gleich. Bezeichne der Index I die 
linke, der Index II die rechte Seite zur Durchlaufungsrichtung, so ist infolge der Bernoullischen 
Gleichung fir instationare Strémung 


1 a 
Gx (% — %) = aL (®; — Py) 


1 [,, Anton, Ing.-Arch. 10 (1939) S. 411. 
13* 
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und da vy n = Uitn ist, so folgt 
1 2) 
cy (vie =e v1) (vn: — v1) = at (2, — @,,) . 


Nun ist 0, —@,, = Ip die Zirkulation um das durch P abgetrennte »freie™* Ende der Unstetig- 
keitsfliche. Bezeichnet man mit v, den Mittelwert der tangentialen Geschwindigkeit und mit y 
die Belegungsdichte der Zirkulation, so ist 


bg 
m4 =4—F, Ory == Octo 
und man hat 
oly 
Yee 


wobei der Differentialquotient 0©/dt und damit auch 0/',/0t an festgehaltener Stelle 2 zu nehmen 
ist. Das bedeutet, daB Ay = y v, At den Unterschied der Zirkulation zweier solcher Abschnitte 
darstellt, die durch Verschiebung von P senkrecht zur Unstetigkeitsflache nach der Zeit At aus- 
einander hervorgehen. Lat man dagegen, wie es bei ahnlich wachsender Unstetigkeitsflache 
zweckmafig ist, den Schnittpunkt P auSerdem in tangentialer Richtung mit der der abnlichen 
VergréBerung entsprechenden Geschwindigkeit V, wandern und vergleicht nunmehr die Zirku- 
lationen J",(t) und J',(t + At) zweier durch P abgetrennter Abschnitte, so findet man 


P(t + At) —I'p(t) = AD y = y (uy, — V;) At 


oder 
al’ 
y (4 — VY) = SF (2) 
t 
(siehe Abb. 2). 
uy 
Plt+At) E 
Plt Vy At 
Vy -A t 
Abb, 2. Skizze zur anschaulichen Herleitung von Formel 2. Abb. 3. Spiralférmige Unstetigkeitsfliche in der ¢-Ebene. 


2. Theorie der unendlich breiten Platte. — a) Ahnlichkeitsgesetz. Durch die konforme Ab- 
bildung € =& + in =, j2 i z/H, wird die z-Ebene in die positiv imaginare Halfte der C-Ebene 
abgebildet. Das Bild der unendlich breiten Platte geht dabei iiber in die reelle Achse der ¢-Ebene 
(yn = 0), wobei das Stiick € < 0 die Vorderseite, & > 0 die Riickseite der Platte darstellt. Der 
Punkt —i H,/2 in der z-Ebene, in dem die Geschwindigkeit -{ v, herrscht, entspricht den willkirlich 
wahlbaren Punkten + €, der reellen Achse der ¢-Ebene. Der Einfachheit halber sei im folgenden 
&, = H, gewahlt. Damit wird die Abbildungsfunktion 


¢=)2i Hz (3) 
und 
ne: ae pee 
TER ale «(8 om aap PaO a (4) 
Zur Berechnung der Geschwindigkeiten in der €-Ebene werde das Bild der Unstetigkeitsflache an 
der £-Achse gespiegelt (Abb. 3). Hierbei bleibt letztere Stromlinie, wenn man zudem y(0) = — ye(C) 


macht. Hier bedeutet y, die Belegungsdichte in der ¢-Ebene, die mit der in der z-Ebene (y) zu- 
sammenhangt durch y = y,|d¢/dz|. Die anfangliche Geschwindigkeitsverteilung in der z-Ebene 
Vy) = dF [dz = v, Vi H,/2z geht in der ¢-Ebene in Wy = Wo = V,/(dl/dz) = v,, also in eine konstante 
Parallelstrémung in Richtung der reellen Achse iiber. Dariiber itberlagert sich nach Ausbildung 
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des Wirbels als dessen Feld eine Stérgeschwindigkeit. Ist die resultierende Geschwindigkeit w, so 
ergibt sich fiir w 


++ 00 
= A 1 
W = w,— iw, = 0, + 5 [ tee (5) 


Dabei ist s die Bogenlange auf der Unstetigkeitsflache in der ¢-Ebene, und es gilt s > 0 fiir die 
Spiralwindung im positiv Imagindren und s < 0 fiir das Spiegelbild im negativ Imaginaren. Die 
Geschwindigkeit v in der z-Ebene ergibt sich hieraus durch die Beziehung v = w dé/dz. 

Da bei der unendlich breiten Platte eine Vergleichslinge fehlt, andern sich Form und Belegungs- 
dichte der Unstetigkeitsflache, wie schon erwahnt, in Form einer ahnlichen VergréBerung, und 
zwar, wie schon vor langer Zeit Prandtl* ganz allgemein gezeigt hat, nach einem Potenzgesetz. Die 
von Anton aufgestellten Ahnlichkeitsgesetze mégen hier noch einmal kurz abgeleitet werden. 

Ist 2“) ein Punkt der Unstetigkeitsfliche zur Zeit «, F® das komplexe Potential in diesem 
Punkt zur Zeit ¢), so hat man fiir einen zu 2) dhnlichen Punkt z zur Zeit t anzusetzen 2/20) 
= (z/t)" sowie fir das komplexe Potential in diesem Punkt F/FO = (1/)". DemgemaB gilt 
dann fiir die Strémungsgeschwindigkeiten 


v ee) v ta. Un UE we dF/dz ae t \m—n 
OO 40 0 ar Oa =(q5) 


Die Zirkulation [‘, um das durch den Punkt z abgetrennte freie Ende der Wirbelflache verhalt 
sich als Potentialdifferenz auf beiden Seiten der Wirbelflache ebenso wie die Potentiale. Es ist also 


TP, (TY = (tf). 


Bezeichnet man die Bogenlange der Unstetigkeitsflache in der z-Ebene mit s,, so ist die Wirbel- 
belegung y = dI'/ds,. Sie ist die Geschwindigkeitsdifferenz auf beiden Seiten der Wirbelflache, 
verhalt sich also wie die Geschwindigkeiten 
= rie = (Ee 
Fur die Wachstumsgeschwindigkeit V = dz/dt bei der ahnlichen VergréBerung gilt 
. PRS Vas. Viera 
VQ VO yO fan 


== (e/))s—) 


Aus der Bedingung (1) v, = V,, ergibt sich nun 


oS ee a 
(io) “=(o) 


Diese Gleichung kann nur bestehen, wenn m = 2 n — 1 ist. Das gleiche Ergebnis liefert auch die 
Bedingung (2) 
y(v, — V;) = dI'p|de 


Fiir die €-Ebene gilt nun auf Grund der Abbildungsfunktion (3) 
j ‘ E 
£/CM = (¢/e) "2 F/FQ) == (t/t) - ia = = (t/e))m—n/2 ; 


Da aber hier die Grundgeschwindigkeit w, = v, zeitlich konstant ist, erfordert die Ahnlichkeit der 
resultierenden Strémung, da auch die Stérgeschwindigkeit (w— w,) in ahnlichen Punkten zeitlich 
konstant ist. Es mu8 demnach 
w/w) = (t/t) ™—"/? eles 
also m = n/2 sein. Zusammen mit der obigen Forderung m = 2 n — | ergibt sich 
m=1/3,- n=2/3 
und damit fiir jeweils ahnliche Punkte 
2/2) = (O28, C/O =f), FIFO STIL = (ef) 
in Ubereinstimmung mit Anton. 


Jber di i i iissigkei a Gebiet der 
Delis dil, Uber die Entstehung von Wirbeln in der idealen Fliissigkeit. Vortrage aus dem 
Hydro. A acadm ance. eae 1922, Berlin 1924. (Herausgegeben von Th. v. Karman und T. Levi- 


Civita), S. 19—34. 
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Da die Stérungsrechnungen zum grofen Teil in der ¢-Ebene durchzufiihren sind, ist es zweck- 
mifig, die Bedingungen (1) und (2) durch die in der ¢-Ebene gegebenen Groen auszudriicken. Fir 
die Geschwindigkeit der ahnlichen VergréBerung in der ¢-Ebene gilt 


ager 


i= ae nae 


mit den Komponenten normal und tangential im Punkt ¢ 


__ $n Bs 
Va W,= = 


Dabei bedeuten C, und ¢, die Komponenten des Vektors ¢ zu einem Punkt der Wirbelflache normal 
und parailel zur Tangente an die Wirbelflache im Punkt ¢. Weiterhin ist 


i dz\ |dz| /df Gy dz| Cn 
ere (ze), ete Gai ae ies 
und 
dé 
Uv, =U, ae . 
Damit nimmt die Bedingung (1) v, = V, die Form an 
dtl? on 
n dz mF, 3t °. (6) 
Da weiterhin, wenn y die Belegungsdichte der Zirkulation in der z-Ebene und y, die in der ¢-Ebene 
ist, 
or ahd pee _ [ds\ _ |dz| (de Tea 
Veh aah? cf Saeed 9S a v= (a) = la (ae) we 
ist, nimmt die Bedingung (2) y(v, — V,) = dl’,/dt die Form an 
ENE 1G ts 
ve (m di Syl see (7) 


Es kommt nun darauf an, Form und Belegungsdichte der Unstetigkeitsflache zu einer beliebigen 
Zeit so zu bestimmen, da (6) und (7) erfillt sind und daB ferner w(0) = 0, wenn anders die Ge- 
schwindigkeit an der Plattenkante (in der z-Ebene) endlich bleiben soll. 


b) Genaherte Berechnung der Unstetigkeitsflache. Wie Kaden} und fiir den vorliegen- 
den Fall Anton gezeigt haben, nehmen die inneren Spiralwindungen der Unstetigkeitsflache immer 
mehr die Form konzentrischer Kreise um den Spiralmittelpunkt an, und zwar gilt mit Polar- 
koordinaten (r,q~) um den Mittelpunkt der Spirale in der z-Ebene) 


2/3 = 
re (Eye F=2xyr 
up 


mit noch naher zu bestimmender Konstante x. 

Die Wirkung des Spiralkerns mit der Zirkulation J’; auf die duBeren Teile der Windung kann 
deshalb naherungsweise ersetzt werden durch die Wirkung eines punktférmigen Wirbels mit der 
gleichen Zirkulation im Mittelpunkt der Spirale (z,). In diesem Sinne wird die wirkliche Unstetig- 
keitsflache ersetzt durch eine auBere Windung (etwa durch den ersten Umlauf der Spirale gegeben) 
und einen Wirbelpunkt der Zirkulation J’, an der Stelle z, (bzw. Cs in der ¢-Ebene), und es wird 
verlangt, daB die Bedingungen (6) und (7) auf der AuBenwindung erfiillt sind, wahrend fiir den 
Wirbelpunkt gelten soll 


tS) a? 
d.h. die Geschwindigkeit, mit der sich der Mittelpunkt der Spirale auf Grund der Ahnlichkeitsgesetze 


bewegt, mu iibereinstimmen mit der Strémungsgeschwindigkeit an dieser Stelle. In Abb. 4 und 5 
sind schematisch AuSenwindung und Lage des Zentralwirbels (I's) in der z- und ¢-Ebene dargestellt. 


1 H. Kaden, Ing.-Arch. 2 (1931) S. 140. 
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Im folgenden wurde ein Naherungsverfahren angewandt, das von dem Antonschen verschieden 
ist und deshalb ausfiihrlicher beschrieben werden soll. 

In einer nullten, groben Naherung wird das von der Wirbelflache erzeugte Geschwindigkeitsfeld 
durch das Geschwindigkeitsfeld eines einzelnen Punktwirbels angenahert, der im Mittelpunkt der 
Spirale zu denken ist. In der ¢-Ebene hat man dann ein Geschwindigkeitsfeld, das sich zusammen- 


wy v7] : 
7) S 
x E =a E 
— 
Abb. 4. AuBenwindung und Zentral- Abb. 5. Au®enwindung und Zentral- Abb. 6. Zentralwirbel und Spiegelbild 
wirbel in der z-Ebene. wirbel in der ¢-Ebene. in der €-Ebene. 


setzt aus der Grundstrémung (w: = v,, w, = 0) und den durch einen Potentialwirbel und dessen 
Spiegelbild erzeugten Geschwindigkeiten (Abb. 6). 
Das komplexe Potential des Potentialwirbels ist nun 


und mit der Abbildungsfunktion (3) ¢ = /2iH,z, 6, = y2iH, zs 


220 t 


F = In (2 i 2/H, —)2ia/H,) = 3 lm 2 if | —In (/2iz/H, +/2i s/Mh,) ; 


woraus sich fiir die Geschwindigkeit ergibt 


a2) = = = i. —{ ik ) i |. 


z 2nilz—2sg \2iz+ 2izs/ y2iz 


Wahrend der erste Term (I°./2 2 i - 1/(z — zs)) dem Wirbel (im Punkt zg) keine Geschwindigkeit 


_ erteilt, liefert der zweite, mit z —> 2, 


come 
25a 4 rere 2.202. 


Dazu kommt noch die durch die Grundstr6mung und das Spiegelbild erzeugte Geschwindigkeit. 


Diese ist in der 6-Ebene 


cae igh 
CS 92 1g’ 
. in der z-Ebene 
Seater 
v=w os = V, sf ai 
dz ne 2% 2s} Cs 
Insgesamt ist also 
é iPS 1 
see Sh v+ —-| H, 
5 wae 2227s} Cs 


— 


bzw. 


f foo] 1 tice 1 
s 
Nun soll sein (zs) = dz,/dt, woraus fiir d,/dt folgt 


d dt\ d de Ae Sree san rai) tae| 
(2) SE - (Ee) =| ra oer + (*+ra8%, esP| 


Die Geschwindigkeit im Nullpunkt der ¢-Ebene ist 


Dielice 
140) = (+ sete 
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Damit diese Null wird, mu sein ; 
=. Pega 2 (8) 


Oben eingesetzt erhalt man 


dis _[_ tbs a oe 2 
dt =e ies? +(1 in) [ese | 2 Mr 


Die Lésung dieser Gleichung mit der Anfangsbedingung C, = 0 fiir t = 0 ist 
3 v, t\1/3 3 
f= H(z ez) a+i) 
oder 


en, Ge) * (0,72 + 0,72 i). 
1 


Fiir die Lage des Punktwirbels in der z-Ebene, die auf Grund Gleichung (3) zg = ¢/2 + Hy 
ist, und fiir dessen Zirkulation, die nach (8) —I°,/2 a = (|fs|?/2 ys) v, ist, erhalt man hieraus 


a: 3 u, 02/3 v, t\2/3 Pass 3 ys vy 11/3 
5 =H (5 i] = 1 (FF) 0,52, 7 =1Alg ee) =aAh(g) -87- 


Diese nullte Naherung ist bereits ziemlich gut. In dritter und letzter Naherung ergibt sich z. B. 
fiir die Lage des Zentralwirbels in der ¢-Ebene 


oS ss ‘aM (0,73 + i+ 0,56) 
Nea 


sowie fiir die Gesamtzirkulation der Unstetigkeitsflache 


aS. (i y 50/6358 


220,H, \H, 


Das durch die Grundstrémung v, sowie durch den Punktwirbel und dessen Spiegelbild erzeugte 
Geschwindigkeitsfeld in der ¢-Ebene ist 


a 1 A 1 1 
Pt ies Ee) 


Dieses Geschwindigkeitsfeld der nullten Naherung wird nun in erster Naherung beibehalten 
und eine Unstetigkeitsflache zur Zeit t = 8/3 (H,/v,) konstruiert, auf der die Bedingung (6) 
w,, |dC/dz? =¢,,/3 t mit einer Normalgeschwindigkeit w, aus der nullten Naherung iiberall erfiillt ist. 
‘ ea Zeit t = 8/3 (H,/v,) ist Cg/H, = (1 + i) und —I’,,/2 7 =v, H, und das Geschwindigkeits- 

e 


w 1 1 
—_—- = 1 So tar 2 
noes) Fd 
Nun ist 
dé dn dé dy 
ripe ide rt ds? > Nae a 
Aus 
a Sas SY 
"| dz 36058 
und 
dt _ |i 
rp is ae 
folgt 
ea ee 
v, | ¢ Sarit Hy ° (9) 
Die mit v, bzw. H, dimensionslos gemachten GréBen 
es: ele or: sowie é I ke 
Uy ct Ure = Nee Hy, A, i, 


\ 
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seien im folgenden der Einfachheit halber mit 


a ae ee Ww oak @ é’ Oe) a 
aS == i — oS Wie = = = ~——_ 
g n v1 Lay) H, 
bezeichnet. Desgleichen sei 
tee Ww , W. e/ Gs / @ 
Ww, = —2 = sik — —n — L — : 
v1 ° W; v, ? Gn H, ? Cr a H, und a He ~ 
Gleichung (9) lautet dann 
Se: eee lovey 
n [’/? — 8 Gn 
oder mit 
w, =w, e peal 
GE 5 ds’ 
; , dé’ dn’ 
d Lit 4 
und ¢,, W ds’ § ds’ 
, d&' , dn’ i i i , dé’ , dn’ 
(w; ds’ ee 7] E/2 1 y7?2 8 (n ds’ g a] ; (10) 


Ferner gilt 


idé'\2 dn’\2 

(3) | (a) pei 
Dies sind zwei Differentialgleichungen fiir die dimensionslosen Koordinaten é’(s’), 7/(s’) der AuBen- 
windung. Die Gleichungen wurden mit der Anfangsbedingung é’(0) = 7/(0) = 0 im Bereich der 
AuBenwindung graphisch gelést. 

Nun bestimmt man fiir die so gewonnene AuSenwindung in erster Naherung eine Belegungs- 

dichte y, aus der Bedingung (7) 
de 
dz 


e Ct Ip 
v2 (o SCs 3t 


mit einem w, aus der nullten Naherung. J’, = /';(s) war erklart als Zirkulation um das freie Ende 
der Unstetigkeitsflache 


(rls) = fre ds). 


Wenn nun letztere, im Sinne der gemachten Naherungsannahme, zerlegt ist in eine AuBenwindung 
(deren Lange s, sei) und einen Zentralwirbel, so bedeutet J", jetzt natiirlich die Zirkulation um das 
freie Ende der AuBenwindung plus die Zirkulation J’, des Zentralwirbels 


(0 — Ire ds + i) ° 
Sei ferner 


P=P(s) =fyds 


0 
und J',, die Gesamtzirkulation; so ist also 


aT yl oder Do ieee 
Nun erhalt man aus (7) mit I’, = (I, — J’) sowie mit y, = dI’/ds 


ar 
ToT +310 |dtids?—C, | ds In (L—L}1'e0). 
Fur 
=3> 


194 E. Wedemeyer: Ausbildung eines Wirbelpaares an den Kanten einer Platte —_Ingenieur-Archiv 


Dieses integriert ergibt 


Daraus folgt 


LNG Digg tees : (12) 


ZweckmaBigerweise fiihrt man hier eine dimensionslose Zirkulation 


; ip 
[= v, A, 
ein. Ebenfalls sei 
Do ee aa 
v, A, es Pico 


sowie die dimensionslose Belegungsdichte der Zirkulation 
i (je ieee 
ve vy «9, ds as” 


Gleichung (12) lautet dann auch 


=U] , (13) 


Si x. 14 
Or (14) 
[C’/2 >t 
Die Zirkulation des Zentralwirbels ist [’; = [’,, — I(s,) oder in dimensionsloser Form 
9 
— ds’ 
eet 
I's , , ae ek 
aH Es = x 1 CS a (15) 


Mit; 55 ==.Sp)/ tae 

Nach (14) und (15) kann die Belegungsdichte der AuBenwindung und die Zirkulation des Zen- 
tralwirbels bis auf den Faktor /’%, berechnet werden. Damit kénnen auch die von der AuBenwin- 
dung und dem Zentralwirbel sowie deren Spiegelbildern erzeugten Stérgeschwindigkeiten in der 
¢’-Ebene bis auf den Faktor J’;, berechnet werden. Fiir die Gesamtgeschwindigkeit, die sich aus 


der Grundstrémung und den Stérgeschwindigkeiten zusammensetzt, erhalt man in dimensionsloser 
Darstellung 


56 ~ 
LP ee Sree Pan be Pr! ad 1 GAS) ae 
V1 == \eaaek T taille a Ini (ae ds . (16) 
~y 
Dabei ist y¢(— s‘) = — yz(s’). Die Werte yz(s’) und I's sind durch (14) und (15) bis auf den 


Faktor I’, gegeben. Der Faktor [°;,, d.h. die Gesamtzirkulation, wird nun bestimmt durch die 
Forderung w’(0) = 0, d.h. durch die Forderung, daB die Geschwindigkeit an der Plattenkante 
(in der z-Ebene) endlich bleiben soll. Damit ist endgiiltig Form und Belegungsdichte der AuBen- 
windung sowie Lage und Starke des Zentralwirbels in erster Naherung festgelegt. 
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Fiir die Berechnung der zweiten Naherung wird jetzt das von der ersten Naherung erzeugte 
Geschwindigkeitsfeld (16) anstelle des vom Zentralwirbel allein erzeugten verwandt und mit (10) 
und (11) eine neue AuBenwindung konstruiert. Die neue Lage des Zentralwirbels €, (baw. zg in der 
z-Ebene ergibt sich wieder aus der Forderung 


dzs 2 2g 
EES dk St 


mit einer Strémungsgeschwindigkeit v(z,) aus der ersten Naherung. 

Danach wird, ebenso wie bei der ersten Naherung, die Belegungsdichte der AufSenwindung und 
die Zirkulation des Zentralwirbels in zweiter Naherung bis auf den Faktor J°3, sowie J’; aus der 
Forderung w(0) = 0 bestimmt usw. 

Die dritte Naherung, die im Rahmen der Zeichengenauigkeit fast vollkommen mit der zweiten 


Naherung iibereinstimmt, wird endgiiltig beibehalten und noch in die z-Ebene iibertragen auf 
Grund der Beziehungen 


A Soe de 
6 = 2 Hees Veto ae 


Abb. 7 zeigt die erste bis dritte Naherung der Unstetigkeitsflache in der ¢-Ebene. Die Abb. 8 und 9 
zeigen Form und Wirbelbelegung der ahnlich veranderlichen Unstetigkeitsflache in der z-Ebene 
zur Zeit v,-t/H, = 1. Fir beliebige Zeiten ergeben sich Form und Belegungsdichte sowie Gesamt- 
zirkulation der Unstetigkeitsflache aus den Ahnlichkeitsgesetzen 


z  [/t\23 pe ef tie ae tie 
ct =) d V1 Sh) Tt : ‘ 
3. Theorie der endlich breiten Platte. Die Anfangsstromung um die Platte der Breite 2 H war 


durch das Potential F = v,, y2 tHz + 27 mit v5 /H =D; /H, gegeben. In der Umgebung der 
oberen Plattenkante, d. h. fiir |2| < _H stimmt diese Strémung mit der an der Kante der unendlich 


ah 
i 


za 
oe iy 
HH, > 


0 LT Ae ae a ae OS Vga pemeig OF 


Abb. 7. Ahnlich veranderliche Unstetigkeitsfliche in erster, zweiter Abb. 8. Ahnlich veranderliche Unstetigkeitsflache in der z-Ebene 
und dritter Naherung in der ¢-Ebene (v,¢/H, = 8/3) fir v,t/H, = 1. 


zweite Nahg. | dritte Nahg. 


nullte Nahg. | erste Nahg. 


227755 
v, H, 2x 


1,00 | 0,83 | 0,852 | 0,867 

breiten Platte iiberein. Aus diesem Grund stimmen auch Form und Belegung der von der Kante 
auslaufenden Unstetigkeitsflache in beiden Fallen tiberein, solange die Wirbel noch sehr klein 
gegeniiber der Plattenbreite sind. Die Ubereinstimmung der beiden Strémungsformen hangt nun 
offenbar nur davon ab, daB vx, / H =v, /H, ist, wobei die Lange H, aut der unendlich breiten Platte 
willkiirlich vorgegeben werden kann. Zur Vereinfachung der Schreibweise sei deshalb im folgenden 
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H, = H und v, = vq gesetzt, ohne daf dies eine Einschrankung bedeutet. Die Platte, die durch 
das Stiick — 2 H< y < 0 der y-Achse dargestellt wird (Abb. 10), wird durch die konforme Abbil- 
dung €=£+in=é, y2 i 2/H + 2?/H? auf das Stick —&&)< g< + &) der §-Achse abgebildet. !) 
Die Symmetrieachse der Strémung, namlich die Gerade z = —1t Het x, geht in die Abschnitte ) 
&< —é, fir x< 0 und &>+&, firx>0 iiber. Zwei Punkte, die zur Achse z = —iH + % | 
symmetrisch liegen, gehen in solche tiber, die zur §-Achse symmetrisch sind, d. h. wenn das Stro-  |fj 
mungsfeld in der ¢-Ebene zur &-Achse symmetrisch ist, so ist es in der z-Ebene zur Achse 


+ wy 
‘8 


Zz 
Be 


2H— 


Abb. 10. Spiralformige Unstetigkeitsflache an den 
Kanten der endlich breiten Platte. 


OY 
Abb. 9. Belegungsdichte der Zirkulation y/v, als 
Funktion der Bogenlange s/H, in der z-Ebene fiir 
0, t/H, = 1. 
va — (1/v, - Hy) Poo = 3,95 (Gesamtzirkulation), 
ee — (1/v, - H,) Ig = 2,70 (Zirkulation des Zen- 
oO OF 08 42 16 60 tralwirbels). 


=—tH-+ x symmetrisch. Die Bilder der Wirbelflachen, die von den Kanten z=0 und 
z = — 2iH ausgehen, liegen also jetzt symmetrisch zur €-Achse und entsprechen dem Bild der 
einen Wirbelflache der unendlich breiten Platte und deren Spiegelbild (Abb. 3). Die an sich will- 
kirlich wahlbare Lange £) wird nun zweckmaBigerweise so bestimmt, daB der Abbildungsmafistab 
in der Nahe der Plattenkanten der gleiche ist, wie bei der unendlich breiten Platte. Dies ist der 
Fall, wenn €) = H ist. Die Abbildungsfunktion wird damit 


; = aces (17) 


Bei Verwendung dieses Abbildungsmafstabes erhalt man fiir den Anfang der Bewegung in der 
¢-Ebene das gleiche Bild fiir die Unstetigkeitsflache und das gleiche Strémungsbild wie bei der 
unendlich breiten Platte. Insbesondere erhalt man fiir die Anstrémgeschwindigkeit in der C-Ebene 


dF 


e = Vo0 — Vy > 
wie bei der unendlich breiten Platte. 


Das zeitliche Anwachsen der Wirbel erfolgt wie frither nach den Gleichungen 
| dC |2 dl dl 
Se ia ; Yew, 


dz dz 
(OL /ot jetzt aber genommen an Punkten, die durch senkrechte Verschiebung auseinander hervor- 
gehen; w,, w,, ye sind die Geschwindigkeiten bzw. Wirbelbelegung in der £-Ebene). Nur ist jetat 
mit € =/2iHz+2 


3 ol 


“ot ? 


n 


2 
=i. : 
was fiir kleine ¢-Werte tibergeht in i H/C. 
Man bekommt also zunachst eine Abweichung vom dhnlichen Wachstum dadurch, da® dé [dz | 
mit wachsendem ¢ mehr und mehr von i H/¢ in 1 — (H/C)? ttbergeht. Dann wird man als Folge | 


hiervon auch Abweichungen von w,, w, erhalten. 


\i 


NY 
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; Fir die Weiterverfolgung der Lésung bei der endlich breiten Platte wurde die Unstetigkeits- 
flache fiir Uco t/H = 0,1 (zu welcher Zeit die Abweichungen vom ahnlichen Wachstum noch duferst 
gering sind) als Ausgangszustand genommen. 

Die Gestalt der Unstetigkeitsflache zu einer Zeit t + At (mit nicht zu groBem At) ergibt sich 


aus der Gestalt zur Zeit t durch Verschiebungen um 


dg |? ‘ 
AG, = w, 5 At (in der ¢-Ebene). 
n 
{7 
16 = ae 
ae 
a 
1 } he 
a Uso bt =F5 
Za H \200..-H 
Zz 010 | 0294 
42 4 Z 
age 9133 | 9324 
A 9763 9365 
10 Le 025 | 0404 
f 7 033 | 0436 
/ 9 | 2 gst | 547 
Vo 031 | g7or 
98 ; 5 ° 151 | 0888 
/ kK an | 1033 
267 
06 2 
yj f 
y] 5 ° 
OY z) 
W \ti\ \ 
| 
4 
iG 
Srey ae SEAN bet 
0 STE ae aay a | 7 ja a ee: 


Abb, 11. Unstetigkeitsflachen in der ¢-Ebene. Gestrichelt: Bild der Helmholtzschen Unstetigkeitsflache. 
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Abb. 12. Belegungsdichten der Zirkulation in der ¢-Ebene. 


Ebenso ist fiir einen Abschnitt der Unstetigkeitsflache zur Zeit t + At, der aus einem Abschnitt 
zur Zeit t durch senkrechte Verschiebung hervorgegangen ist, 


d 2 
T(t + At) = I'p(t) + yz, Z| At, 


woraus man y;(é + At) bestimmen kann. Dann berechnet man wieder die Geschwindigkeiten 
w(t + At) (mit Hilfe graphischer Integration) usw. Auf diese Weise wurde im Verlauf von zehn 
Verschiebungen die Entwicklung der Unstetigkeitsflache bis zur Zeit Vo. t/H = 8/3 verfolgt. Die 
graphischen Darstellungen Abb. 11 und 12 zeigen die verschiedenen Stadien der Wirbelflache und 


Ingenieur-Archiv 
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der Belegungsdichte in der ¢-Ebene sowie Abb. 13 acht dazwischen interpolierte Zustande mit 
gleichmafigen Zeitabstanden v,, At/H = 1/3 in der z-Ebene. Die Gesamtzirkulation der Unstetig- 
keitsflache ist in Abb. 14 in Abhangigkeit von der Zeit dargestellt. Zum Vergleich ist in Abb. 14 
die Zirkulation bei ahnlich veranderlicher Unstetigkeitsflache mit eingezeichnet. 


Es zeigt sich, daB der Anfang der Unstetigkeitsflache mit fortschreitender Zeit mehr und mehr 
in den Anfang der Helmholtzschen Unstetigkeitsflache asymptotisch iibergeht (Abb. 11). Vermutlich 
stellte letztere auch den Endzustand dar, wenn die Strémung nicht vorher instabil wurde. 


0 04 08 12 16 20 24 


Abb. 14, Gesamtzirkulation der Unstetigkeitsflache bei unendlich breiter und endlich breiter Platte als Funktion der Zeit 


: =F Yeo t\1/3 
estrichelt: 2s es) 
(« che lave H 0,635 ( 7] ) . ) 


4, Experimentelle Priifung und Erliuterung der Versuchseinrichtung. Die im Vorhergehenden 
behandelte Strémungsform, nach welcher die Platte von der Zeit t = 0 ab mit konstanter Conant 
digkeit v., angestrémt wird, laBt sich schwer im Strémungskanal ausfiihren. Die Platte miiBte 
Sa schon sehr schnell in den Kanal geschoben werden, wobei wohl Stérungen mit in die Strémung 

amen. 


Man erhalt die gleiche Strémungsform (vom Standpunkt ei i i 
punkt eines mit der Platte mitb t 
Beobachters), wenn man umgekehrt die Platte durch die ruhende Fliissigkeit zieht. In aves: Yall 
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hat man die zunachst ruhende Platte zur Zeit t = 0 augenblicklich auf eine alsdann konstant 
bleibende Geschwindigkeit zu bringen. Das ist genau genommen natiirlich unméglich. Trotz der 
damit verbundenen Bedenken wurde dennoch dieser Weg gegangen, wobei die Abmessungen der 


Apparatur so gehalten waren, daB die Beschleunigungsstrecke kurz gegen die erste Beobachtungs- 
strecke blieb. 


Abb. 15 zeigt eine Skizze der Versuchseinrichtung. Zwischen zwei Gleitschienen G ist ein Schlit- 
ten Saufgehangt. Mit Hilfe eines Bandes B, das ein Motor M langsam auf eine Rolle R aufwiecklt, 
wird der Schlitten in Pfeilrichtung gezogen. An der Riickseite des Schlittens ist eine 6 cm breite 
Platte P montiert, die bis dicht iiber den Boden in ein wassergefiilltes Bassin eintaucht (Eintauch- 
tiefe 36cm). Dadurch ist eine nahezu zweidimensionale Strémung um die Platte gewahrleistet. 
Unter dem Schlitten ist ein Photoapparat Ph aufgehangt, dessen Objektiv sich dicht iiber der 


—— oe Wasserobertliche 


Seitenansicht Riickansicht 


| Ausiritt der Farbtlissigkeit 


| 


Abb. 15. Versuchseinrichtung, 


Wasseroberflache befindet. Bei Annaherung des Schlittens an den Dorn D st6Bt letzterer durch. 
ein Loch in der Platte auf den Ausléser A des Photoapparates. Zu halber Eintauchtiefe tritt nach 
Offnen eines Regulierhahnes an beiden Plattenkanten eine Farbflissigkeit aus, die durch ein Réhr- 
chen an der Plattenvorderseite herab bis dorthin geleitet wird. 


Die an den Kanten vorbeistreichenden Flissigkeitsteilchen werden also eben dort angefarbt, 
so daB die zu irgendeinem Zeitpunkt vorgefundene Farbverteilung Streichlinien darstellt. Die- 
jenigen Streichlinien aber, die an den Plattenkanten vorbeifiihren, markieren die Unstetigkeits- 
flachen selbst (genauer einen durch sie gefiihrten Schnitt mit der Strémungsebene). 


Die Gestalt der Wirbelflichen hangt (bei gegebener Plattenbreite 2 H) nur vonSder GréBe 
Vx, t = 1 ab, wenn I die Strecke ist, die die Platte seit Bewegungsbeginn zuriickgelegt hat. Die ver- 
schiedenen Strecken, die die Platte bis zur Auslésung der photographischen Aufnahme durchlaufen 
soll, sind durch Marken auf den Gleitschienen abgetragen. 


Mit der geschilderten Versuchsanordnung ist es freilich nicht méglich, dieselbe Wirbelflache in 
verschiedenen Stadien ihrer Entwicklung zu photographieren. Vielmehr muf fiir jedes weitere 
Stadium der ganze Vorgang von neuem eingeleitet werden. Auf diese Weise entstanden die bei- 
gefiigten zehn Bilder (Abb. 16). Die auf der oberen rechten Bildecke notierte Zahl bedeutet die 
zum Zeitpunkt der Aufnahme von der Platte zuriickgelegte Strecke / = v,, tin Zentimetern. | lauft 
von | bis 10. 


Die zugehérigen ,,dimensionslosen™ Zeiten v., t/H sind (mit H = 3 cm) der Reihe nach 1/3, 
2/3,..., 10/3. Fiir die gleichen Zeitpunkte, von 1/3 bis 8/3, wurden acht errechnete bzw. inter- 
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Abb. 16. Photographien der mittels Farbfaden sichtbar gemachten 
Unstetigkeitsflachen, 


polierte Zustande der theoretischen Un- 
stetigkeitsflache im gleichen Mafstab dar- 
gestellt, was einen unmittelbaren Vergleich 
der Experimente mit dem theoretischen 
Resultat gestattet (Abb. 13). Fiir v. t/H 
>6/3 wird die Ubereinstimmung merk- 
lich schlechter. Dies ist evtl. darauf 
zurickzufiihren, daB die wirkliche Wirbel- 
flache allmahlich instabil wird, worauf 
auch die photographischen Aufnahmen 
(Abb. 16) hindeuten. 


5. Zusammenfassung und Schlubbemer- 
kung. Die vorliegende Arbeit behandelt 
die Ausbildung eines symmetrischen Paares 
von Wirbelflachen an den Kanten einer 
senkrecht angestrémten Platte wahrend 
einer Zeitdauer, in der diese Strémungs- 
form mdglicherweise noch als stabil an- 
gesehen werden darf. Als Ausgangspunkt 
hierzu dient die bereits von Anton durch- 
gefiihrte hier aber nochmal auf andere 
Weise wiederholte Berechnung einer Wir- 
belflache an der Kante einer einseitig 
unendlich breiten Platte. Als wesentliches 
Ergebnis ist hervorzuheben, daB die Wirbel 
sich von der endlich breiten Platte schneller 
entfernen und ihre Zirkulation schneller 
zunimmt als beider unendlich breiten Platte 
und daf ferner die von der Platte abge- 
hende Wirbelschicht sich in ihrer Form 
mehr und mehr der Helmholtzschen Unste- 
tigkeitsflache nahert. 

Die anschlieBend durchgefiihrte experi- 
mentelle Priifung zeigt im stabilen Bergich 
gute Ubereinstimmung mit den theoreti- 
schen Resultaten. Die Frage nach dem 
Instabilwerden der Wirbelflache wurde 
hier nicht naher untersucht und bleibt 
einer weiteren Behandlung des Problems 
uberlassen. Bekanntlich wird die Strémung 
mit der Zeit unsymmetrisch; die Wirbel 
lésen sich dann abwechselnd von den 
Kanten ab, wo wieder neue entstehen und 
formieren sich zu einer WirbelstraBe, 
wie sie von v. Kdrman beschrieben wurde. 

Zum Schlu8B méchte ich Herrn Professor 
Betz fiir die Anregung zu dieser Arbeit und 
fiir wertvolle Hinweise danken. 


(Eingegangen am 25. Mai 1960.) 


Anschrift des Verfassers: Dipl.-Phys. E. Wedemeyer, Gottingen, BunsenstraBe 10. 
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Zur Berechnung der Unterwassertragfliigel hei wellenférmiger Anstrémung 
Von W. H. Isay : 


1. Einleitung. In zwei vorangehenden Arbeiten! habe ich Unterwassertragfliigel als ebenes 
Strémungsproblem (unendliche Spannweite) behandelt. Und zwar enthalt Arbeit I eine weitgehend 
vollstandige Untersuchung des stationdren Falles bei homogener Anstrémung insbesondere auch 
fiir zwei hintereinander fahrende Tragflichen (Tragfliigelbboot). Die Arbeit II enthalt die 
Theorie eines Unterwassertragfliigels bei instationarer, wellenférmiger Anstrémung. In der vor- 
liegenden Arbeit wird zunachst die Theorie der Arbeit II auch auf den Fall zweier hintereinander 
fahrenden Tragflachen ausgedehnt, und anschlieBend werden Zahlenbeispiele durchgerechnet und 
diskutiert. Dabei wird die Kenntnis der Arbeiten I und II vorausgesetzt, so daf} wir uns hier teil- 
weise etwas kiirzer fassen kénnen. 


2. Zwei hintereinander fahrende Unterwassertragfliigel. Bei der Behandlung des gleichen 
Problemes in stationarer Anstrémung (vgl. Arbeit I, Ziff. 4 und 10) hat sich gezeigt, daB die 
Wirbel der hinteren Tragfliche 2 das Strémungsfeld an der vorderen Tragflache 1 kaum noch 
beeinflussen (Abb. 1). Das liegt daran, daB weit vor den Tragflachen das von ihren Wirbeln indu- 
zierte Geschwindigkeitsfeld verschwindet; dieses gilt nach Formel (II, 15), (II, 16), (II, 22), (II, 23) 


auch bei instationdrer Anstr6mung. Wir werden deshalb das von den Wirbeln des hinteren Trag- 
fliigels induzierte Geschwindigkeitsfeld in der Randbedingung am vorderen Fliigel vernachlassigen. 
Fiir die Randbedingung am Fliigel J ist das Geschwindigkeitsfeld des Fliigels 1 genau in der Form 
(II, 15), (II, 16), (IT, 22), (II, 23) einzusetzen. GemaB der in Arbeit Il entwickelten Theorie wird 
angenommen, daf der Tragfliigel sich in der mittleren Tiefe h, unter dem mittleren Wasserober- 
flachenniveau befindet. Dabei werden die instationaren Tiefenschwankungen vernachlissigt, die 
durch die wellenférmige Wasseroberflache und die Eigenbewegung des Fliigels bedingt sind. 


Fiir die Randbedingung am Fliigel 2 mu dagegen das vom Fliigel J induzierte Geschwindig- 
keitsfeld beriicksichtigt werden. Aus der stationaéren Theorie der Arbeit I ist bekannt, daB durch 
den vorausfahrenden Fliigel 1 die Wasseroberflache in der Umgebung des Fliigels 2 deformiert 
wird. Genau nach der Methode von Arbeit I, Ziff. 6b werden wir auch jetzt den stationaren Anteil 
dieser Deformation bei der Randbedingung am Fliigel 2 beriicksichtigen; d.h. der mittlere Abstand h, 
vom mittleren Oberflachenniveau muS entsprechend angepafit werden. Vernachlassigt werden 
dagegen die Tiefenschwankungen, die durch die wellenformige Wasseroberflache sowie durch Tauch- 
und Stampfhewegungen des Fliigels bedingt sind.? . 

Mit dem angepaBten h,-Wert ist das Geschwindigkeitsfeld des Fligels 2 in die Randbedingung 
am Fliigel 2 genau in der Form (II, 15), (II, 16), (II, 22), (II, 23) einzusetzen. 

Wir haben ferner das vom Fliigel J] am Ort des Fliigels 2 induzierte Geschwindigkeitsfeld zu 
berechnen. Dabei kann zunachst das Geschwindigkeitsfeld der gebundenen Wirbel des Fliigels 1 
am Ort des Fliigels 2 (— a, + L < %< a, + L) ohne weiteres aus (II, 15), (II, 16) entnommen 


1 ; Ing.-Arch. 27 (1959) S. 295; Arbeit I; Ing.-Arch. 29 (1960) S. 160; Arbeit II. 
e rpeerais Side ae PAS BE Theorie noch wesentlich komplizierter. An Beispielen werden 
wir spater sehen, dafB fiir kleine Wellenamplituden die obige Vernachlassigung gerechtfertigt ist. 
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werden. In Anbetracht der groBen in Frage kommenden %,-Werte ist es fiir die numerische Aus- 
wertung zweckmaBig, an Stelle von (II, 10), (II, 12), (II, 13) die folgenden Formeln zu benutzen: 


.l—e 


g u : Be ed 
ices WHET Ma) a lee o RCL 


T, (xg, — hy, §4, /4) = 2 
a (x%_ — 2o| My 2 (%_ — 4) 
L 


— pak ty = %,— E, . + 0 
per or eee) e up ate (ms (= —£)] 4-019, (1 - (Xp —§)) h iz en) 
0 


& («+i f) (hyp h,—i x.+7%§) : & f : 
I,(% — by &) =— 08 E,((0 + 1p) (hip tim —i6)] 
! 
ee [44 — fer}. (2) 
(a + iB) “(i — iS, — hy — hy) (4 + iB) = (i x, — i &;y — hy — hy) 


2 
5 0 


I 


Entsprechend ist die Formel fiir I,(x,, —hg, €,). Zum Beweis bemerken wir, daB sich (1) aus (II, 10) 
durch asymptotische Auswertung des dortigen Integrals ergibt, wahrend (2) aus der bekannten’ 
asymptotischen Darstellung der E;-Funktion gewonnen wird. 

Wesentlich schwieriger ist die richtige Erfassung des Geschwindigkeitsfeldes, das die freien 
Wirbel des Fliigels 1 am Ort des Fliigels 2 induzieren. Die Formeln (II, 22) und (II, 23) sind unter 
der Voraussetzung abgeleitet, daB die freien Wirbel in unveranderter Starke hinter dem Fliigel ab- 
schwimmen; dieses ist jain Wirklichkeit nicht der Fall. Da nun jeweils stromaufwarts eines Wirbels 
dessen Einflu8 mit zunehmender Entfernung schnell abklingt, und die ersten (und somit entschei- 
denden) freien Wirbel hinter dem Fliigel noch annahernd ihre volle potentialtheoretische Starke 
haben, diirfen Formel (II, 22), (II, 23) sicher fir x-Werte in der unmittelbaren Umgebung des 
Eigenfliigels benutzt werden. Keineswegs kann man jedoch erwarten, dafs diese Formeln weit 
hinter Fliigel J, d.h. am Ort von Fliigel 2, noch richtige Ergebnisse liefern. Es liegt somit nahe, bei 

g 

der Zirkulation der freien Wirbeleinen multiplikativen ,,Abklingfaktor“ e : up aes (6X =o) 
anzubringen, der naherungsweise den Zerfall der freien Wirbel beschreibt. Der Nachteil dieser 
Methode ist jedoch, daB der Wert von € aus der Theorie nicht entnommen werden kann, so daB 
hier eine Unsicherheit besteht.2 Gliicklicherweise hat sich an den Zahlenbeispielen gezeigt, dab 
der Wert von ¢ in der Randbedingungsgleichung keinen starken EinfluB ausiibt. Dieses liegt vor 
allem daran, daf die (in Ziff. 3 noch zu berechnenden) von der Stampf- und Tauchbewegung des 
Fligels stammenden Geschwindigkeitsanteile einen wesentlichen EinfluB auf die Randbedingung 
haben, und der Einflu8 der freien Wirbel demgegeniber zuriicktritt. 

Es geniigt somit, wenn man einen méglichst guten und gegeniiber etwaigen Variationen un- 
empfindlichen Mittelwert von ¢ in die Rechnung einsetzt. Um diesen Mittelwert zu gewinnen, 
gehen wir folgendermaBen vor: Wir wollen annehmen, daB der Ort, an dem die freien Wirbel des 
vorderen Fliigels nur noch etwa 5% ihrer urspriinglichen potentialtheoretischen Starke haben, 
zwischen den Stellen x = L/2 und x = L liege. (L ist der Abstand der beiden Fligel). Wegen 
e 3 ~ 0,05 gilt dann 


2 
Gert Spe" : 


* Vel. Lésch-Schoblik, Die Fakultat und verwandte Funktionen, S. 120, Leipzig 1951. 
* Man konnte daran denken, hier irgendwie von dem bekannten Zerfallsgesetz eines ebenen Einzelwirbels 
in laminarer Stromung 


i Soop OM CE ra Aaa (v = kinemat. Zahigkeit) 


Gebrauch zu machen. Jedoch ist die hier vorliegende Strémung voll turbulent, und dann erfolgt der Zerfall 
der freien Wirbel ganz wesentlich schneller und mit starken Zufallseinwirkungen, so da8 aus diesem Gesetz 
keine quantitativen Hinweise fiir die GroBe von ¢ erhalten werden kénnen. AuSerdem haben wir es ja hier 
nicht mit einem Kinzelwirbel zu tun. — In der Theorie des Voith-Schneider-Propellers (wo ebenfalls die freien 
Wirbel beim Durchflu8 durch die hintere Fliigelreihe schon weitgehend zerfallen sind) ist es dem Verfasser 
gelungen, eine summarische Naherungsdarstellung fiir das Geschwindigkeitsfeld der in der turbulenten Stré- 
mung zerfallenden freien Wirbel abzuleiten; vgl. Ing.-Arch. 24 (1956) S. 148. Diese enthalt keinen freien Para- 
meter mehr, und sie hat sich bei der Durchrechnung von Beispielen bewahrt. Fir ihre Ableitung war jedoch 
die Rotationsbewegung der Propellerfliigel wesentlich, so daB eine ahnliche Darstellung fiir den Fall des Unter- 
wassertragfliigels nicht gefunden werden konnte. 
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Berechnet man nun mit einem noch unbestimmten €-Wert das Geschwindigkeitsfeld ur (2, C) 
v¢,(2,¢) der freien Wirbel des Fliigels 1 am Ort des Flii 
verwendete Mittelwert gegeben durch 

2 


? 
gels 2, so ist der von uns fiir die Rechnung 


u u2 
e l res . zi 
5 a a a a‘ 
Ly(2)m == a | w(2,6) aC, [oy (2)Im =z 5 | 17 (2,6) de. (3) 
ue U2 
3 7 3 +E 


Fur die Berechnung der Geschwindigkeitskomponenten (3) ist es zweckmaBig und auch aus- 
reichend, den Fliigel J durch einen in seinem 1/4-Punkt (d. h. etwa Druckmittelpunkt) é, = — a,/2 
angebrachten Punktwirbel der gleichen Gesamtzirkulation 


Ty) = Al) + aM eror 4 AY ie 


zu ersetzen. Dann haben die freien Wirbel hinter dem Fligel also die Starke 


a, . @ 
L 


| 5 3) & ( ) (x 3) 3 3) 

d i exer t xen ee X+ ¢ x+ 

—i= yer & eo) efeket 7} — AM e—iote age. ; ere 
0 


Um nun die Geschwindigkeiten u,(2,¢), v¢(2,¢) aus der grundlegenden Formel (II, 9) zu be- 
rechnen (analog Arbeit II, Ziff. 4), benétigen wir einige komplizierte Integrale, die wir nachfolgend 
zusammenstellen. In Anbetracht der groBen in Frage kommenden x,-Werte ist es dabei méglich, 
sich auf den .,quasiasymptotischen*‘ Wert dieser Integrale zu beschranken. Dadurch ergeben sich 
betrachtliche Vereinfachungen, die bei unserer Naherungsdarstellung natirlich angenehm sind. 
,,Quasiasymptotisch* soll bedeuten, daB wir beschrankte Integralausdriicke, die mit dem Faktor 

§ ai 
e i 49 (« : : ae sind, vernachlassigen, und da wir ferner in der asymptotischen Entwicklung 
der E,-Funktion die beiden ersten Glieder beriicksichtigen, z. B. 


: (C410) Pa ‘ a, \\ 1 1 
> al ‘)p,(c +i Qs (x, 2) F SA GleT z anit: 
Q C+ 12 %(at a e+ 5 («+ 9) 
0 0 


An Zahlenbeispielen hat sich gezeigt, daB diese quasiasymptotische Darstellung vollsténdig aus- 
reicht. Es gelten dann mit den Abkiirzungen 


Ay, = abe - Z| ? Ay, = a (hy + he) (4) 


die folgenden Integralformeln: 


= lta, (Barat gh 2 Gee Hugs ecu Bia” (5) 
ree (x, — XP (hy + Ay? Q=it AM a Q—it AR, 
—a,/2 
Cc 
—£ (x14) Cee (ee) aX Sh Las 1 1 1 (6) 
or (@—XP+ (thy Qi Ay (@—if)? Ad’ 
—a,/2 
= oo ; 
— 2 (x+ cri — i M3 (x,—X) hee Pir 
: “al ") wie [c. (BE (—%9) + #8,(S3 G —%)) +25 
re iy Hey My 8 
2 Cre gee) | 8 ax 70 fly @ tt An 
a 4g am Oa iS FE aS 
p2 aL in (= —¥%)| 
0 y 
1 eft Fi, i eft A, L i eft A, (7) 
~ Qi Ag, — t Hey Tig (Qt) (Ag, —t Ha)? (Q —il)® (Ag, —i H,)* 
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foe) co 

& 4 i 1 tS Xo—. . 

pevuter anc = as eae ip ena : : (8) | 

— > | 

e—p+2nQ+ @ Be “"R QD (Q=if) (4g +i Aa? | 

—a,/2 0 | 

co 8000 1 

Pee 3 ay i i ee a Hy 

ae (x+ +) C+ BY ey Bye Ut ua’ X) rar 1 1 ij 

eet te Q+ & Mag © ~ OF Q=iD) An Ft A | 

—1a7/2520 | 
gto ees 1 i 1 

: (9) 


TO) Gattig  POAOare ig 


Die Integralformeln (5) bis (9) entsprechen der Reihe nach den Integralformeln (II, 17) bis (II, 21). 
Es sei noch bemerkt, daB Formel (6) auch gilt, wenn h, + h, = 0 ist. Den Beweis der obigen In- | 
tegralformeln bringen wir in Ziff. 6. 


Mit Hilfe von (5) bis (9) kénnen die gesuchten Geschwindigkeitskomponenten (3) leicht an- |} 
gegeben werden. Nach Ausfiihrung der Mittelungsintegration tiber € mit den Formeln 


| 


6 u2 
gL (FE) + 2 | 
1 gl dpe Vg hl (0 tl ee au 1 iy gl i ae 
[Oil oa e(Ft 9) arctg (27 a)+3 “Gays @ a=, | 
3u2 iL i 
Se S | 
g 
6 ue 
gL 6 us 
suze0 — ey 
Dig bier Us 9 8 aay rt age Pe Sata pi) + 4 
Bug f m—Q+il 38 UG aver) arcte (TT a=) aaa a 
are | tae 
3 uz gl 
=“ 
ee 1 
one 
6ue 
GL 
1 gL dt Pig Lk 1 1 1 
=a | joe = Se | es — Ss ) = | I, = Mittelwert 
3u2 gl gl 
“ 


ao) | 


Und einigen elementaren Zwischenrechnungen erhalt man schlieBlich unter Beriicksichtigung von } 
My be = 2? 
[up (2)]n = Uf (%_, — hg, t) = x = AD {---tebot Saal ae A AE. -j}e tot, 


27 Up : 
2% My Sue 1 270 : 1 ; 

{. 7 -} a ea Ay, Cm es Ay, - | Me e— He Aa, es oe a eae 
gies OF TET at : w—O+it),¢ ON} 
ees #) |p =| : | 

21 21 pt | 

1 @ 1 : 1 —1\\———— 

2 = = = — —— (1) eee vu ae WE (1) —iwt. W| 
[vp,(2)]m M4, (e» — hay t) = o> Ar) fs Jokes +55 SADE Fert; 

270 le oa 1 27 | 

{- y + See M1 Has @ im, Ary —_—_—_————— — Me el Ay, tf, Ay, | 1 H 

Ma D4 ily aly : wari, ¢ 9) i 
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3. Die Stampf- und Tauchbewegung der Fliigel. Bevor wir dazu tibergehen kénnen, aus den 
Randbedingungen an den Fligeln die Integralgleichungen zur Bestimmung der Fligelzirkulationen 


Vir» t) = OE) + AM(E,) ef # + SME) eto, 


Ty) = AY + AP ot 4 AM e-tot = f y(é,0) dey, 


a ete (13) 
Yo(Eo, t) = 65) a 6P(E,) ebot 6(é,) em ior, 
== a,+L 
L(t) = Ay? + AP? eto + AP e-iot — J Yo(&, t) dé, 
—a,+L 


abzuleiten, miissen wir noch die Stampf- und Tauchgeschwindigkeiten der Fliigel berechnen. 
Diese diirfen in der Randbedingung am Fiiigel keinesfalls vernachlassigt werden’; im Gegensatz 
zu der bei einem Fliigel allein auftretenden Tauchbewegung (vgl. Arbeit II, Ziff. 5) ist bei einem 
Zwei-Fliigel-System auch die Stampfhewegung zu beriicksichtigen. Es sei M die Masse des Trag- 
fliigelbootes und s sein Schwerpunkt, ferner L — L, + L, (Abb. 2). Wir nehmen an, daB die 


Masse der Fliigel gegeniiber der Bootsmasse vernachlassigbar ist und sehen auch von Dampfungs- 
-einfliissen des Wassers ab. Dann lautet die Bewegungsgleichung des Schwerpunktes (und damit 
auch des ganzen Bootes, wenn Durchbiegungen vernachlassigt werden) fiir das Tauchen 


Mi © —Mg + ou lilt) +e w L(t) =—Mg 
+ Q@ Uy (AS? + Ay eos 1 Ay? e—iot) | ou, (4? 4 AP eiot 4 7) e—iat), 
Da die Masse M im zeitlichen Mittel mit dem Auftrieb im Gleichgewicht ist, gilt 
Mg =o wu (Ay Ae). 
Dann liefert die Integration fiir den Ort des Schwerpunktes unter Beriicksichtigung der Periodizi- 
titsbedi ee) 
atsbedingung a ee 
g APLA® . g AP+ A® Biy.- 
n(*) = oa AD AS ° ‘ADE A® ° oF (14) 
Wie schon in Ziff. 2 gesagt, soll angenommen werden, daB fiir die Theorie 7 ~  gesetzt werden 
kann, d.h., daB auch die Tauchtiefen der Fligel ungefahr konstant bleiben. Dagegen ergibt sich 
die fiir die Randbedingung wesentliche Vertikalgeschwindigkeit 
. gi AP+A® . gi AY + A® ar: 
i) = EAE AE os 5 IE AE set as 
Bezogen auf den Schwerpunkt als Drehpunkt lautet die Bewegungsgleichung fir das Stampfen 
bei kleinem Stampfwinkelgy und mit @ als Trigheitsmoment des Bootes 


Op 1,608 y°Q ty y(t) + Lz 008 pg ty Ta(t) & —@ to (Ab? Ly — Ad” Le) 
a Eso (4. L,— A?” L,) 2?! —@ Uy (A ay L,) earor, 


1 Es sei denn, man beschrankt sich auf die Untersuchung von Fliigeln, die etwa in einem Priifstand fest 
‘montiert sind. Dieser Fall ist in der hier entwickelten Theorie mit enthalten; man braucht nur alle Stampf- 


und Tauchgeschwindigkeiten wegzulassen. 
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: 23 : il 
Auch fiir die Stampfhewegung besteht die Periodizitaétsbedingung ( +) = 900: Weiter i 
22/0 


sel | y(t) dt = qm) =; denn ein von Null verschiedener mittlerer Neigungswinkel gy) kann | 
7 


: . « . | 
im Rahmen der hier vorgelegten Theorie ohne weiteres mit in die Fligelprofilneigungen einbezogen | 
werden. Dann ergibt sich 
1 (Qharaee mage Oe \ eee. al! OE eee it 

AW L, = Ab? Ly; ates ap ae ; fa AD AW} (16) | 


und die Integration liefert mit (16): 


ig ML AP AP —-AP AM. , tg ML A® A@ — AP AY ee 17) | 
Vt) a ore AO AD Binns o O-- [APFAPP$ «aig i 

1 ua ML? [A A —A® AY. AD A® — A® AW ae 
Y(t) = 22 aE ra) ‘ [ap BE AD PO le a [A® + AQP =O ietalt (18) 


Fiir die Randbedingung an den Fliigeln ist die y-Komponente der Stampfgeschwindigkeit 
wesentlich; sie ist gegeben durch —L,pcosy ~ —L,@ bei Fligell, L,pcosy » Lop bei 
Fliigel 2. Dagegen werden wir die entsprechenden x-Komponenten L,gysing ~ L,yqy bzw. 
— L,g~sin y = — L, y¢ als von héherer Ordnung klein vernachlassigen. Dieses ist bei der hier 
entwickelten linearen Theorie selbstverstandlich. Analog haben wir in den Bewegungsgleichungen 
fiir das Stampfen und Tauchen die Fliigelkrafte in y-Richtung in der Form K, = Q uy I’ einge- 
setzt und nicht in der genaueren Form (Ziff. 5) unter Beriicksichtigung der Stérgeschwindigkeiten. 
Noch unwesentlicher ist die Vernachlassigung der Kraftanteile K, sing, die an sich in die Be- 
wegungsgleichung des Stampfens mit eingehen. 

Andernfalls wiirde die Theorie durch nicht lineare Glieder sehr kompliziert, und es miBte 
dann auch die Annahme der konstanten Fligeltauchtiefe aufgegeben werden. Damit wiirde unsere 
grundlegende Ausgangsformel (II, 9) fiir das Geschwindigkeitspotential eines Wirbels ihre Giltig- 
keit verlieren. 


4, Aufstellung und Lésung der Integralgleichungen. Mit f;(x,) und f,(x,) als Neigungen der Fligel- 
skelettlinien lauten die Randbedingungen am Fliigel J und 2 


/ 1) — be — he ) + Ov, hy ) ef,  — by) — 10 + Ly OO) 
fil) is A ) Uy (%, cea hy, t) + Uy (x, ae hy, t) + es (x, a hy, t) + Ug 
Fr(%2) + p(t) = (19) 
0 (%2, — hy, t) + Vy, (%2, — hy, t) + Uf, (%2, — hg, t) + vy, (%o, — hg, ty + YF, (xg, — hg, t) — H(t) — L, p(t) 

Uy (Xo, ayy he, t) ar uy, (Xp, res hg, t) ae uf, (xo, ie hg, t) + uy, (Xo, os hg, t) + uf, (Cor An 3 hg, t) ta Ug ; 


Da nach den Voraussetzungen unserer Theorie sowohl y als auch u,, u,, uf klein sind, kénnen wir 
entsprechend unserer Bemerkung am Schlu8 von Ziff. 3 die Produkte @ u,, y u,, y uy in den Rand- 
bedingungen vernachlassigen und nur @ uy beriicksichtigen. Wie man leicht erkennt, wiirden 
diese Produkte in den Integralgleichungen AnlaB zu nichtlinearen Gliedern geben und damit die — 
Theorie sehr komplizieren. . 

Aus (19) entstehen durch Koeffizientenvergleich in 1, e'®', e—‘”! sechs Integralgleichungen, 
eS papi nur vier wesentlich sind, denn mit 6, ist ja auch 6, bekannt. Der Koeffizientenvergleich 
in 1 letert 


? 


2 mafiles) =z | OME) |g + Kul 2d} dB 
om a,+L 
—Puofile) =z f OG) Knlontd as +> | AMG) [Ae + Kelp Go| tas f O) 
— a, er z 
(-asi sa, +L Se sa+L). 


Dabei sind die Kernanteile K,,, Ky, K,, genau mit denjenigen des Systems (I, 24) identisch: — 
lediglich den Kern K,, haben wir diesesmal weggelassen, da wir aus Arbeit I wissen, dafi der 


hintere Fliigel praktisch keinen Einflu® auf den vorderen austibt. Die Auflésung des System (20) | 


-die stationdren Mittelwerte der Zirkulation u 


\ 
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erfolgt also genau nach den Methoden der stationdren Theorie aus Arbeit I, und alle Formeln 
kénnen wortlich von dort iibernommen werden. Damit sind dann 6N(&), 6(&), AS, AY) als 
bekannt anzusehen. Insbesondere sind also jedenfalls im Rahmen der hier vorgelegten Theorie 
nabhangig von der wellenférmigen Anstrémung und 
der durch sie ausgelisten Tauch- und Stampfbhewegung der Fliigel1 Fiir die effektive Durch- 
fiihrung der Rechnung muB auch der stationdre Anteil des Strémungsfeldes vollstandig bekannt 
sein, ehe mit der Behandlung des instationaren Anteils begonnen werden kann. 


Durch Koeffizientenvergleich in e! ergeben sich aus (19) die beiden folgenden Integral- 
gleichungen: 


g ee Se My a ily af 2g ars —u, —in So 
2 ta ig Bi (1 + tfi(m)) ¢ ee . 1 olan By (1 + i f;(x,)) e See : 
2 ML? 1 uz AY AP — A® Aw gi APL AP. giML A® AD AP — A® AD 
ee Oe sO gL [A@ + A@ @ AD+ A@ | iO 8 [AM + A@ 5 
een 1 CPST: eee 
=1 | HG) ie g + Ful, &) —22 n |S | dé, , (21) 
gh, EX, gh, - & %2 
g Bh My 2 My 7 g oe Ue 2 Ups 2 
2 Hi ag Bs (1 + i fa(%s)) € oe : Pie Us Te Ba(l +t falta) © rik ; 
ML? 1 wy APAP—APAP |, gi AP4 AM 5 gi ML Ae AD A® — A® AP 
— 2 Uy O Qe [AD + Ag lear) AP+42°—“ @ ©O 0 [AP + A@, 
a,.+L 
ee HG a a 1 2 1 iD, |%2— § 
an (E+. —Z)= 2 { Ped et heelee bs) — in| — las, 2) 
—a,+L 


Die Kernanteile k,,(x,,¢,) und kyo(x»,&) sind absolut von der gleichen Form wie (II,29). An sich hatte 


a 


‘ ae 
~ in Gleichung (22) an Stelle von = Ay hay (L to 7) stehen miissen — | 69(E) Bop (605 €1) dE, 


2 
@ 2 
In Anbetracht der in k,, enthaltenen Naherungsdarstellung (11), (12) fiir die von den freien Wirbeln 
des Fliigels 1 am Ort des Fliigels 2 induzierte Geschwindigkeit haben wir uns bei dem Ausdruck kp, 
mit der +/,—%/,-Punkt-Methode begniigt. Es ist 

C7 


a a g (4o,—f3 Ho, ake ea Smee AVE Ps) 
Bi(h+3-— 9) =a ane ta ths (1-44) 4+ 4.9 


ease 2 
uo 


SS Tae 
marae re 


1 [Re RS a al eee 1 | | 
@ Ms #421 @_* Ha 401 > 
acornme aera fearesrie 


(23) 


—¢+fs(1 +iG)L(L+9s—he—F) + i—fo)| 


b me oe OT Wee ere at pd 
Me 2 — bart ae eo e M1 Aa; @ Upc Se eG Hy 4te1 @ Hag 
4 (L 2° fo» 2 Hs) (a bs Mos 


. ON ELEVA eine 
- 2(1 a. Vesta Vent TIS e— Hite @ My Ag, 
( ifs) uo ae 


@ (2 porte TE ie it 9,2 1 1 
rei erree saat la=vel,,+ * iy Az, |((Q—ioP 


In (23) sind J, I, I, nach Formel (1), (2) zu berechnen. Entsprechend Gleichung (4) ist jetzt 


m 


Ay, =2(L +242). SchlieBlich bedeutet z. B. 
u 2 2 


ey EN ee en et ae oe cde at arr 


2 
0 


1 2 bg ed Bharata Es 
eer ae ut | my —Q+iC 


ae den Mittelwert gema Formel (10). 


1 Dieses ware nicht mehr der Fall, wenn in den Randbedingungen nichtlineare Glieder wie y wu, usw. be- 


‘riicksichtigt wiirden. 
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Die in Arbeit II, Ziff. 6 fiir die Integralgleichung (II, 27) mitgeteilte Auflésungstheorie kann 
fiir die beiden Gleichungen (21), (22) ohne weiteres iibernommen werden. Die durch die etwas 
anderen Glieder auf den linken Seiten von (21), (22) bedingten Anderungen sind leicht zu ther- 
sehen, so da® wir darauf verzichten, sie hier noch besonders anzugeben. Durch die AbfluBbedin- 
gungen an den Profilhinterkanten 

é(a,) =0, dP(a, + L) =0 
ergibt sich schlieBlich ein lineares Gleichungssystem fiir die Berechnung der Zirkulationskoeffi- 
zienten A$? und A, 

Damit ist nach den Gleichungen (14), (18) auch die Tauch- und Stampfhewegung des Trag- 
fliigelbootes bekannt. 

Genau wie in der stationéren Theorie Arbeit I, Ziff.7 kann fiir Naherungsrechnungen auch 
jetzt die 1/,—3/,-Punkt-Methode verwendet werden; d.h. die Fliigel werden durch Punktwirbel 
in ihrem 1/,-Punkt (~ Druckmittelpunkt) ersetzt, und die Randbedingung wird im °/,-Punkt 
erfiillt. Man vermeidet damit die tatsachlich sehr langwierige Berechnung der Fourierkoeffizien- 
ten b,, gemaB (II, 32) fiir k,, und hyo. 


5. Berechnung der Fliigelkriifte und Form der Wasseroberfliche. Die resultierenden Fliigel- 
krifte in «- und y-Richtung K, und K, werden wie iiblich nach dem Kutta-Joukowskischen Satz 
berechnet. Fiir Einzelheiten dieser Rechnung verweisen wir auf Arbeit I, Ziff.5 und Arbeit IT 
Ziff. 7. Es wird mit 9 als Wasserdichte fiir Fligel J 


Ky = 0 (4 ileres — h, j + uy, Fae —h, ) + uy, (—F— hy | L(t) » 
K! =—e| Hl) + E,90) +%(—$— hat) +8, (—$— hat) +e (—F— he HJ 


(24) 
und fiir Fliigel 2 


Fo 25 (¥ uy, (EF se ) 4 us,( ihe 7 + uy, (U3 she ) 
a, (tS ) + ug (E—$ he | T(t), 
Ko =p [ilo — 1, 9(t) + (t—%, ape ) a B(L—9, ai. ; + ap, Gas he ) 
+v,, ( —3, ps ) + v,, (Z —3, Pass 55 i Lay; 
(25) 


¢ Die Berechnung der Form der Wasseroberflache kann nach Gleichung (II, 1) erfolgen, d.h. 
ier 


1 /a® a®. ) o®@ a 
y ee 1 ta fr V2 Df, 
(*, 1) g ( Gere bo er te oe ) 5 
a (u(x, 0, t) ae Uy, (%> 0, t) 3 Uy, (x, 0, t) = Uy,(%, 0, t) =e Uf,(x, 0, t)) (26) 


Der zeitunabhangige Teil von (26) ist natiirlich genau mit der Oberflichenform der stationaren _ |} 
Theorie der Arbeit I, Ziff. 6b) identisch. Ein fiir alle x-Werte giiltiger formelmafiger Ausdruck 
1aBt sich fiir (26) wegen des Zerfalls der freien Wirbel schwer angeben. (Vgl. die Ausfiihrungen 
in Ziff. 2.) Jedoch laBt sich fiir die Umgebung der beiden Fliigel folgendes sagen: In der Um- 
gebung des Fliigels 1 gilt Formel (II, 39) unverandert; Fligel 2 iibt keinen Einflu8® aus. In der 
Umgebung des Fliigels 2 gibt Formel (II, 39) (natiirlich jetzt bezogen auf Fligel 2) den EinfluB 
or Fliigels 2 sen eres genoa , Anstrémung auf die Form der Wasseroberflache wieder. 
inzu tritt noch additiv ein Ausdruck AY, der d z ul i a i 
BOUCHE Mig ote cc en Emflu8 des Fliigels 1 auf die Oberflache in 


Wir benutzen die Formel (II, 15) fiir wu, (x, 0, t) und Formel (11) amt), == _ fiir uf (x, 0, t); 


ee ody, : | 
entsprechende Ausdriicke fiir ea ( a ; kénnen in analoger Weise aus (II, 9) unter Be- |] 
= = z ot | 
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nutzung der Integralformeln (5), (6), (7), (9) erhalten werden. Dann ergibt sich nach einigen 
elementaren Zwischenrechnungen 


1 A® gh 


A*Y (Xp, t) = 5— am E (x, 0,— 3, C “5 (x2, 0 23 1)—20 a “5 sin £ (x, + “1 | 
A® . 


1 a 
ie a Up Saute J+y u ett EY, 
mit 
1 a. 1 a Q 
{ eee iF) (m0, ae 5) t T ig) Ii (m 0, mq ce Iy(e,0 a 
Q —Q 
DB: I, (» 0 a= | (> 0, see 1 fe 0, Hs) (27) 
7 hom" e— 4 A, eit Any Fr f,— 2 e— Hs Aa; e— tH, Ay 
aA Mo— by 
Q(4— ; 1 
6 ae My Ay, g@— ty An anes et 
: My — be e- : fy — 2+ ale 
2 (Ug — 2) ; 1 
a a t Se SET RRs ita lorie ie pega g7r | eee en 
Pomme . - M,— 2+ a 


or 20 < 2 Bs i eet 2): 1 
(a Q ii) bear AE, \o— 7, 


Fir die Berechnung von I,, I, I, dienen wieder Formel (1) und (2). 


6. Beweis der Integralformeln. Wir haben nunmehr den Beweis fiir die Integralformeln (5) 
bis (9) nachzutragen. Hierzu verweisen wir noch ausdriicklich auf unsere Definition des quasi- 
symptotischen Wertes auf S. 203. 

a) Beweis von (5), (6). Wir fiihren den Beweis fiir Formel (6); ganz analog ist er fir (5). 
Mit den Abkiirzungen (4) und der Substitution X = x, + (h, + h,) @ geht das Integral (6) in die 


Form iiber 


oO a he 0 db 
— p— 401 (+i 2 — HH}, (€+1 2) — A,, (€+1 2 A, (+i Q) 8 ° 
e ( ) e ( ) rv) we ( ) é ( ) rv) Ly r) 


— A,,/Fp, 0 
denn das Integral iiber ? von 0 bis 00 wird in der quasiasymptotischen Darstellung vernachlassigt. 
Durch direkte Integration und Vergleich! mit der komplexen Reihenentwicklung fiir die E,- 
Funktion findet man 

Ay, /H; 


odd 
—An+i2) | eHaC+ia)e 
e€ Je e fPao1 ae 
0 
Ay, /Hoy (a iQ) (8 i) 
B54 Het ee ote AO aka, Te = 
a5 \(C-+4Q) ES -- a do 
0 
A,,/Hoy H (4 iQ) (9 + i) 
- A,,+iH,,) (C+i2 al = = ; = a 
a Measteg e+ 1 aes 


== = Cr arn) [E; ((¢ + i Q) (Ag, —t H,,)) =H, (Hoy (Q—i ¢))| 


fo ent thy) 6452) [TE ((C + 4 Q) (Ags + i Hay) —E; (Hes (— 2 + 46))| 

eae 1 1 ttre 1 i 1 

Se Daas e mee Fea ae on 2 (Q—itp (Cressy (4;,+ i H,)? 
- oe ees 
Oe Ae (Oily 4’ 


1 Vel. z. B. Lésch-Schoblik, FuBnote S. 2., S. 118. 


. 
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denn A3, > H3,. Im Spezialfall H,, =0 verlauft der Beweis ahnlich mit Benutzung der Re- 


lationen 


tT — Ap, 


co 
e— 4 (6-419) fe c+iayo ee 2 es — g—An E+i 2) [E; (Ay, (C +7 2)) + iz. 
ae 


b) Beweis von (7). Mit der Abkiirzung (4) und der Substitution t = as ater %_) geht das 
Integral (7) in die Form iiber 


fp — 426 —im tia) 
e— 4 (F+i 2) | Ca C,(t )—i Sr) +i > 


—p,As, — yA, 


S. Ra crim ee 
dt + | ae i e een 
6 


Wir betrachten zunachst den ersten Anteil; es ist 


1 
ee 1 : : eo —— (—ip,tiQ)r © (eo) 
me aS Sea nee as My fod] 
= fy Agy Tv 2 =F a ae Tt ? — My Ag, 
f ~( 2) (4 +iQ)A 
—— 6+i2)7 a > ye tp, tis 21 i 7 
Pre ee e Fn [Cults Ans) +i Sit Ana) $i FS 
—h1A 
7 (€—ip,+i 2) A), 
ot Seti [eee GRO en ; we 
ge SPITE if ie O— Ay, fy — Q+ iC liz 4 Ei Uy Ao3)] 


0 


a airae bin + EC An +12 ny)] 5 


dabei wurde von der Formel (II,10) und der Reihenentwicklung der Ei-Funktion Gebrauch 
gemacht. 
Untersuchen wir nun den zweiten Anteil des Integrals (7), so kann das Integral iiber t von 0 


bis co in der quasiasymptotischen Darstellung vernachlassigt werden, und es bleibt 


Hy Ay, 


" = +9)" sat 
Aatesioy | ef ie e Jone d} dr 
6 


191 as i tT—i 
i @4i0)"7 tT+id eft Sa lah ee 
eee * -- drt dd 
Tat t™—iO 
0 


me lz, (( ne) (40, we 2))_ E, i (Os 0)| do 


fs) 


pe a | es E,((¢ +4.9)(4,, +2) a0 
0 vf 


My, HL, 
= —j e— 4 (6+i2) 


— 7 e— 4a E412) 


Hy By 
a—i | | My 1 wi 1 | a = i py ( gts Hay 6) 
d 


 [Q=it OF in, 4, | Q=ibt (9 + i py Ag,)* (Q—it}? \Ha +1433 i Ay 


Hy Hy, 


ee te O—ipy Ag; ee titan ae 
zl +a=R aa) | pase +S a oUt Be 


Bea fa HE, i if Le mA : . : 
e(OLtOP ae iA +7)—p4( TH ae ei An TE (Uy Hy, +1 My Ao) —E; (i wy Ag,)) . 
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Setzat man nun die beiden Anteile des Integrals (7) wieder zusammen, so erhalt man in quasi- 
asymptotischer Form nach kurzer, elementarer Zwischenrechnung das behauptete Ergebnis. 


c) Beweis von (8), (9). Wir geben den Beweis fiir Formel (9 i it 

a a)? $ log ist er fiir (8). 
Zunachst erhalten wir unter Vertauschung der Int tionsrei 1 eee oe ss = 
Integration iiber X fiir Formel (9) Oem entolee mack “Ausfibtung dés 


ee) 


—if eo An eit As du 
| Hut BEDE DWH O= iO) 


Nun ist 


(W—p+ 2424 O)U+ Q—it)” 28 Q—i—2\np+atip ppa—ip 


ae i 1 a+ if (Seek ee ‘ 1 1 
268 Q—if—@\uta+tip eae) 2S? eee 


1 i pear | 1 1 


Benutzt man ferner die quasiasymptotische Form des Integrals (IT, 12) 


oe HA +i Ag, ; : : 
| pane ap Ue HP) Bs (Ag — Hs) (0 +t 18) 
0 
1 1 


! oh (¢ Agy— Hgy) (& + iB) (i Ag, — Hy)? (4 + if)?’ 
so wird mit 2? = 9? + f2 und 246 =2 2—1 


co 


; | eH ix eittAn du 
<0 
; (W— pt 2pQ2+ 2) (u+ Q—if) 


A 1 Ue es re 1 1 
ol Q +o=%) Bs, Hi, 


i (30 =) (13 ry 1 1 1 
Gasig sats Qt to + (Qe (Age Hey 
Fal o 1 8M Q— 2iFQ+ it 1 
OO lh ee 24 (Q—iop (Age 


Dieses ist bereits das zu beweisende Ergebnis. 

7. Zahlenbeispiele von Einzelfliigeln. Bei den vier folgenden Beispielen handelt es sich um 
Fliigel mit dem Plattenprofil f’(«) = konst. = — 0,1. Ferner sei fiir 

Beispiel la: ag/ui=0,1, eee 051. ; 2 = ugie ==2or, 

Beispiel lb: ag/ui=—0,1, hg/uz=0,1, Q=owlg = 20, 

Beispiel 3a: ag/u2=0,05, hg/uz=0,05, Y2=wul/g=4z, 

Beispiel 3b: ag/uis=0,05, hg/uy=0,05, 2 =wu/g = = huge 
Die Berechnung des stationiren Anteils der Beispiele la, 1b baw. 3a, 3b wurde bereits als 
Beispiel 1 bzw. 3 in Arbeit I, Ziff. 9 durchgefiihrt, und alle wesentlichen Gré8en, insbesondere 


die stationadren Mittelwerte der Fliigelzirkulation A, (in Arbeit I mit J’ bezeichnet) und der Wasser- 
oberflache, kénnen von dort unverindert ibernommen werden. 


Der Frequenzzahlwert Q = 22 baw. Q =4 bedeutet, daB der Tragfliigel in der Zeit = 
einer Wellenperiode mit der Geschwindigkeit uy die Strecke 10 @ zuriicklegt; bei den niedrigeren 


Werten 2 = = und 2 = e hetrigt diese Strecke 30 a. Bei jeder Frequenz w der wellenférmigen 
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Anstrémung ist dann gemaB (II, 3) eine kiirzere (Wellenzahl ,) und eine langere Welle (Wellen- 
zahl (2) méglich. Es ergibt sich fir* 

Beispiel la: uw, = 9,339, py = 4,227, == _ D183 ot ee otO0n 

Beispiel lb: uw, = 4,126, uw, =—1,063, 0. ==(_ Lobe. Bisa goes 

Beispiel 3a: uw, = 16,646, [lg = 9,486 , ow = 12,060, =p = 3,009; 

Beispiel 3b: uw, = 6,796, fly = 2,982, = 3,689, 6 = 1,985. 

Fiir die Auswertung der Theorie und die Diskussion unserer Beispiele wollen wir davon aus- 
gehen, da8 die von der Wellenanstrémung geformte Wasseroberflache weit vor dem Tragfliigel 
(x > — 00) gegeben ist durch 


+ ¥(x, 1) = By sin (» ii, af) + BY sin ( $y =) (28) 
0 40, 
GemaB (II, 39) ist dann 
B B 
BE =—2(u,—2) 7» By = —2(w,— 22) 7 


Da die Berechnung der Fourierkoeffizienten b,, fiir den Kernanteil k(x, €) der Integralgleichung 
(II, 27) sehr miihsam und langwierig ist, wurde fiir die Berechnung des instationaren Zirkulations- 
koeffizienten A, hauptsichlich die 1/,—*/,-Punkt-Methode verwendet. Um eine Vergleichsméglich- 
keit zu haben, wurde das Beispiel3b auch nach der genaueren Theorie behandelt, d.h. durch 
Auflésung der Integralgleichung (II, 27). Wie schon friiher in der stationaren Theorie hat sich 
gezeigt, daB die !/,—*/,-Punkt-Methode recht brauchbare Naherungswerte liefert; diese sind besser, 
wenn die Tauchbewegung des Fliigels mit beriicksichtigt wird, als bei festgehaltenem Fiigel. 

In Tabelle 1 sind die nach der !/,—*/,-Punkt-Methode berechneten Zirkulationskoeffizienten 4, 
in Einheiten von au, zusammengestellt. Fiir Beispiel 3b sind die Fourierkoeffizienten b,, der 
Reihenentwicklung (II, 32) in Tabelle 2 enthalten. Da die b,, wesentlich kleiner als 1 sind, kann 
das Gleichungssystem (II, 36) fiir die Koeffizienten P, der Lésung (II, 35) leicht durch Iteration 
aufgelést werden. Das Ergebnis ist in Tabelle 3 enthalten. (P, in Einheiten von uy.) Wie zu 
erwarten war, sind die 4,-Werte mit Beriicksichtigung der Tauchbewegung merklich kleiner als 
diejenigen ohne Tauchbewegung bei festgehaltenem Fligel; denn die Tauchgeschwindigkeit tibt 
ja in der Randbedingung einen stabilisierenden Einflu8 aus. 


Tabelle 1. 
Beispiel A, mit Tauchbewegung A, ohne Tauchbewegung 
la | —(0,072 + 0,038 i) BY + (0,067 + 0,061 i) BX | — (0,146 — 0,019 i) B¥ + (0,164 + 0,022 i) BE 
lb | — (0,022 + 0,040 i) B¥ + (0,011 + 0,030 1) B* | — (0,223 — 0,056 i) B¥ + (0,157 — 0,015 i) B¥ 
3a — (0,080 + 0,022 7) B* + (0,081 + 0,039 iz) BF — (0,106 — 0,009 7) B¥ +. (0,114 + 0,011 i) BF 
3b | — (0,046 + 0,052 i) B¥ + (0,030 + 0,043 i) Bk | — (0,152 — 0,042 i) B*¥ + (0,119 — 0,020 i) Be 


Tabelle 2. bay. 


2 
1 + 0,575 + 0,1707| + 0,275 — 0,154 7 | — 0,016 — 0,017i | — 0,013 + 0,001i 
2 — 0,122 + 0,076 i + 0,032 + 0,0347 | + 0,044 — 0,002 7 | — 0,001 — 0,002 i | — 0,003 0,000 
3 — 0,007 — 0,007 i | — 0,038 + 0,002 i | ++ 0,001 + 0,003 i + 0,009 0,000 | —0,001 
4, + 0,003 — 0,001 — 0,0017 — 0,007 + 0,001 7% + 0,003 0,000 
5 0,000 + 0,001 0,000 — 0,002 0,000 | + 0,001 
Tabelle 3. 
mit Tauchbewegung ohne Tauchbewegung 
A, |— (0,0451 + 0,0511 i) BY + (0,0305 + 0,0438 i) BY | — (0,1526 — 0,0453 i) BY + (0,1232 — 0,0215 i) BE 
P, |—(0,0275 + 0,0360 i) BF + (0,0175 + 0,0214 i) B¥ | —(0,0657 — 0,0131 i) B¥ + (0,0517 —0,0135 i) Bt 
P, |— (0,0094 + 0,0163 i) BF + (0,0056 + 0,0043 i) B* | — (0,0053 + 0,0048 i) B¥ + (0,0035 — 0,0051 i) Bs 
Ps |+ (0,0002 — 0,0008 i) BF + (0,0004-+ 0,0003 i) BF | —(0,0005 — 0,001 i) Bt + (0,0009 — 0,0002 i) Bt 
P, 0,000 BY + 0,000 BY 0,000 B¥ + 0,000 Bs 
P; 0,000 B¥ + 0,000 BF 0,000 B¥ + 0,000 Bs 


1 Fiir die Durchfihrung der in dieser Arbeit enthalten ischen R i 
Aare oe ae eee altenen numerischen Rechnungen danke ich Frau M. 
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Fir die weitere Diskussion und die zahlenmaBige Darstellung der Ergebnisse ist es zweck- 
maBig, bei jedem Beispiel die zwei charakteristischen Hauptfalle zu betrachten: 


Fall 1) 


Be er, 


Be =0, 


Fall 2) 


BY =0, 


fea. 
Bye e. 


Aus 1) und 2) kénnen dann durch lineare Superposition beliebige weitere Falle zusammengesetzt 


werden. 


Tabelle 4 zeigt die Funktion [(t)— A, d.h. den instationaren Anteil der Fligelzirkulation 
in Einheiten von au,e. Bei den w t-Werten kann man sich auf den Bereich bis beschranken, 


denn nach unseren Ansatzen fiir die Fligelzirkulation ist ja ( + ae Ay =— I(t) +A). 


Tabelle 4. I(t) — Ay. 


Beispiel la Beispiel 1b Beispiel 3a Beispiel 3b 
Fall 1 Fall 2 Fall 1 Fall 2 Falll |  Fall2 Fall 1 Fall 2 
+ 0,022 
eaat n/6 | —0,087 | + 0,055 | + 0,002 | —0,011 | —0,117 0,101 | —0,028 | + 0009 
a ache m/3 | —0,006 | —0,039 | + 0,047 | —0,041 | —0,042 | + 0,013 | + 0,044 | —0.044 
bewegung m/2 | + 0,076 | —0,122 | + 0,080 | —0,060 | + 0,044 | —0,078 | + 0,104 | — 0,086 
2 22/3 | + 0,138 | —0,173 | + 0,091 | — 0,063 | + 0,118 | —0,149 | + 0,136 | —0,104 
52/6 | + 0,163 | —0,177 | + 0,078 | —0,049 | + 0,161 | —0,179 | + 0,132 | —0,095 
0 — 0,292 | + 0,328 | —0,446 | + 0,314 | —0,212 | + 0,228 | —0,304 | + 0,238 
ae m/6 | —0,272 | + 0,262 | —0,442 | + 0,287 | —0,193 | + 0,186 | —0,305 | + 0,226 
Pavchi: m/{3 | —0,179 | + 0,126 | —0,320 | + 0,183 | —0,122 | + 0,095 | —0,225 | + 0,154 
bewegung m/2 | —0,038 | —0,044 | —0,112 | + 0,030 | —0,018 | —0,022 | — 0,084 | + 0,040 
22/3 | + 0,113 | —0,202 | + 0,126 | —0,131 | + 0,090 | —0,133 | + 0,079 | —0,084 
52/6 | + 0,234 | —0,306 | + 0,330 | —0,257 | + 0,175 | —0,208 | + 0,221 | —0,186 


Tabelle 5. 0,Ky. 


Beispiel la 


Beispiel 1b 
Fall 1 Fall 2 


Beispiel 3a 
Fall 1 Fall 2 


Beispiel 3b 
Fall 1 Fall 2 


Die Kraftkomponenten K, und K,, nach Formel (II, 37) setzen sich aus drei Anteilen zusammen. 
Zunachst ein von ¢ unabhangiger stationarer Anteil, der zahlenmafig genau mit dem in Arbeit I, 
Ziff. 9 berechneten K,- bzw. K,-Wert tibereinstimmt. Sodann ein in € linearer Anteil 0,K, bzw. 

-0,K,; dieser ist in Einheiten von 9 u? a ¢ in den Tabellen 5 und 6 enthalten. Hierbei ist die Super- 
position der Falle 1) und 2) méglich. SchlieBlich tritt ein in ¢ quadratischer Anteil auf. Wie wir 
noch sehen werden, muS im Giiltigkeitshereich unserer Theorie stets e* < 0,1 sein. Es hat sich 
gezeigt, daB dann der in ¢ quadratische Anteil von K, etwa von der GroSenordnung 0,001 @ us a 
ist und somit in Anbetracht des grofen stationaren Hauptwertes von K, vernachlassigt werden 
kann. 
 Dagegen diirfte wenigstens der in ¢ quadratische und zeitunabhingige Zusatzanteil 0,K, zum 
Wellenwiderstand K, trotz seiner Kleinheit von gewissem Interesse sein, denn der Hauptwert 
des Wellenwiderstandes ist ja selbst klein. (vgl. Arbeit I, Ziff.9). 0,K, ist in allgemeiner Form 
[Fall 1) und 2) nicht superponierbar] in Einheiten von @Q uj a in Tabelle 7 enthalten. 
Tabelle 8 zeigt den Verlauf der fiir die Tauchbewegung bzw. die Tiefenlage des Fligels maB- 


gebenden Funktion = n(t) + 1 gemaB Gleichung (II, 25) in Einheiten von «. 
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Tabelle 6. 0,K,,. 


Beispiel 3b 


Beispiel la Beispiel 1b Beispiel 3a 
~ Fall 1 | Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fal) 1 Fall 2 
| 
o | + 0,013 | —0,016 | + 0,006 | —0,002 | + 0,015 | —0,019 | + 0,010 | —0,007 
x/6 | —0,004 —0,006 | —0,008 | + 0,003 | + 0,004 | —0,013 | —0,001 | —0,001 
mat x/3-| —0,019 + 0,006 | —0,019 | + 0,007 | —0,008 | —0,003 | — 0,012 | + 0,004 
Tauch- n/2 | —0,030 | + 0,016 | —0,025 | + 0,009 | —0,018 | + 0,008 | —0,020 | + 0,009 
bewegung | 97/3 | —0,033 | + 0,022 | —0,025 | + 0,009 | —0,024 | + 0,017 | —0,022 | + 0,012 
52/6 | —0,026 | + 0,022 | —0,018 | + 0,007 | —0,022 | + 0,021 | —0,019 | + 0,011 
o | + 0,008 | —0,016- | + 0,016 | —0,012 | +.0,011 | —0,017 | + 0,017 | — 0,016 
1/6 0,000 | —0,016 | + 0,011 | —0,009 | + 0,006 | —0,016 | + 0,013 | —0,008 
ae n/3 | — 0,007 | —0,011 | + 0,003 | —0,009 | —0,001 | —0,011 | + 0,006 | + 0,002 
tener. n/2 | —0,013 | —0,004 | —0,006 | —0,003 | —0,008 | —0,002 | —0,003 | + 0,011 
bewegung | 97/3 | —0,015 | + 0,005 | —0,013 | + 0,003 | —0,013 | + 0,007 | —0,011 | + 0,018 
52/6 | —0,014 | + 0,016 | —0,017 | + 0,011 | —0,014 | + 0,014 | —0,016 | + 0,020 
Tabelle 7. 
Beispiel 0,.Ky mit Tauchbewegung | 0,Kz ohne Tauchbewegung 
la | — 0,007 B¥2 — 0,007 B32 + 0,015 BY B* | — 0,006 Bt? — 0,012 BY? + 0,020 B* Bt 
1b — 0,004 B¥® — 0,001 B32 + 0,004 B¥ Bs | — 0,018 B¥2 — 0,010 B%? + 0,027 B* B¥ 
3a — 0,005 B¥® — 0,005 B%2-+ 0,010 B* Bt | — 0,004 B¥2 — 0,006 BS2-+ 0,010 B* Bs 
3b — 0,004 B%? — 0,002 B¥2 + 0,006 B¥ BY — 0,009 B¥? — 0,005 B¥? + 0,013 Bt Bt 
: 
Tabelle 8. a nt) +1. 
: Beispiel la Beispiel 1b Beispiel 3a Beispiel 3b 
@ 
Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 | Fall 2 Fall 2 
0 + 0,107 | — 0,099 10) 293-' | <= .0,147 40,080" e008 + 0,260 | —0,170 
nl6} + 0,064 | —0,041 —0,013 | + 0,073 + 0,036 | — 0,032 + 0,078 -| — 0,025 
n/3| + 0,004 | + 0,028 =, 315¢ |. ---0,272 + 0,013  —0,005 0,125 | 40,125 
n/2| —0,058 | + 0,090 — 0,533 | + 0,400 — 0,014" | ~--0,024 = 01294 2b B48 
27/3| —0,104 | + 0,128 —0,609 | + 0,420 — 0,037 ~| + 0,047 —0,385 | + 0,295 
52/6| —0,122 | + 0,131 — 0,521 | = 0,327 —0,050 | + 0,056 =3 0372-1) 0.269 


SchlieBlich ist in Tabelle9 und 10 der instationadre Anteil an Y (x, ros o(x) der freien 
a a 


Wasseroberflache berechnet fiir die drei Punkte x =—2a, 0, +2a in der Umgebung des 


Fligels enthalten. (In Einheiten von ¢.) Dabei ist = Y,(x) der bereits in Arbeit I berechnete sta- 


tiondre Mittelwert der freien Wasseroberflache (vgl. dort Tabelle 5, Beispiel 1 und 3). 


i 
Tabelle 9. Te: Y(x, t) — a Y,(x) (mit Tauchbewegung). 


Beispiel la 


Fall 1 


+ 0,954 
+ 0,050 
— 0,884 
+ 0,667 
a0. 581 
— 0,850 
++ 0,202 
++ 0,870 
— 0,587 
— 0,318 
+ 0,976 
— 0,168 
— 0,752 
+ 0,820 
+ 0,297 
— 0,985 
+ 0,445 
+ 0,682 


Fall 2 


+ 0,758 
— 0,048 


Beispiel 1b 


Fall 1 Fall 2 
+ 0,718 | + 0,224 
— 0,014 | + 0,012 
— 0,696 | — 0,232 
+ 0,954 | + 0,687 
+ 0,474 | + 0,519 
— 0,293 | + 0,310 
+ 0,934 | + 0,966 
+ 0,835 | + 0,887 
+ 0,189 | + 0,769 
+ 0,664 | + 0,986 
+ 0,972 | + 1,018 
+ 0,620 | + 1,022 
+ 0,216 | + 0,742 
+ 0,849.| = 0,875 
+ 0,885 | + 1,001 
— 0,290 | + 0,299 
+ 0,498 | + 0,499 
+ 0,913 | + 0,712 


Beispiel 3a 


Fall 1 Fall 2 
+ 0,998 | + 0,814 
+ 0,056 | —0,056 
— 0,938 | —0,902 
+ 0,804 | + 1,010 
+ 0,540 | + 0,466 
— 0,787 | —0,561 
+ 0,396 | + 0,935 
+ 0,880 | + 0,862 
— 0,426 | —0,070 
— 0,119: | +. 05610 
+ 0,984 | + 1.028 
+ 0,050 | + 0,440 


Beispiel 3b 
Fall 1 Fall 2 
+ 0,610 | + 0,272 
— 0,014 | + 0,014 
— 0,560 | — 0,300 
+ 0,902 | + 0,730 
+ 0,463 | + 0,530 
— 0,138 | + 0,258 
+ 0,954 | + 0,991 
+ 0,816 | + 0,904 
+ 0,321 | + 0,747 
+ 0,749 | + 0,987 
+ 0,950 | + 1,036 
+ 0,694 | + 1,036 
+ 0,344 | + 0,719 
+ 0,830 | + 0,890 
+ 0,881 | + 1,047 
— 0,154 | + 0,258 
+ 0,487 | + 0,506 
-+ 0,832 | + 0,778 
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Tabelle 10. = Nico — oa o(x) (ohne Tauchbewegung). 


Beispiel 1b 
Fall 1 


Beispiel 3a 
Fall 1 Fall 2 


i 


Beispiel 3b 
Fall 1 Fall 2 


+ 0,909 | + 0,928 | + 0,755 | + 0,855 


Fir eine kurze Diskussion der erhaltenen Ergebnisse ist es zunachst notwendig, sich wenigstens 
ungefahr uber den Giiltigkeitsbereich der Theorie klar zu werden. Die einschneidendste Verein- 
_fachung bei der vorgelegten Theorie liegt in der Tatsache, daB wir die Randbedingung der freien 
Wasseroberflache nur lings der Ebene y = 0 erfiillen konnten und nicht lings der genauen wellen- 
férmigen Oberflache. Durch diese Wellenform und die Tauchbewegung des Fligels. schwankt 
der Abstand des Fliigels von der freien Wasseroberflache periodisch, und diese Abstandsanderungen 
_ wurden in unserer Theorie vernachlassigt. Bei Untersuchungen im Rahmen der static naren Theorie 
der Arbeit [-hat sich gezeigt, daB (bei den Froudeschen Zahlen unserer Beispiele) eine Halbierung 
baw. Verdoppelung des Fliigelabstandes von der freien Wasseroberflache eine Zirkulationsanderung 
-der GréB®enordnung von etwa 15% bis 18% bedingt. Wir wollen daher als ungefahre Grenze fiir 
den Giiltigkeitsbereich der hier mitgeteilten Ergebnisse folgendes festsetzen: Ist h der mittlere 
Abstand des Fliigels von der freien Wasseroberflache, so soll der Maximalabstand bei den peri- 
_ odischen Schwankungen nicht gréSer als 1,25 h und der Minimalabstand nicht Kleiner als 0,75 h 
sein. Dann ist nach den in der stationéren Theorie gemachten Erfahrungen zu vermuten, dab 
der Fehler bei den Ergebnissen des instationaren Anteils der Theorie in der GroBenordnung von 
10% bleibt. Eine genaue Aussage iiber diesen Fehler la8t sich natiirlich erst machen, wenn es 
gelingt, eine im obigen Sinne exakte Theorie zu entwickeln. Dieses ist bisher nicht gelungen, 
da die mathematischen Schwierigkeiten zu groB waren. 

Aus Tabelle 9 und 10 erkennt man, daB die Amplituden des instationaéren Anteils der freien 
‘Wasseroberflache tiber dem Fliigel bei allen Beispielen etwa von der GréBenordnung ¢ h sind 
(a =h). Die Amplitude der Tauchbewegung (Tabelle 8) fallt demgegeniiber bei den Beispielen la 

und 3a nicht ins Gewicht; fiir Beispiel 1b und 3b ist sie jedoch mit zu beriicksichtigen; und zwar 
- wirkt die Tauchbewegung im Fall 1) im Sinne einer Verstarkung der zeitlichen Abstandsschwankung 

Fliigel-Wasseroberflache und im Fall2) im Sinne einer Verminderung. Der Giitigkeitsbereich 
-der erhaltenen Ergebnisse liegt also gemaf unserer obigen Festsetzung fiir die Beispiele ohne 
‘Tauchbewegung und Beispiel la und 3a mit Tauchbewegung bei ¢ < 0,25. Dagegen kann fiir 
Beispiel 1b und 3b mit Tauchbewegung im Fall 2 (d.h. gréBere Wellenlainge der Anstrémung) 
noch ¢ = 0,35 zugelassen werden, im Fall 1 (d. h. kKleinere Wellenlange der Anstrémung) héchstens 
be — 0,18. 

Entsprechend der wellenformigen Anstrémung ist fiir den gleichen Zeitpunkt (d.h. bei einer 
Momentaufnahme) die Verformung der Wasseroberfliche im Fall 1) stets steiler als im Fall 2). 
Auch hier sind natiirlich der Giiltigkeit der Theorie insofern Grenzen gesetzt, als Anstroémungs- 

-wellen, bei denen das Verhaltnis Mb von Wellenlange zu Fliigeltiefe zu klein wird’ nur mit 
1,2 


1 Auch hier ist es natiirlich sehr schwer, eine genaue Grenze anzugeben, da es eine Theorie mit Beriick- 

sichtigung der Kriimmung der Wasseroberflache in der Randbedingung bisher nicht gibt. Auf jeden Fall muB 
man jedoch prinzipiell die genannte Einschrankung im Auge behalten, wenn man die Zuverlassigkeit der 
_ Ergebnisse beurteilen will. Vgl. auch Arbeit I, Ziff. 6. 
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noch verkleinerten Amplituden behandelt werden kénnen. Hier ist fir 


2 2 

Beispiel la: “2 2 — 3,36, “% 7 — 7,43; . Beispiel lb: “2 — 7,62, = — = 29,55, af 

SUL YEA ff [bp § a My ECD [eb | 

2 2 

Beispiel 3a: “4 — 3,77, % ~— 6,62;  Beispiel 3b: “° = 9,25," — = 2430.8) 
a fy BO hg Ba fy § & Me 


Bezogen auf den stationiren Mittelwert 4, betragen die instationdren Schwankungen der Fliigel- 
zirkulation gemaB Tabelle 4 (dabei wurde iiberall ¢ = 0,25 gesetzt): 


mit Tauchbewegung 


Beispiel la: 1) 12%, 2) 13%; Beispiel lb: 1) 7%, 2) 5%, 
Beispiel 3a: 1) 10%, 2) 11%; Beispiehi3b<. 1) 11%, 2)01%s 


ohne Tauchbewegung 


Beispiel la: 1) 22%, 2) 24%; Beispiel lb: 1) 33%, 2) 23%, 
Beispiel 3a: 1) 13%, 2) 14%; Beispiel 3b: 1) 19%, 2) 15%. 


Man erkennt daraus wieder den stabilisierenden Einflu8 der Tauchbewegung; und zwar ist dieser 
um so ausgepragter je kleiner die Frequenzzahi (2 ist. Entsprechend ist ja auch die Tauchbewegung 
gema8 Tabelle 8 am starksten fiir Beispiel 1b und am schwachsten fir Beispiel 3a. 

Ein analoges Verhalten wie die Fliigelzirkulation zeigt gem4B Tabelle 5 auch die Auftriebs- 
kraft K,. Kleine Unterschiede zwischen J’ und K, sind durch die zur Hauptanstrémung uy hinzu- 
tretenden Stérgeschwindigkeiten bedingt [Vgl. (II, 37)]. Auch der Wellenwiderstand K, hat gema8 
Tabelle 6 gewisse periodische Schwankungen, auf deren GréBe die Tauchbewegung aber keinen 
merklichen Einflu8 ausibt. 


8. Zahlenbeispiele von zwei hintereinander fahrenden Fliigeln. Auch bei den vier folgenden 
Beispielen handelt es sich um Fligel mit dem Plattenprofil f{(x,) = f3(%.) =— 0,1 = konst. und 
a, =a, =a. Ferner sei fiir 


Beispiel 4a: ag/ui—0,1, h,g/ue—0,1, hag/uz=—0,05, Q=22, L=20a, MB/O—10. 1 
Beispiel 4b: ag/uz=0,1, hyg/ug—01, hyg/u?=0,05, Q—=lx, L=20a, MLO =10. 


Beispiel 5a: ag/u}=—0,05, h, g/ug =0,05, h,g/u{=0,025, Q=4n2, L=20a, ML-JO=10. 
Beispiel 5b: ag/ug =0,05, h, g/u2 =0,05, hog/u2 = 0,025, Q=4n, L=20a; ML?/O=10. 


Die Berechnung des stationaren Anteils der Beispiele 4a, 4b baw. 5a, 5b wurde bereits als Beispiel 4 
bzw. 5in Arbeit I, Ziff. 10 durchgefiihrt, und alle wesentlichen GréBen, insbesondere die stationaren 
Mittelwerte der Fliigelzirkulationen Af), AP (in Arbeit I mit [,, [, bezeichnet) und der Wasser- 
oberflache, konnen von dort unverandert ibernommen werden. 

Die Werte von (4, lg, %, 6 sind fiir die Beispiele 4a, 4b, 5a, 5b die gleichen wie fiir die Beispiele la, 
lb, 3a, 3b aus Ziff. 7. 

Bei der Auswertung der Theorie gehen wir auch von der Formel (28) aus, die die Form der |} 
freien Wasseroberflache weit vor den Fliigeln wiedergibt, und betrachten bei jedem Beispiel die | | 
zwei charakteristischen Hauptfalle 


Falll):. Bi =e, Br =0, Fall2):. BY = 0 BP se. 


Aus 1) und 2) kénnen durch lineare Superposition beliebige weitere Falle zusammen gesetzt werden. | 

Fir die Berechnung der instationdren Zirkulationskoeffizienten wurde die 1 /4—3/4-Punkt- | 
Methode verwendet. Dabei haben wir die Zwei-Fliigel-Systeme nur mit Beriicksichtigung der 
Stampf- und Tauchbewegung durchgerechnet, da diese bei Tragfliigelbooten gerade wesentlich ist. 
Bei festgehaltenen Fligeln hatte sich ohnehin fiir Fliigel ] genau das gleiche Resultat ergeben wie || 
bei den entsprechenden Einzelfliigeln ohne Tauchbewegung. | 

Tabelle 11 zeigt die Funktionen [j(t)— A® bzw. I(t) — A® d.h. den instationa’ren Anteil |} 
der Fligelzirkulationen in Einheiten von a uy ¢. Die Tabellen 12 und 13 geben den Verlauf der fiir | 


die Tauch- und Stampfhewegung der Fliigel maBgebenden Funktionen — 7(1) — 2 d in #| | 
Einheiten von ¢ wieder. Vgl. Formeln (14) und (18). ag ME) a Mound oe) na 


1 Der Wert M L?/© = 10 entspricht einem mittleren bei Tragfliigelbooten vorkommenden Wert. 
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: peer x) al I 
Der instationare Anteil = Y(x,t)—— Y,(x) der freien Wasseroberflache ist in Tabelle 14 ent- 


_halten berechnet fiir die drei Punkte « — —2 a, 0, + 2 ain der Umgebung des Fliigels J, und in 
Tabelle 15 fiir die drei Punkte x = 18 a, 20 a, 22 ain der Umgebung des Fliigels 2. (In Einheiten 


von €.) Dabeiist = Y,(x) der bereits in Arbeit I berechnete stationare Mittelwert der freien Wasser- 
oberflache (vgl. dort Tabelle 7*, Beispiele 4 und a) 

Fiir den Geltungsbereich unserer Ergebnisse erhalten wir bedingt durch die Strémungsverhalt- 
nisse bei Fliigel 2 ein ungiinstigeres Bild als bei den Einzelfliigeln. Zundchst ergibt sich aus Ta- 
belle 14 und 15, daB die Amplituden des instationdren Anteils der freien Wasseroberflache iiber 
Fliigel 1 baw. Fliigel 2 von der GréBenordnung ¢ h, baw. 2 ¢ hy sind (a = h, =2h,). Weiter sind die 
Amplituden der Tauch- und Stampfhewegung zu beriicksichtigen. Die durch die Tauch- und 
Stampfhewegung bewirkte Auslenkung gegeniiber der stationaren Mittellage ist bei Fliigel J gegeben 


durch = [n(t) — yo — L, p(t)] und bei Fliigel 2 durch = [n(t) —. + Ly v(t)]. 

Fir Beispiel 4a und 4b ist nach (16) L, =10,5a, L, =9,5 a, und fiir Beispiel 5a und 5b ist 
i= 31 a,b, = 11,9 a. 

Man erkennt: Relativ unerheblich fiir die zeitlichen Schwankungen des Abstandes Fliigel- 
Wasseroberflache ist die Tauch- und Stampfbewegung bei den Beispielen 4a und 5a mit den groBen 


Q-Werten, wahrend sie bei den anderen Beispielen teilweise groBe Bedeutung hat. Und zwar wirkt 
die Bewegung des Fliigels J im Fall 1 verstarkend und im Fall 2 teils verstarkend teils abschwachend 


Tabelle 11. 13) — AY. 


Beispiel 4a 


Beispiel 4b Beispiel 5a Beispiel 5b 


Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 
—0,088 | + 0,052 | 

n/6 | —0,084 | + 0,063 | —0,063 | + 0,058 | —0,056 | + 0,056 | —0,052 | + 0,060 

Fliigel 1 n/3 | + 0,001 | + 0,016 | —0,020 | + 0,048 | + 0,009 | —0,002 | + 0,026 | + 0,024 
(Gj=1) n/2 | + 0,086 | —0,035 | + 0,027 | + 0,025 | + 0,072 | —0,060 | + 0,098 | —0,018 
2n/3 | + 0,148 | —0,078 | + 0,068 | —0,004 | + 0,115 | —0,101 | + 0,143 | —0,055 

52/6 | + 0,170 | —0,099 | + 0,090 | —0,032 | + 0,128 | —0,115 | + 0,150 | —0,077 

0 £20,157 | 4° 0,064 | —0,113 | + 0,069 | +.0,119 | — 0,086 4 0,097 |--+ 0,047 

af6— | 0,077 | + 0,140 | —0,122 |°+-0,072 | + 0,129 | —0,026 | —0,057 | + 0,068 

Fliigel 2 nfs | + 0,024 | + 0,179 | —0,098 | + 0,055 | + 0,104 | + 0,041 | —0,002 |. + 0,070 
(j = 2) nj2 | +0,119 | + 0,170 | —0,048 | + 0,024 | +.0,051 | + 0,098 | + 0,053 | + 0,054 
2n/3 | + 0,182 | + 0,115 | + 0,015 | —0,014 | —0,016 | + 0,128 | + 0,095 | + 0,024 

52/6 | + 0,196 | + 0,030 | + 0,074 | —0,048 | —0,078 | + 0,124 | + 0,111 | —0,013 


1 1 
Tabelle 12. ae n(t) — ae 


Beispiel 4a 


Beispiel 4b Beispiel 5a Beispiel 5b 


oh Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 
— 0,382 — 0,002 
m/ 6 + 0,057 — 0,072 0,583 — 0,409 — 0,014 — 0,006 + 0,184 — 0,213 
m/3 — 0,009 — 0,069 + 0,374 — 0,326 — 0,021 — 0,007 — 0,040 — 0,159 
0/2 — 0,072 — 0,047 + 0,065 — 0,156 — 0,023 — 0,007 — 0,253 — 0,062 
2 2/3 — 0,116 — 0,013 — 0,261 + 0,056 — 0,019 — 0,005 — 0,399 + 0,052 
5 7/6 — 0,129 + 0,024 — 0,518 + 0,253 — 0,009 — 0,002 — 0,437 + 0,152 


Tabelle 13. g(é). 


Beispiel 4a Beispiel 4b Beispiel 5a Beispiel 5b 
ae Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 
90128 — 0,0091 —0,0106 | + 0,0085 + 0,0086 + 0,0036 
0405 — 0,0063 — 0,0092 | + 0,0036 + 0,0107 — 0,0134 
0574 | —0,0019 — 0.0054 | — 0,0024 + 0,0099 — 0,0268 
0589 | + 0,0030 — 0,0001 — 0,0077 + 0,0064 |. — 0,0330 
0446 | + 0,0072 + 0,0052 — 0,0109 + 0,0012 — 0,0304 
,0184 | + 0,0094 + 0,0092 — 0,0112 — 0,0043 — 0,0196 


15 
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1 1 
Tabelle 14. = Y(x, t) — a Yj(z)ic 


Beispiel 4b Beispiel 5a Beispiel 5b 


Beispiel 4a 


gs i Fall 1 | Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 | Fall 2 
, ; 
| + 0,748 hse 0,202 + 0,992 | + 0,815 + 0,613 + 0,254 
| — 0,032 + 0,010 | — 0,018 + 0,038 | —0,035 —0,007 | — 0,011 
+ 2 0 — 0,875 — 0,783 — 0,666 | — 0,249 —0,918  — 0,878 — 0,541 | — 0,316 
ae 7/6 + 0,666 + 0,987 | + 0,963 | + 0,669 | + 0,802 + 1,006 | + 0,903 + 0,713 
0 1/6 + 0,531 + 0,477 + 0,493 | + 0,490 + 0,520 + 0,480 -- 0,466 | + 0,505 
+ 2 /6 — 0,841 | — 0,387 — 0,246 | + 0,262 — 0,796 | — 0,512 —0,119 | + 0,218 
—2 n/[3 + 0,202 | + 0,962 -- 0,935 | + 0,957 + 0,398 | + 0,928 + 0,950 | + 0,982 
0 r/3 + 0,868 | -+ 0,859 --+ 0,844 | + 0,866 + 0,863 | + 0,867 + 0,814 | + 0,885 
+ 2 m/3 — 0,581 | + 0,113 + 0,241 | + 0,703 — 0,461 | — 0,008 + 0,335 | + 0,693 
—2 m/ 2 — 0,319 | + 0,679 + 0,656 | + 0,988 —0,113 | + 0,600 + 0,743 | + 0,987 
0 me/ 2 + 0,972 | + 1,011 + 0,969 | + 1,010 + 0,975 | + 1,021 + 0,943 | + 1,028 
j59 n/2 — 0,166 | + 0,582 | + 0,663 | + 0,955 | —0,002 | + 0,497 | + 0,699 | + 0,983 
—2 2 1/3 — 0,752 | + 0,214 + 0,202 | + 0,755 — 0,594 | + 0,113 + 0,338 + 0,728 
0 2 7/3 + 0,817 | + 0,892 + 0,834 | + 0,884 + 0,826 | + 0,901 + 0,820 | + 0,896 
ais9 22/3 | + 0,294 | + 0,895 | + 0,907 | + 0,952 | + 0,457 | + 0,870 | + 0,875 | + 1,009 
—2 5 7/6 — 0,985 | — 0,308 — 0,306 | + 0,319 — 0,916 | —0,406 — 0,159 | + 0,273 
0 5 2/6 + 0,442 | + 0,534 + 0,476 | + 0,521 + 0,455 | + 0,540 + 0,478 | + 0,523 
+ 2 5 70/6 + 0,675 | + 0,969 + 0,908 | + 0,693 + 0,794 | + 1,009 + 0,818 | + 0,765 


Beispiel 4a Beispiel 4b Beispiel 5a Beispiel 5b 
Fall 1 | Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 Fall 1 Fall 2 


auf die durch die wellenférmige Wasseroberfliche bedingten Abstandsschwankungen ein. Die Be- 
wegung des Fliigels 2 verstarkt im Fall 1 und vermindert im Fall 2 die Abstandsschwankungen. 

Das fiir den Giiltigkeitsbereich der vorliegenden Theorie ungiinstigste Verhalten zeigt Beispiel 4b 
Fall1: Hier betragen die Abstandsschwankungen zwischen Fliigel 2 und der Wasseroberflache 
maximal etwa + 4h,e. Schon fiir ¢ = 0,25 wiirde der Fliigel also die Wasseroberflache beriihren, 
und dieses ist nach den unserer Theorie zugrunde liegenden Voraussetzungen natiirlich ganz unzu- 
lassig. 

Wenn wir fir den Giiltigkeitsbereich der Ergebnisse die gleichen Grenzen festlegen wie bei den 
Einzelfliigeln (vgl. unsere Uberlegungen auf S. 215), so erhalten wir ungefahr folgende obere Schran- 
ken fiir die e-Werte: Beispiele 4a und 5a¢ w 1/8, Beispiel 4b Falll ¢ ~ 0,065, Fall2¢ @ 0,21. 
Beispiel 5b Falll ¢ 0,10, Fall2 ¢ ~ 0,20. 

Dabei miissen wir uns an den Verhiiltnissen bei Fligel 2 als den ungiinstigeren orientieren. 

Tabelle 11 zeigt, da8 die Amplituden der instationaren Schwankungen der Fligelzirkulation 
ungefahr von der gleichen Gré®enordnung sind wie bei den entsprechenden Einzelfligeln mit Tauch- 
bewegung. Das gleiche gilt fiir die Tauchbewegung selbst; diese ist wieder bei Beispiel 4b (mit der 
kleinsten Frequenzzahl Q) am stirksten ausgepragt, wahrend sie fir Beispiel 5a (gréBter Q-Wert) 
fast ganz verschwindet. Die durch Tabelle 13 gegebene Stampfbewegung hat nur scheinbar eine 
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kleinere GréBenordnung als die Tauchbewegung; denn die g-Werte miissen ja noch mit den Hebel- 


_armen L, bzw. L, multipliziert werden, die von dex GréBenordnung 10 a sind. 


Trotz ihres bereits recht komplizierten Aufbaues stellt die instationire Theorie der vorliegenden 
Arbeit sowie der Arbeit II noch eine erhebliche Vereinfachung des Problemes dar. Fiir eine genaue 
Behandlung mii®te die Randbedingung (IT, 2) lings der wellenférmigen Wasseroberflache erfillt 
und die Schwankung des Fliigelabstandes von der Wasseroberflache beriicksichtigt werden. Au®er- 
dem ware unter Umstinden eine genauere Erfassung der Tauch- und Stampfhewegung der Fliigel 
notwendig, z. B. mit Beriicksichtigung der Dampfungswirkung des Wassers. Hierbei sei jedoch 
darauf hingewiesen, daB Weinblum’ in einer tiberwiegend schwingungsmechanischen Untersuchung 
(ohne hydrodynamische Beriicksichtigung der freien Wasseroberflache und der freien Wirbel) tiber 
das Tauchen und Stampfen von Tragflachensystemen in regelmaSigem Seegang im wesentlichen 
dieselben Bewegungsgleichungen verwendet wie wir. Diese haben sich fiir Stabilitatsrechnungen 
bewahrt, wie eine neue Arbeit von Kaplan und Jacobs? zeigt, die auf den Weinblumschen Glei- 
chungen beruht. In einer weiteren Arbeit? werden diese Gleichungen dadurch erginzt, daB die 
instationaéren Kraftanteile nicht mehr quasistationar nach der Methode der Auftriebsbeiwerte 
berechnet werden}, sondern mit den Formeln der instationaren Tragfliigeltheorie, allerdings ohne 
jede Beriicksichtigung der freien Wasseroberflache. 


Aus dem Institut fiir Angewandte Mathematik und Mechanik der Deutschen Akademie 
der Wissenschaften zu Berlin 


(Eingegangen am 14. Juni 1960.) 
Anschrift des Verfassers: Prof. Dr.-Ing. habil. W. H. Isay, Berlin-Dahlem, Schweinfurthstr. 90 


1 Nach G. Weinblum, Schiffstechnik 5 (1958) S.2 (aber in unserer Bezeichnungsweise) lauten die Bewe- 
gungsgleichungen fiir das Tauchen und Stampfen, wenn ein Austauchen der Tragflachen-aus dem Wasser 
nicht zugelassen wird, 


my = K,+ K,+ K,4+ K,, 09=M,+ M,+ M,+ M,; 
mit (a, Ss sind die Auftriebsbeiwerte der Fliigel) 


Q v Paw Q oh pad 
Ka pubes Bb ac, i), Kamp (ine, Trea, 2), 
> Ls¢ , Le4 , , 
Ky= gus (2a, 2 — 2a, ¢, 20), K,=— = uj (2 4, 67,9 + 242 ¢,, 9), 
Mie? (2a eT, — 2 a, ta, Dp Mo= ila, 1, = eb 
1 0 1%, 11 a fel oe 2— 9 Uo 1&q, $4 2Cq, Le > 
L236 , L234 é 
My=— Fu (Bae, e+ Bac, oe M,=+u§ (2a, ¢, Lip —2a,¢, 129) . 


Dabei ist K, die anregende Kraft der Wellen mit v,/u, und v,/uy als den an Fliigel 1 und 2 wirksamen Anstrém- 
winkeln der Wellengeschwindigkeit. Kg ist die stabilisierende Kraft, die durch die Tauchgeschwindigkeit 
bedingt ist. K,ist die durch die Stampfgeschwindigkeit an den Fliigeln induzierte Kraft. Die Kraft K, schlief- 
lich ist durch den infolge der Stampfbewegung geanderten Neigungswinkel der Fliigel bedingt; dabei soll fur 
die zeitliche Mittellage p = 0 angenommen werden. M, bis M, sind die den Kraften K, bis K, entsprechenden 
Drehmomente. Man erhalt so fur die Tauch- und Stampfbewegung die folgenden simultanen Differential- 
gleichungen: 
mj + Q Uo (41 4, + 42 6,) 1+ Q Uo (— 46, Dy + a,c, Ls) Y 

+ UG (464, + Ae C1.) P = 0 Uy (% C, U1 ++ gc), V2) » 


OP + Ou (4 cy, Li + age, L3) + QU (—% ly + dy c, L2) a7] 
+ 9 ug (— a6, Ly + a,c, Le) p = @ Uo (— a, 6, Ly 0 + az ¢, Le v2) « 
Fiihren wir nun fiir den instationaren Anteil der Fligelzirkulation die Bezeichung gema8 Formel (13) ein, so 


wird entsprechend dem aus der Methode der Auftriebsbeiwerte bekannten Zusammenhang zwischen Zirku- 
Jation und Anstrémwinkel 4 


AD Oho AW e108 ay C (vy — 4% + L, @ — % ) ? 
MAY FOF a AP) +? — ay cf, (¥3—) — Lap — U9) - 


Durch Einsetzen dieser Relationen in das obige Differentialgleichungssystem ergeben sich genau unsere in 


Ziff. 3 abgeleiteten Bewe sgleichungen. ; 
= P. Kaplan u. W. Toa Stevens Institute of Technology, Davidson Laboratory Report 704 (1959), 


Hoboken, N.J. : 
3 F. Ogilvie, J. ship research 3 (1959/60), Heft 3. 
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Ein neues Variationsprinzip in der Elastizitatstheorie* 
Von D. Riidiger 


1. Einleitung. Im folgenden wird ein Extremalprinzip fiir einen anisotropen, elastischen Kérper 
vom Volumen R und der Oberfliche F (Abb. 1) angegeben. Dieser wird auf allgemeine Koordinaten 


x* baw. x#, x?, x¢ (A, 4, ¥,0 = 1, 2, 3) 


bezogen. Normalenvektor der Oberflache in kontra- und kovarianter Darstellung ist n’, n,. Das 
anisotrope Medium hat den Elastizitatstensor 


E4uve — Ereku — ForvAn — Forres 


A 
70; Ty, 


Volumen: R 
Volumenelement: av 


Qberfache : FF 
Obertltichenelement: af 


Abb, 1, Elastischer Kérper. 
mit den angegebenen Symmetrieeigenschaften. Das Verschiebungsfeld in kontra- und kovarianter 
Darstellung ist 
VA (ot, 22, 03) Va x1, x3, x9).. 
Es gelten das Summationsiibereinkommen 


ae 
2,4 B= Ai B, 


und die kovarianten Ableitungen 


av? Hoe ahy aVa , 
fa eae peers aa ee 2 iesrareed le 


Hierin sind J” uy die Christoffelschen Dreiindizessymbole zweiter Art. Als Abkiirzungen werden 
der Spannungstensor PSS BEY | 

der Spannungsvektor der Oberflache # = Ee V,|, nj, 

die Divergenz des Spannungstensors p! = E4uve Valea 


und die Formanderungsarbeit A= 2 | | | Eure Vala Vol, dt 
R 


eingefiihrt. AuSerdem wird der Gaufsche Integralsatz 
SS hae Jf mat 
verwendet. 


* Vortrag, der auf dem Zehnten Intern. Kongr. f. Angew. Mechanik in Stresa (1960) gehalten wurde. 
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2. Prinzipe von Green und Castigliano. Es sei eine allgemeine Randwertaufgabe der Grund- 
gleichungen in der Elastizitatstheorie kleiner Verschiebungen 


E¥ere Vijjq + pt = 0 


vorgegeben; oe fen! sollen im betrachteten Raumteil R die Neg TE a) den 1 


oberflachen G, CG, G die Flachenkrifte A, B, # und auf den Restflachen ee if — G G = ff — C 


Z 2 3 
Teilfltiche G leilfliche G Teilflache G 


heilfliiche G 


Teilfliche G 


v v 
Abb. 2. Teilflachen G und G. 


a. a == = 
G = F—G die Verschiebungen V,, V,, V; vorgegeben sein (Abb. 2). Dieser Integrationsaufgabe 
sind bekanntlich die Variationsprinzipe von Green! 


U=A—Jffpe V,dr—ffev, df—ffev, df—Jff ® V,df = Extrem. 
; G G c 
[6Vj]a =0, [OV.]2 =0, Tale 2 
G 

und von Castigliano 

eee eee” df = Extrem. 

a 
[or], =0, [d?7]2 =0, [do], =0, [dp"]z = 0 
G G G 


aquivalent. 


3. Variationsprinzip. Es wird behauptet, da8 durch Hinzufiigen von 


2 — [he — aio. —Ma— (fe Aa 


zum Greenschen Hee und von 


—Z={ f (e—a)na+ [ {ea nat f {nat f [lop vad, 
é : : : 


zum Castiglianoschen Prinzip ein Variationsprinzip 
3 * * 
M20+Z=—(U +2) Extrem. 


der Elastizitatstheorie erhalten wird, an dessen benachbarte Verschiebungsfelder hinsichtlich der 
_ Randbedingungen keine F ord career mehr gestellt zu werden brauchen. Das Extremalprinzip 


lautet, wenn die Werte U, Z bzw. U, Z eingesetzt werden, 


M=U+Z=A— Dr 9 ATS a 
Shee V f= “100 LP) af als Pg =Baeem. 


1 Vgl. z. B. C. B. Biezeno u. R. Grammel, Technische Dynamik, 2. Aufl. Bd. 1 S. 80. Berlin/Géttingen/ 
Heidelberg 1953. 
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bzw. 


m=—(U+2)=—4+/ [fo —Pyy a+ 
R 


+{{e—ane+ffe-anag+] f@—8) nat 


G 


+ ffeva+fferat | [ ° Veaf = Extrem. (2) 


Die Funktion WV, Vos Ves A, 2, 8, Vij ee V) wird als potentielle Energie bezeichnet. Die 
Identitat der beiden Formulierungen des Extremalprinzips ist mittels der Formanderungsarbeit 


tad [[frraces [fora 


F 


R 
=z [[[mudts[[anars|fenars | [ene+ 
R 


ee 
G 


vs[fevarcs [formes [Jena 


Q|+ 
ale 


Allee 
lee 


a : 
ohne weiteres zu bestatigen. 
* * 2 
Setzt man diese namlich in JJ = U+ Z baw. IT = — (U o Z) ein, so ergibt sich als dritte 


Formulierung des Variationsprinzips 


T2037 Ss ye7) == ||| oP) V, dt + 


R 


+z [[@-2)na+ 5 [[@-2Anaets 
é é 


[@—2e) = 


— 


a|° 


-1ffatn—afa—2 fatn—aha—} [ [e(—2 Rapataren. 6) 


ale-— 


Ol 
Ol 


Sie ist fiir eine Anwendung der direkten Methoden der Variationsrechnung besonders gut ge- 
eignet’. Wenn man an die benachbarten Verschiebungszustande die Bedingung stellt, daB sie 
die geometrischen Randbedingungen erfiillen, dann entfallen in (1) die Integrale iiber die Teil- 


flachen G, und das neue Variationsprinzip geht in das Greensche iiber. Fordert man jedoch, daB 


die benachbarten Verschiebungsfelder die dynamischen Randbedingungen erfillen und mit der 
v 


Volumenkraft p” im Gleichgewicht sind, dann entfallen in (2) die Integrale iiber die Teilflachen G 
und tiber das Kérpervolumen R und // ist mit der Erganzungsarheit des Castiglianoschen Prinzips — 
identisch. 


z 2 3 dats 2 Se = = 
In dem Sonderfall G —G =G =G und G=G=G=F—G6 =GE geht (1) in ein von 
E. Reissner® angegebenes Extremalprinzip fiir eine spezielle dritte Randwertaufgabe iiber. 


1 D. Riidiger, Ing.-Arch. 27 (1960) S. 421. 

2 E. Reissner, Sympos. on the Calculus of Variations and its Applications of the American Mathematical 
Society, April 12, 1956. Vel. hierzu auch die im Falle groBer Verschiebungen giiltigen Prinzipe fiir ein Rand- 
wertproblem, bei dem auf der festen Teilflache G samtliche drei Komponenten des Spannungsvektors t? und auf | 


der festen Restflache G alle drei Komponenten des Verschiebungsvektors V’ vorgegeben sind: E. Reissner, — | 
4 5 sath. Phys. 32 (1953) S. 129; E. Koppe, Nachr. Akad. d. Wiss. Gottingen, math.-phys. Klasse, II a, (1956) 
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4. Variation der potentiellen Energie. Die Variation der potentiellen Energie JT nach (1) liefert 
zunachst 


oll = 0 = 64 —Jff pedV, de — 


ae = Sf PV, Often’ ae df— 


eo 


—1so%(r, —¥,) y—ffee(v.— Ve) df — ff 6 (V_—V,) af — 
¢ C 
_pfear, df —[f 25V, df —JSf ®6V, df. 
c a i 
Mit der ersten Variation der Formanderungsarbeit 
OA iif Etwe V),0V,|, dt 
oe lea OV, de tL Eture V,), OV,)|, de, 
die bei Anwendung des Gaufschen Integralsatzes auf das zweite Integral in 
0A = — fff Bere Vly, 0Vy de + [f EXHre Vy 3V, md 
=— fff Breve V,,, OV, dx + ff w ov, df 
R F 
_tibergeht, und mit der Identitat 
— ffeov, df—ff Pov, df — ff B® dV, df 
=— ff woV, df+ ff eov, df + ff ®ov, df + ff ® ov, df 
: $ é : 
-folgt daraus die verschwindende erste Variation zu 
df =0 =— fff (E*#"e VI, + p4) OV, de + 
R 


ad st) OV Be Ys aL) 


— [fo (V,— Vi) of — Jf 62 (V.— V,) af— oe — ¥;) df. 


Thr entnimmt man die Differentialgleichungen 
. E+we Vi, + pt =0 
und die Randbedingungen 


: [oJ = [E4uve VI, ma] — th; (ie saa a 
[Vile =V, 3 ( = 1,2, 3), 
G 


ohne daB an die variierten GréBen Forderungen zu stellen sind. 


5. Hamiltonsches Prinzip. Bezeichnen g die Dichte des Materials und ¢ die Zeit, so gilt bekannt- 
lich fir die kinetische Energie des elastischen Kérpers 


: eatin: 
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1 


a 
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Mit dieser lautet das Hamiltonsche Prinzip in der Kinetik elastischer Korper’ 
I =f (T— U) dt = Extrem. 


[6Y]i = 0, 


G 


[6V,.1=9, 


sche Prinzip in der modifizierten Form 


[dVaJ2 =0, 
G 


[dV,]3 =0. 
G 


[OV Jp = 9, 
Unter Verwendung der potentiellen Energie // (VY; Vo, Vay t, 0, 08, Vis Vo, V3) kann das Hamilton- 
a 
S =f (T—JZ) dt = Extrem. [Ov a =O: 
1 


[dV le = 0 
bzw 


geschrieben werden. Dabei sind auch hier die Bedingungen fiir die variierten Verschiebungen 
[dV,,] a (lee 
2 


(u = 1,2, 3) 
entfallen, und das Variationsprinzip hat die gleiche Form wie in der Kinetik der Massenpunkt- 
systeme. Wegen JJ = U+Z=— (u +- Z) erhalt man die beiden Formulierungen 
2 
S=f(T— U—Z) d& = Extrem. 
1 
* * 
S={(T4+U+42Z)d = Extrem. 
i 


[OV], =0, 


[Vile = 9, 
[0 Vil =0, 


(Eingegangen am 26. Juli 1960.) 


[d Vile => 0 ° 
Anschrift des Verfassers: Professor Dr.-Ing. D. Riidiger, Freiberg (Sachsen), StraBe der Einheit 12. 
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